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Abstract 



In this thesis, we consider the suitability of using the charged cold fluid model 
in the description of ultra-relativistic beams. The method that we have used is the 
following. Firstly, the necessary notions of kinetic theory and differential geometry 
of second order differential equations are explained. Then an averaging procedure 
is applied to a connection associated with the Lorentz force equation. The result of 
this averaging is an affine connection on the space-time manifold. The corresponding 
geodesic equation defines the averaged Lorentz force equation. We prove that for 
ultra-relativistic beams described by narrow distribution functions, the solutions of 
both equations are similar. This fact justifies the replacement of the Lorentz force 
equation by the simpler averaged Lorentz force equation. After this, for each of these 
models we associate the corresponding kinetic model, which are based on the Vlasov 
equation and averaged Vlasov equation respectively. The averaged Vlasov equation 
is simpler than the original Vlasov equation. This fact allows us to prove that the 
differential operation defining the averaged charged cold fluid equation is controlled 
by the diameter of the distribution function, by powers of the energy of the beam 
and by the time of evolution t. We show that the Vlasov equation and the averaged 
Vlasov equation have similar solutions, when the initial conditions are the same. 
Finally, as an application of the averaged Lorentz force equation we re-derive the 
beam dynamics formalism used in accelerator physics from the Jacobi equation of 
the averaged Lorentz force equation. 
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Chapter 1 



Introduction 

1.1 Motivation of the thesis 

Current models of classical electrodynamics of charged point particles contain logical 
inconsistencies that arise when back-reaction effects are considered. For example, 
the standard theory of back-reaction is based on the Lorentz-Dirac equation [1-5]. 
However, it is well known that the Lorentz-Dirac equation is problematic from a 
physical point of view: some of its solutions contain pre-acceleration effects; others 
are run-away solutions. This is the case for a large class of initial conditions. This 
peculiarity of the Lorentz-Dirac equation is due to the fact that it is a third order 
differential equation. 

A possible solution to the problems of the Lorentz-Dirac equation is the theory 
proposed by Landau and Lifshitz in [2] (recently reviewed for instance in reference 
[5]). From the analysis of this question performed in reference [5], one extracts the 
following conclusion: 

The charged point particle description is valid iff the changes in the acceleration 
of the particle occur over time scales longer than the characteristic time parameter 

The parameter q is the charge and m is the mass of the point particle and one 
has assumed units such that the permeability of the vacuum eo is 1 and the speed of 
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light c is also set equal to 1. If the point particle approximation condition holds, the 
Landau-Lifshitz reduction of order procedure can be applied to the Lorentz-Dirac 
equation to obtain a second order differential equation as an approximation, free 
from pathological solutions. In this regime, the Landau-Lifshitz equation can be 
considered an appropriate approximation of the Lorentz-Dirac equation. 

Despite solving the problem of the Lorentz-Dirac equation, there are several rea- 
sons why the solution proposed by Landau and Lifshitz is not completely satisfac- 
tory: 

1. The order reduction procedure is an ad-hoc procedure (although consistent 
with the point particle picture). 

2. The Landau-Lifshitz equation is the leading order term approximation of the 
Lorentz-Dirac equation. Therefore, it is not a fundamental equation. 

3. The characteristic time to is proportional to ^q. Therefore, let us consider a 
physical system with a large number of identical charged particles performing 
a collective motion. The prototype example is the motion of a bunch of parti- 
cles in an accelerator machine. It can happen that the behavior of the system 
is coherent and that one has to read the factor q as the total charge of the 
bunch and m as the total mass of the whole bunch. Under these conditions, 
it is natural to consider that the factor -2- remains the same as for an individ- 

m 

ual charged point particle, but q increases proportionally with the number of 
particles. Then for intense beams of particles, the point charge approximation 
and the reduction of order procedure will break down. 

Even if the first two points can be covered under the interpretation of classical 
electrodynamics as the limit of the fundamental quantum electrodynamics, the third 
point has relevance for us. The energies and luminosity achieved in modern particle 
accelerators can push to the validity of present models of electrodynamics its limits. 
This is basically because one is dealing with bunches containing a large number of 
charged particles, which can reach 10 9 — 10 11 particles per bunch, moving together 
in a small phase-space domain (all the particles are concentrated around a center of 
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mass, in position and velocity). 

Since the possible effect discussed in point 3 is additive, for modeling systems like 
those bunches of particles, one needs an alternative description to Lorentz-Dirac and 
Landau-Lifshitz models. 

In this context, fluid models have been used to study the dynamics of ultra- 
relativistic beams of charged particles. One of these models is the proposal contained 
in [6]. In that work, it was shown how to do an asymptotic analysis of the charged 
cold fluid model. The main claim in [6] was that the model proposed provides a self- 
consistent description of the fully coupled dynamics of a bunch of particles with the 
electromagnetic field. The reason for this is the smoothness properties of the fields, 
compared with the discrete and singular character of the point particle description 
behind the Lorentz-Dirac equation. 

However, the use of the charged cold fluid model was not justified in [6]. This 
justification is necessary, because of the discrete nature of a bunch of particles. 
Therefore, prior to the use of this model, one has to address the following question: 

When is it a good approximation, in the regime of ultra-relativistic dynamics, to 
describe the interaction of a large number of charged point particles with the total 
(external and associated) electromagnetic field by a charged cold fluid model? 

A simplified version, is the following question: 

When is it a good approximation, in the regime of ultra-relativistic dynamics, to 
describe the interaction of a large number of charged point particles with an external 
electromagnetic field by a charged cold fluid model? 

In the present thesis we address this second question. In particular, we present 
an averaged description of the collection of charged point particles, defining a mean 
velocity vector field. The averaging operation is interpreted from a kinetic theory 
point of view, introducing the one-particle distribution function as a solution of 
the Vlasov equation [7,8] and the associated averaged Vlasov equation, that we will 
introduce later. Our final result in this direction is contained in theorem 5.3.7, which 
can be stated in words in the following way: 
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For narrow distribution functions and in the ultra-relativistic regime, one can 
usehe t charged cold fluid model as a good approximation to the Vlasov model, in the 
dynamical description of ultra-relativistic bunches of charged particles. The error of 
the approximation is of the same order as the area of the support of the distribution 
function in velocity space. 

What this result means is that, in this regime, if the Vlasov equation holds, the 
differential equation defining the charged cold fluid equation holds approximately. 
Also, we note that a precise statement is involved, requiring some technical assump- 
tions which we will discuss in the appropriate place. 

In order to achieve the above results, a connection associated with the Lorentz 
force connection will be introduced (the Lorentz connection). We introduce an 
averaged version of this connection (the averaged Lorentz connection). The main 
advantage of this technique is that the averaged connection is simpler than the 
original one. This allows us to perform calculations whose results are not easily 
obtainable using to do in any other way. 

1.1.1 Other results of the thesis 

Another application of the theory of the averaged Lorentz model is the following. 
After introducing the Jacobi equation of an affme connection, we discuss the Jacobi 
equation associated with the averaged Lorentz connection. Then we prove that the 
linear dynamics used in accelerator physics [9-11] is an approximation to the Jacobi 
equation of the averaged Lorentz connection. Based on this interpretation, we define 
a notion of reference trajectory in beam dynamics in terms of observable quantities. 
Also, using the averaged Lorentz equation, we provide observable consequences of 
the collective nature of the bunch of particles. 

We have also considered the question of how the gauge invariance principle affects 
the interpretation of the Lorentz force equation as an Euler-Lagrange equation of a 
functional. This led us to a precise definition of semi-Randers space. 

Finally, we should mention that during the analysis of the main problem consid- 
ered in this thesis and its mathematical formalization we found a generalization of 
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the notion of connection in differential geometry. We call this object almost (pro- 
jective) connection. This is described in the appendix. 

1.2 Structure of the thesis 

In this chapter we introduce some notation and conventions that we will follow 
through this thesis. 

In chapter 2, an introduction to relativistic kinetic theory is provided, following 
reference [7] . Then we define the charged cold fluid model and consider the asymp- 
totic method developed in [6]. We will state the main problem considered in this 
thesis and give a short out-line of the strategy to solve it. 

In chapter 3, we introduce the theory of non-linear connections defined by a second 
order differential system [12,13]. We also introduce the formalism and the notion 
of averaged connection, following the method already contained in ref. [22]. In 
particular we define the average of linear connections acting on sections of some 
relevant bundles (the pull-back bundles 7r*T( p ' 9 )M). The data that we need to 
determine these connections is a system of second order differential equations called 
semi-spray. Those connections are obtained basically from the structure of the 
corresponding differential equations. 

The original content of the Thesis constitutes chapters 4, 5, 6, 7 and 
appendix 4. 

In chapter 4, the notion of semi-Randers space is introduced and discussed as a 
geometric description of the interaction of charged point particles with an external 
electromagnetic field [14, 15]. Then the Lorentz connection is obtained. Following 
the theory described in chapter 3, the corresponding averaged Lorentz connection is 
determined. It turns out that, if the dynamics happen in the ultra-relativistic limit 
and the support of the probability distribution function / is narrow (in a sense to 
be specified), the solutions of the Lorentz force equation can be approximated by 
the solutions of the averaged Lorentz force equation. We give an estimate of the 
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approximation as a function of the time of evolution, the energy of the system and 
the diameter of the distribution. 

In chapter 5 it is proved that under the same assumptions as in theorem 4-6-6, the 
relativistic charged cold fluid model can be obtained as an approximation from a 
kinetic model. The method that we follow to obtain this conclusion is the following. 
First, we introduce the averaged Vlasov equation and compare it with the original 
Vlasov equation. In particular we prove that, for the same initial conditions, both 
models have similar solutions in the ultra-relativistic regime when the distributions 
functions are narrow. After this, we use the averaged Vlasov model to give a bound 
on the acceleration of the main velocity vector field of the averaged Vlasov model. 
This bound is given in terms of the diameter of the distribution function, the energy 
and the time evolution. Then we prove that the mean velocity field associated with 
the solution of the Vlasov equation is similar to the mean velocity field obtained 
from the solution of the averaged Vlasov equation. This fact finishes the proof of 
theorem 5.3.7, which is the answer to the main problem considered in this thesis. 

In chapter 6 we use the Jacobi equation of the averaged connection to provide 
a geometric formulation of the transverse and longitudinal linear beam dynamics. 
Particular examples illustrate the general formalism. We obtain from the Jacobi 
equation of the averaged Lorentz connection the equations of motion of the trans- 
verse dynamics in magnetic dipole and quadrupole fields. In a similar way, the lon- 
gitudinal dynamics in a constant and alternating electric field are obtained from the 
Jacobi equation. Although these are known examples, they illustrate the usefulness 
of the Jacobi equation of the averaged connection in beam dynamics. Corrections 
to the ordinary dynamics coming from collective effects are considered. As an ap- 
plication of the formalism we can provide a definition of reference trajectory that by 
construction is given in terms of observable quantities. 

In chapter 7 we discuss some of the results presented in this thesis as well as 
perspectives for further developments. 

The present thesis work has produced the following articles and pre-prints [14], 
[15] and [16]: 
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1. R. Gallego Torrome, On the Notion of Semi-Randers Spaces, arXiv:0906.1940 , 

2. R. Gallego Torrome, Geometric Formulation of the Classical Dynamics of 
Charged Particles in a External Electromagnetic Field, arXiv:0905.2060, sub- 
mitted. 

3. R. Gallego Torrome, Fluid Models from Kinetic Theory using Geometric Av- 
eraging, arXiv:0912.2767, submitted. 

4. R. Gallego Torrome, Averaged Lorentz Dynamics and an Application in Plasma 
Dynamics, arXiv:0912.0183 , accepted to publish in the Proceedings of the 
XVIII Fall Meeting in Geometry and Physics, American Physics Society. 

1.3 General conventions used in the thesis 

For the main physical applications in this thesis, the space-time structure will be a 
flat four dimensional manifold endowed with a Lorentzian metric n with signature 
(+, —,—,—). However, some results and techniques are valid for arbitrary dimension, 
signature and curvature. In these cases, it is explicitly stated. Sometimes we will 
require, to simplify the calculations, that the metric n is flat. When this is the 
case, it will be indicated. In all cases we assume that the space-time manifold M is 
time-oriented. 

Einstein summation convention is considered for any identical and repeated co- 
variant and contravariant indices, if the contrary is not stated. All Latin indices 
run from to n — 1, where n is the dimension of the space-time manifold. Vector 
notation is used for the spatial components (with respect to a given frame) of a 
vector. Indices are lowered using the metric ijij and raised using the inverse metric 
rfi , unless anything else is stated. The exterior product and the exterior derivative 
are normalized as in reference [17]. 

We have adopted the following convention for the physical parameters and con- 
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stants appearing in the models, 



q = 1, m = 1, e = 1, fio = 1, c 



where eo are the dielectric and magnetic permeability constants of the vacuum; 
m is the mass and q the charge of the species of particles that we are considering. 

We also use the following convention [3, pg 618]: 



D = E + P, H = B — M. 



(1.3.1) 



P is the polarization vector and M is the magnetization of the medium. Since we 
are considering that the bunch of particle propagates in the vacuum, we have that 
P = 0, M = 0. Hence, one has E = D, H = B. 

The Maxwell equations are 



V-E = p, VxB=(/+^), 



(1.3.2) 



dB_ 

dt ' 



V-B = 0, VxE = - — . (1.3.3) 
The Maxwell equations can also be written in a covariant form in the following way: 



d l F jk + dkFji + dj F ik = 0, ''V.F 1 j = r) k3 J k 



(1.3.4) 



where Wj is the covariant derivative associated to the Levi-Civita connection along 
the direction e^. 

The electromagnetic tensor is described by a 2-form, that in a local frame deter- 
mines the following matrix: 



( E x (x) E 2 (x) E 3 (x) \ 

-E^x) -B 3 {x) B 2 (x) 

-E 2 (x) B 3 {x) -Bxix) 

\ -E 3 (x) -B 2 (x) B x {x) / 
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The Lorentz force is written as 



^ / a\ da 

F:=q(E + vxB), v= — , 



m- 



IF' 



(1.3.5) 



The parameter r is the proper-time along a associated with the metric r\. In covari- 
ant formalism, the Electromagnetic field is described by a 2-form F = F^- dx l A dxK 
The Lorentz force equation for a mass m = 1 and a charge q = — 1 is 



dr 2 



d7 



(1.3.6) 



There are several categories of metric structures where we will work. The most 
general is semi-Riemannian category [47]. The results stated in this category refer 
to structures {rj(x)} which are non-degenerate, symmetric bilinear forms for each 
fixed x € M, and are smoothly defined on the n-dimensional manifold M. 

The second most general category refers to Lorentzian manifolds. In this case 
the metrics {t](x)} are billinear, symmetric forms with signature (+,—,—,—) for 
fixed x £ M and are smoothly defined on the manifold is four dimensional space- 
time because physical reasons. Those results be defined may also be defined for 
n-dimensional space-times with signature (+, — , — ). 

The third category of geometry is Minkowski geometry. In this case the metric 
is the Minkowski metric but the manifold considered is a domain inside the space 
R 4 . In this category it makes sense to speak of global inertial frames: a global 
inertial frame (eo,ei,e2,es) on M such that the metric r\ has the following metric 
components: 



(Vij) = 



( 1 o \ 
0-10 

0-10 

-1 



The contraction operation is defined on the following way: given a tangent vector 
W G T X M, there is an homomorphism on the space of covariant tensors over x 
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denoted by: 

l w : T^M — > T^'f-^M 

given by 

tw r(Xi,...,x I ,_i):=r(w;xi,...,x p ), vi^t.m, feT.M. 

A similar definition applies pointwise to sections of J 7 (M)-multilinear maps of vector 
fields contracted with a given vector field given by 

L W T(X 1 ,...,X p _ 1 )(x):=T(x)(W,X 1 ,...,X p ), VJQ G T(TM), W £ r(TM). 

When we write down the results, we try to formalize in the largest category 
possible. Generally speaking, results from chapter 3 fall into the category of semi- 
Riemannian metrics (indeed, some of them are even more generic that for metric 
structures) . Results in chapter 4 fall into this category of semi-Riemannian category, 
except for the main comparison results, where we explicitly use the flatness property 
of the metric r\ in some of the calculations. 

In chapter 5, the results depend on the results of chapter 4- Therefore, although 
some of them are formulated for semi-Riemannian manifolds, the main results are 
formulated for compact domains of the Minkowski space. 

In chapter 6, the main results are formulated for Minkowski space, since we have 
in mind to apply the geometric formalism to describe the behavior of beams of 
particles in accelerators, where gravitational effects are usually neglected. 

In chapter 7, we point out the general conclusions of this thesis as well as open 
problems proposed. 
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Chapter 2 

Fluid and kinetic models for 
ultra-relativistic beams 

In this chapter we consider some basic notions that we will use later. In the same 
way, we introduce additional notation. There is also a short introduction to fluid 
models, kinetic models and to the asymptotic model described in [6]. 

2.1 Basic relativistic kinetic theory 

In this section we review some elementary notions of the covariant kinetic theory 
which are relevant for our work. We mainly follow the notation of reference [7]. We 
will consider collision-less processes and detailed balance processes. 

2.1.1 Intrinsic covariant formalism for relativistic 
kinetic theory 

In this thesis, the kinetic models are based on the following general assumptions: 

1. The space-time manifold M is 4-dimensional and it is endowed with a Lorentzian 
metric rj. The signature of the metric is (1, —1, —1, —1), and the space-time is 
time orientable [7] and oriented. In general, the metric rj is not flat. However 
we will use the flatness condition on the metric i] in some of the calculations. 
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The Lorentzian metric r\ has an associated Levi-Civita connection W. Also, 
r\ determines the Hodge star operator 



★ : TA P M — ► rA 4_p M 
0Jn... ip e h A ... A e 1 " ^ uj h ... ip e^ +1 A ... A e u -" e* 1 -** ip+1 ...U- p - 

TA P M := {uj : M — > A P M} is the set of smooth sections of the vector bundle 
A P M ->■ M, with A P M the bundle of smooth p-forms over M; e ll - lp i p ...i 4 _ p 
is the total skew-symmetric symbol, where the indices are raised using the 
Lorentzian metric r] and with €0123 = 1- ^i 1 ...i p e tl A ... A e tp is an arbitrary p- 
form expressed in a dual basis {e°, e 3 } of an orthonormal basis {eo, 63}. 

The dual of a vector field is the 1-form defined pointwise by the relation 

(Vf(x)=r ] (V(x),-). 
Similarly, the dual of a 1-form is a vector field defined pointwise by the relation 

J(x)=r ] - 1 (u,.), 
where rj~ l is the bilinear form 

v -i . x*M x T*M — ► R 

2. The electromagnetic field is encoded in the 2-form F, while the excitation field 
is encoded in the 2-form G and the current density is a 3- form J, all living 
on M. They satisfy Maxwell's equations, which can be written in terms of 
differential forms as 

dF = 0, d*G = J; (2.1.1) 

d : /\ p M — > A P+1 M is the exterior derivative operator acting on forms. This 
is a coordinate free form of the equations (1.3.4). 
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3. The relation between F and G is given by the constitutive relations. We 
assume that these relations are linear and in particular we put F = G, since 
the electromagnetic medium that we are considering is the vacuum. 

4. The matter content of the models consists of a collection of identical charged 
point particles. The trajectory a(s) of each particle follows the Lorentz force 
equation 

f = F M r W , ✓ = (2.1.2) 

where r is the proper-time associated with the trajectory. The parameter r is 
such that V (*p,*p) = l. 

5. The support of the one-particle distribution function f(x, y) is in the 7-dimensional 
unit hyperboloid bundle, 

£ := {(x, y), x 6 M, y G T,M, V (y, y) = 1, y° > }. (2.1.3) 

A global coordinate system on ~E is (x a ,y l ), a = 0,1,2, 3, i = 1,2,3, in- 
duced from any natural coordinate system on TM. In the unit hyperboloid, 
y° is given as a function of y 1 , y 2 , y 3 and x a . The manifold := {y £ 
T X M, rj(y,y) = 1, y° > 0} is called the unit hyperboloid over x. 

6. There is defined a volume form on the unit hyperboloid bundle S. This volume 
form is obtained in terms of the metric rj. On the tangent space TM there is 
a volume 8-form 

^\detri\ dy° A • • -dy 3 A dx° A ■ ■ ■ A dx 3 , d 4 (x) = dx° A • • -dx 3 , 

with y/rj the determinant of the matrix associated to the metric in a given 
coordinate system. The isometric embedding e : ~E ■— > M induces a volume 
form on the manifold We denote this volume form by dvol(x, y)Ad 4 x. Since 
the space-time manifold M is 4-dimensional, the volume form dvol(x,y) Ad x 
is a 7-form. 

The volume form dvol(x,y) on ~E X is obtained contraction of dvol(x,y) A d A x 
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on the orthogonal frame {eo, 63}: d 4 x(eo, e%) = 1. 

7. The Liouville vector field L \ is tangent to the hyperboloid Using the 
conventions of section 1.3, the Liouville vector field L \ can be written using 
local coordinates as 

L X = y% + (F i j yi - r jk yiy k )— t , i,j,k = 0, 1,2,3. (2.1.4) 

Remark Note that we have adopted the extrinsic formalism, where y° is con- 
sidered an independent coordinate. Later we will explain the relation between 
the intrinsic and extrinsic formalism, and that they are equivalent for our 
purposes. 

8. The one-particle distribution function f(x,y) is defined over £ and satisfies 
the equation 

L x(f) = 0- (2.1.5) 

Equation (2.1.5) corresponds to the Vlasov equation in plasma physics and 
kinetic theory. The one-particle distribution function f(x,y) is introduced as 
the probability density of finding a particle at the point x € M with velocity 
vector y G T X M. [7]. This interpretation was supported in [7] using balance 
arguments and assuming that f(x,y) is continuous. We will also assume ad- 
ditional smoothness and regularity conditions for f(x,y). 

2.1.2 Extrinsic formulation of the kinetic model 

There is an alternative description of a kinetic model to the intrinsic one. In this 
alternative description the calculations are performed on the whole tensor bundle 
TM and then the results are restricted to the unit hyperboloid. One uses the 
constraint rj(y, y) = 1 when it is necessary. Note that the action of L \ is on the 
ring of smooth functions of the hyperboloid ~E, since L \ is a tangent vector to the 
unit hyperboloid ~E. This follows from the fact that L x|s ■ {viUiU)) = an d that 
the function rj(y, y) = r}ij(x)y % yi generates a foliation of TM. A formal proof of this 
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fact can be found for instance in [18]. 

We will also define later an averaged Liouville vector field < L \ >• This vector 
field lives on M rather than on the tangent bundle TM. However, it defines a second 
order differential equation and therefore a Liouville equation. The flow of < L \ > 
does not preserve the function r)(y,y)(x) := r)ij(x)y l y^. Indeed, it is not guaranteed 
that the flow will preserve a structure. Therefore for the study of these kind of flows, 
it is more convenient to adopt the external formalism. 

Using the volume form dvol(x,y) one can obtain the velocity moments of the 
distribution f(x,y). Therefore one can define moments of the distribution function, 
which are the expectation values of polynomials on y. With these moments, one can 
define the mean velocity field, the covariant kinetic energy-momentum tensor and 
the covariant energy-momentum flux tensor: 

^ - ^w*,*) L Mnx ' v)iai{ '' v) ' (2X8) 

The balance equation for the number of particles implies the relations [7] 

"V^(x) = F'j ;V* + - r — l —— -[ y lL X (f)dvol(x,y). (2.1.10) 

}sJ{x,y)dvol{x 1 y) J Sx 

Since / follows the Liouville equation (2.1.5), one obtains 

"ViV^x) = 0, "^^(i) = ¥ l j V j . 
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2.2 Relativist ic charged cold fluid model 



We introduce some geometric and physical objects that we need in the description 
of the asymptotic expansion of the relativistic cold fluid model proposed in [6] . The 
electromagnetic field is encoded in the 2-form F, which is a solution of the Maxwell 
equations (2.1.1). The external electromagnetic field F is created by the external 
current density J such that in the space time regions that we will consider, one has 
that J(x) = 0. The current density J describes a system of charged point particles 
which also contributes to the total electromagnetic field. The whole dynamics is non- 
linear and one needs additional information to completely determine the dynamics. 
There are two additional pieces of information: 

1. One has to postulate the dynamic equation for the current density J . Exam- 
ples for these joint dynamics are the Maxwell-Lorentz system, Maxwell- Vlasov, 
Klimontovich-MaxwelPs system [8, section 2.5]) and Maxwell-Lorentz-Dirac 
system [5]. 

2. In order to completely determine the system, constitutive relations between F 
and G are needed. We assume that the constitutive relations are G = eoF. 
We have adopted units such that eo = 1. 

In flat regions, the metric rj admits a set of translational Killing vectors {iQ, i = 
0,1,2,3}, C Ki T] = 0. 

Using differential forms, one can write conservation laws in a geometric way. For 
any vector field W on M there is an associated drive 3-form [17]: 

Tw = 2(twFA*F-t w *FAF). 

For the case of Killing vector fields, the exterior derivative of t^J 1 is 

dr^P = -l w Y A J, 
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where lw~F is the contraction of the vector field W with the 2-form F. In the region 
outside of the sources J = 0, 

dr^ 1 = 0. (2.2.1) 

In the presence of matter, equation (2.2.1) has to be generalized. For example, let 
us consider a model for matter described by a time- like vector field V. Then for a 
dust, the stress-energy tensor is 

T(x) = Mv\x) ®V\x). (2.2.2) 

M is a regular scalar density field and the velocity field is normalized, r](V, V) = 1. 
The current density J is proportional to the velocity field: 

J = N*(V) b . (2.2.3) 

and then dj = 0. Combined with the assumption of the total momentum conser- 
vation dijyj 1 + *iwT) = implies the field equation of motion for the fluid: 

r 'V v V{x) = (i V Ff (2.2.4) 

as a balance equation [17, pg 242-243], [19]. 

The dynamics of the relativistic charged cold fluid model is described by the 
following coupled system of differential equations, 

dF = 0, d*F = -p*V\ r iV v V{x) = {i V ¥)\ t](V,V) = 1. (2.2.5) 

Although this is a complete model, in this thesis we will work with external electro- 
magnetic fields. In this context, the mathematical and physical analysis are highly 
simplified. 
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2.3 Asymptotic expansion of the relativistic charged cold 
fluid model 

Let us consider the following 1-parameter family of differential forms and vector 
fields: 

+00 +00 +00 

n=— 1 n=l n=—l 

where 

T4 G TTM, p n G r A M, F„ G T A 2 M (2.3.2) 

and e is a small parameter. Substituting these expansions in equation (2.2.5) and 
equating terms of equal in e, one obtains enough conditions to determine the fields 
(2.3.2) inductively [6]. For instance, the leading order terms are the vector field V-\ 
and the 2-form F_i such that: 

"Vy.^-i = tv-iF-i, dF_ 1= 0, d*F = 0, 7 7 (F_ 1 ,F_ 1 ) = 0. 

Given initial data for V-\ and F_i, these equations are compatible. Note that they 
describe a charged mass-less fluid interacting with an external electromagnetic fluid. 

In general the equations for the fields appearing in the expansion (2.3.1) are 
obtained from the equations (2.2.5). The procedure for the higher orders is as 
follows: 

1. Consider a given electromagnetic field F_i, solution of the differential equa- 
tions 

dF_i = 0, cUF_i=0 (2.3.3) 

for some initial value of F_i on a space-like hypersurface. Physically F_i is 
interpreted as the external electromagnetic field. 

2. Then one has to solve the equation 

"Vy.^.i = (tv^F.i)" (2.3.4) 
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subject to the condition 

r,(V- 1 ,V- 1 ) = 

and for given initial data in the space-like hypersurface t = to. This is possible 
because equation (2.3.4) can be re-written as an ordinary differential equation 
and one can apply standard results (see for instance the appendix). 

3. Then one solves the equation for pi, which is 

d*(p 1 Vi 1 )=0, (2.3.5) 

for given initial values in a space-like hypersurface. 

4. The 2-form Fo is a solution to the Maxwell equation, that can be written as 

dF = 0, d*F = -*/oi Vli, (2.3.6) 

where one has to specify the initial values on a space-like hypersurface. 

5. Vo is the solution of the equation 

"Vk.^o + = (iv^Fo + ty F_i)" (2.3.7) 

subject to the requirement that r/(V_i,Vb) = and for fixed initial values of 
the vector field Vq in a space-like hypersurface. 

6. The density p 2 is defined as the solution of 

d*{p2V\) + d*(p 2 V\) =0, (2.3.8) 

again after given the initial values of p2 on a space-like hypersurface. 

7. The equation 

d*Fi = -■kp 2 Vt 1 -*PiVq (2.3.9) 
can be solved for Fi, once the initial values for Fi are specified. 
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8. V\ is a solution of the equation 

"Vy.^i + v Vv V-i + "V yi y_i = (ty^Fo + l Vo F-i + ty^-if- (2.3.10) 

We need to specify the initial values on a space-like hypersurface. 

Through a generalization of this procedure, the fields (2.3.2) can be solved order 
by order in e. The only non-linear differential equation to be solved is for V-\. 
Indeed it can be written as an ordinary second order differential equation for the 
integral curves of V-\. These properties make it easier to solve both the analytical 
and numerical treatment of the problem than to solve the original equations (2.2.5), 
which are a system of non-linear and coupled partial differential equations. 

The vector field V_i has a difficult physical interpretation, because it corresponds 
to a charged cold fluid composed of mass-less particles and at the same time inter- 
acting with an external electromagnetic field. There is no known classical physical 
system (in vacuum) with such characteristics (in quantum physics, the low energy 
limit of graphene admits states which are mass-less and interact with the electro- 
magnetic field [56]). 

2.4 Statement of the main problem considered in this 
thesis and out-line of the strategy to solve it 

It was claimed that the model introduced in [6] is able to provide a consistent 
treatment of the back reaction and self-force problems that appear in classical elec- 
trodynamics, for some situations which are of practical interest like ultra-relativistic 
plasmas. This claim is based on the assumption that the charged cold fluid model 
is an acceptable description of the dynamics of bunches of particles in the ultra- 
relativistic regime. 

On the other hand, fluid models have been used intensively in the description 
of the dynamics of plasmas [8, 37-39]. However, these usual models are based on 
assumptions on the moments of the distribution function, which are difficult to check 
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in experimental conditions. 

These motivate the main problem considered in this thesis: 

7s it mathematical justified to use the charged cold fluid model (2.2.5) under the 
conditions present in the currently used particle accelerators? 

We will estimate the value of the differential operators appearing in the equation 
v V v V(x) = (iyF)H, where V is the mean field (2.1.6) for a given distribution func- 
tion /. Then we will show in chapter 5 that the differential expression for the charged 
cold fluid model equation is bounded and controlled by powers of the diameter a of 
the distribution function f(x, y), powers of the energy}] of the system and powers of 
the coordinate time evolution t. The relation is such that for narrow distribution 
functions and in the ultra-relativistic regime, the charged cold fluid model is a good 
approximation of the kinetic model. 

The strategy that we will follow is the following. Since the fluid model V{x) is an 
approximate description of the system, we interpret V(x) as an averaged quantity. 
On the other hand, given a dynamical system, we can associate a non-linear connec- 
tion. That connection can be averaged, using a distribution function. If in addition, 
the difference between the original connection and the averaged connection is small, 
one can substitute the original one by the averaged connection in the description of 
the dynamics. 

Any dynamical system described by a connection has an associated kinetic model. 
In particular, this is true for the Lorentz force equation and the averaged Lorentz 
force equation, that we will define. Since the Lorentz connection is similar to the 
averaged Lorentz connection in a sense that we will explain later, one can also 
substitute the associated kinetic models. 

Working with the averaged model has technical advantages. In particular one can 
give estimates of the value of some differential operators which appear in the fluid 
models. This will be done for the charged cold fluid model. 



The notion of energy of a bunch of particles that we will use in those bounds is not trivial and 
will be introduced in chapter 4- 
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Chapter 3 

The averaged connection 



In this chapter we introduce the notions of non-linear connections and the associated 
averaged connections, before applying the method to the connection associated with 
the Lorentz force equation in the next chapter. The construction is adapted from 
reference [22]. This chapter explains the mathematical theory that we will use in 

chapters 4, 5 and 6. 

3.1 Non-Linear connection associated with a second or- 
der differential equation 

3.1.1 Second order differential equations and the associated non- 
linear Berwald-type connection 

Let M be an n-dimensional smooth manifold. A natural coordinate system on the 
tangent bundle ir : TM — > M is constructed in the following way. Let (x, U) 
be a local coordinate system on M , where U C M is an open sub-set of M and 
x : U — > R™ a local coordinate system. An arbitrary tangent vector at the point 
p 6 U is of the form X p := X = X k -^\ p . The local coordinates associated with 
the tangent vector X p <E T X M. C TM are (x k ,y k ). We will identify the point 
x £ M with its coordinates, by notational convenience. N is a sub-bundle of the 
tangent bundle TM. From the imbedding e : N > TM, e(N) acquires the induced 
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differential structure from TM; e(N) is denoted by N. 

We recall the following notion of connection [20, pg 314]. Let 7r : N — > M be a 
bundle over M and consider the differential function dn : TN — > TM. Then the 
vertical bundle is V = ker(dir) C TN. 

Definition 3.1.1 A connection in the sense of Ehresmann is a distribution % C 
TN such that 

1. There is a decomposition at each point u G N ; T U N = % u © V u . 

2. The horizontal lift exists for any curve t \-¥ a(t) G M, t\ < t < ti and is 
defined for each £ G T M M and u G 7r~ 1 (x). 

Let us consider a set of n second order differential equations, with n the dimension 
of M. The solutions are parameterized curves on M. Assume that the system of 
differential equations describes the flow of a vector field G x G TTN. In particular, 
the system of differential equations has the following form 

i = !,...,„. (3.1.1) 

This system of differential equations is equivalent to the following system of first 
order differential equations on N, 

d \ K ' (3.1.2) 
% = y\ i = l,..., n. 

The coefficients G' l (x,y) are called spray coefficients if they are homogeneous func- 
tions of degree one on the coordinate y; in the general case where they are not 
homogeneous those coefficients are called semi-spray coefficients. G % {x, y) transform 
under a change of natural local coordinates on N, induced from changes of coor- 
dinates on M, in such a way that the system of differential equations (3.1.1) is 
covariant. Explicitly, if the change of local natural coordinates on the manifold N 



32 



IS 



x 1 = x l (x), 

I •'/• = fry/'- 



then the associated co-frame transforms as 



~i — d 2 x l „,k 



dy 



dx k dxi 



yk dx j + grfyi. 



The induced transformation in the associated system of differential equations is 



where the coefficients G l (x,y) are 

& ^y) = E ^)y ] ^)y k ^L - E (f^M^)- 



The vertical distribution V admits a local holonomic basis given by 



(3.1.3) 



Using these spray coefficients it is possible to define a horizontal n-dimensional 
distribution of the fiber bundle TN — > N. The basis for the distribution is 



d dG i d 



5 5 5 

dx 1 ' ' 5x n ' 5x k dx k dy k dy' 1 



i,j = l,...,n. 



(3.1.4) 



It generates a supplementary distribution to the vertical distribution. The non-linear 
connection coefficients N l k{x,y) are defined by the relation 



G\x,y) :=y k N\(x,y). 



For a spray, the connection coefficients of the non-linear connection are: 
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Since the spray coefficients G % are transformed under the a change in natural co- 
ordinates in a well defined way, the non-linear connection coefficients N l j(x,y) are 
also transformed in a characteristic form [24, 34], 

Given a spray G l (x,y), we define the connection coefficients such that the only 
non-zero coefficients correspond to the covariant derivative of horizontal sections of 
TTN along horizontal sections of TN and such that they are given by the Hessian 
of the spray: 

Jk[ ,V) 2 dyidy k ■ 
All the other coefficients are zero. This type of connection resembles the so-called 
Berwald connection used in Finsler geometry [45]. From the point of view of the 
geometry of sprays, it is a natural connection. Note that, while P jk(x,y) can be 
associated with a linear connection on TN, the connection coefficients N* j(x,y) 
cannot (this is why they are called non-linear connection coefficients); N l j(x,y) 
determines a connection which is non-linear in the direction of the derivation. 

Given the non-linear connection, one can define the horizontal lift of the tangent 
vectors; the horizontal lift of X = X % di G T X .M to the space T U N is defined 
by h{X) = X 1 ^. This lift is defined here using local coordinates. However, an 
intrinsic definition can be found in [13]. After introducing this lift, one can define 
the horizontal lift of vector fields and tensor fields of the corresponding bundles. 

3.1.2 The pull-back bundle 

Let us consider the product N x TM and the canonical projections 

7Ti : 7r*TM — »• N, (u, — »• u, 
vr 2 : 7r*TM — »• TM, (u, f ) — »• f . 
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The pull-back bundle 7r*TM — > N of the bundle TM is the minimal sub-bundle of 
the cartesian product N x TM such that the following equivalence relation holds: 
for every u G N and (u,£) G ^(u), (u,g) G 7r*TM iff it o ir 2 (u,£,) = n(u); The 
pull-back bundle 7r*TM — > N is such that the following diagram commutes, 

ttTM^TM 

7T1 7T 

N M. 

7r*TM — > N is a real vector bundle with fibers diffeomorphic to T^M. For instance, 
let {ei,i = 0, ...,n — 1} be a local frame for the sections of the tangent bundle 
TM. Then {ir*ei,i = 0, ...,n — 1} is a local frame for the sections of the pullback 
bundle 7r*TM. Let £ l (x, y)^* x y \ e i( x ) be an arbitrary element in the fiber over 
(x,y) G N; the element ^ xy ^ e i{ x ) is the unique element in 7r*TM such that (it o 

Tl) eiW) = (^2 7r )( 7r (x, 2/ ) e i( x )) and that ° ("'(x.y) 6 *^)) = e «( X )- 

Another way to visualize this pull-back bundle is the following. Let us consider 
the bundle it : N — > M and a fiber 7r _1 (a;) C N. On each point u G 7r _1 (x) we 
attach a copy of the vector space T X M.. This assignment is done by the definition 
of 7r* on a local frame: it* : {ei(x), e n (x)} — > {7r*| u ei(x), 

if* |w^n(*^)} and 

taking linear combinations of the elements of this local frame. When we consider 
sections of the bundle T7r*TM these linear combinations are u-dependent, instead 
of x-dependent. 

Similarly, other pull-back bundles can be constructed from other tensor bundles 
over M, for instance vr*T*M — ► N and v^T^'^M — ► N, with N C TM a 
sub-bundle, T*M the vector bundle of 1-form over M and T^M the bundle of 
(p, q) -tensors over M. 

Given a non-linear connection on the bundle TN — > N, there are several related 
linear connections on the pull-back bundle 7r*TM — > N. 

Let x a semi-spray defined on N. We stipulate the following connection on 7r*TM, 
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defined by the conditions 



V^vr*Z := x T(x, y y jk Z k n*e u V v tt*Z := 0, V G V. (3.1.5) 

Here {7r*ej, i = 0, n — 1 } is a local frame for sections r(7r*TM). This connection 
can be generalized to general tensor bundles over M. 

3.2 The average operator associated with a family of 
automorphisms 

3.2.1 Average of a family of automorphisms 

The averaged connection was introduced in the context of positive definite Finsler 
geometry in [22]. However, in this thesis we need to formulate the theory for arbi- 
trary linear connections on the bundle 7r*TM — > N, where N — > M is a sub-bundle 
of the tangent bundle TM — > M. 

Let tt*, 7Ti, 7T2 be the canonical projections of the pull-back bundle ^T^^M — > N, 
T^'^M being the tensor bundle of type (p,q) over M, tt*T(p>«)M the fiber over 
u G N of 7r*T (p ' 9) M, Tl P ' 9) M the tensor space over x G M S x a generic element 
of t£ p,,?) M and S u is the evaluation of the section S G T(tt*T^M) at the point 
u G N. 

For each tensor S z G Ti P '^M and v G 7r _1 (z), z G U C M the following isomor- 
phisms are defined: 

tt 2 \ v : <T^M — > T^M, S v ^ 

< : Tf^M — > <T^')M, S 2 m- 
To define the averaging operation we need two type of structures: 
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1. A family of non- intersecting, oriented sub-manifolds 

Nu := [J N x , N s C T,M. 

2. A measure at each point x G M, which is an element f(x,y)uj x (y) G /\ m ~N x , 
where m is the dimension of and f x : N x — > [0, oo], f x := /(x, •) is 
required to have compact support on N x . 

Consider a family of endomorphisms, {A w : 7r*TM — > 7r*TM, u G 7r _1 (x)}. Let us 
consider the integral operations 

0K2\uAiK) ■ S(x) := / (vr 2 | u yl u 7r* S(x,u)) f(x,u) u x (u), 

The volume function is defined as 

x i-)- vol(N x ) := / co x (u)f(x,u). 
Jn x 

Definition 3.2.1 Consider a family of endomorphisms, 

{A w : tt*TM — ► tt*TM, w G tt _1 (x)}. 
TTie average endomorphism of this family is the endomorphism: 

< A > x : T X M — ► T.M 
^ — TT^TT / k 2 \ u A u it* u ) 

uo/(e x ) yj^ J 

u G 7T _1 (x), S x G r x M. 

We denote the averaged endomorphisms by symbols between brackets. 

Remark. There is a similar notion which applies to families of homomorphisms, 
instead of endomorphisms between different vector bundles. 

The averaging operation has the following effect. Let us consider an arbitrary 
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tensor S of a tangent space T x J M. Then the action of the integrand on S is 
obtained as follows: 



N} of the bundle 7r*Tfe«)M. 

2. On this image, the operator A u acts: A u : 7r ] ^ 1 (u) — > ■n^ 1 (u). 

3. The second projection again changes the fiber from 7r ] ^ 1 (n) to the fiber tt' 1 (x). 
However, repeating this procedure for each ti£U, being U an open set. The 
result is not an element of T^f'^M, since there is a dependence on u G T, x . 

4. The integration of this variable provides the desired element, eliminating the 
dependence on u. 

From this short discussion we observe that the geometric interpretation of the aver- 
age operation is quite subtle. We are actually seeking an intrinsic definition, besides 
the general one in [50]. 

One can prove the following fact. The averaging operator acting on a element 
S l ei is the following: 



where [AJ 1 j is the coordinate representation on a given basis of the linear operator 
A at the point u G 7r _1 (x). 

3.2.2 Examples of geometric structures which provide an averaging 
procedure 

1. Lorentzian structures [23]. The geometric data is a Lorentzian metric r] defined 
on M. The disjoint union of the family of sub-manifolds H x C T^M defines 
the fibre bundle tt : S — > M, which we called the unit hyperboloid bundle 



1. First, 7r*5(x) moves S from the fiber tt 1 (x) to the fiber {ir 1 1 (u), u 6 ir 1 (x) C 
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over x £ M, 



£:= [_\{S X CT X M}, H x :={y€T x M\r](y,y) = l} 

The manifold S x is non-compact and oriented. The measure on Sl x is given 
by the following (n — l)-form 

/(x, := f(x, y)^j± dy l A-Ady"" 1 , y° = x\ a;"" 1 , y\ y"" 1 ), 

The function y° defines the parameterized hypersurface S x C T X M, since y° 
can be solved from the condition rj i j^y l y^ = 1. This equation can be expanded 

n— 1 n— 1 

V + 2 ^ V + ( E WV) = L 

a=l a, 6=1 

There are two type of solutions for y°: 
(a) If 7]oq / 0, one obtains a two-fold hyperboloid 



Too v ^ 



a=l 



n— 1 n— 1 



a=l a, 6=1 



(b) If 7700 = 0, the solution for y° is: 

, _i-(EI,hwV) 



2 Ea=l V0aV a 



Note that the Lorentzian metric r\ does not determine the manifolds {~E Xl x £ 
M}. For instance, one can consider £ to be the collection of null cones over 
M: 

£ := [J {NC X C T-jM}, NC, := {y € T,M \ {0} | V (y, y) = 0}. 
xeM 

7r : NC — > M is the null cone bundle over M and NC X is the null cone over 
x; on the other hand, e : NC <^ TM is a sub-bundle of TM — > M. 
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2. Finsler structures [24, 45]. In this case, the Finsler function F(x, y) defines the 
fundamental tensor gij(x,y) = \ which is positive definite, homoge- 
neous of degree zero on y, smooth and lives on the sub-bundle N := TM\{0}. 
The bundle 5] is defined as the disjoint union, 

£:= U^CT.M}, l x :={yeT x M\F(x,y) = l}. 

rrSM 

£ is the indicatrix bundle over M. The manifold I x is compact and strictly 
convex for each x € M and is the indicatrix at x. The volume form is 

u x (u) = dvol(x,y) := ^g^dy 1 A • • • A dy n -\ y° = y°(y\ ...,y n ~ 1 ), 

where the function y° is a solution of the implicit equation F(x, y) = 1. We can 
see that locally this equation has a solution using the implicit function theorem 
and the homogeneous properties of the function F (in particular, using Euler's 
theorem of homogeneous functions). If we take the derivative respect to y° of 
the function (p(x, y) = F(x, y) — 1 and we put it equal to zero, we then get the 
condition of vanishing jacobian: 

£s(F(x,v) - 1) = ^o(fly(x,y)yV - l) = 0. 
Using Euler's theorem one obtains: 

= Qjjo(F(z,v) ~ 1) = (^9o j (x,y)y j +2^^(x,i/))i/V) = 

/ i \ i d .d 2 F 2 (x,y) s , k s 
= (2^(x, y V+^(^-i^) y V). 

Commuting the derivatives and considering Euler's theorem for homogenous 
functions, since F is homogeneous on y, the above expression is 

= (2 90j (x,y)y> + y>" A ( *£<™) )/) = 2 g 0j (x, y)y* = 0. 

Since the metric g is positive definite, the only solution is y = 0, which is 
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outside the indicatrix I x . Therefore, we can apply the hypothesis of the implicit 
function theorem and the equation <p(x,y) = can be solved for y°. 

3. Symplectic structures [25]. In this case, there is defined on T*M a non- 
degenerate, closed 2-form oj. Due to Darboux's theorem [25, pg 246], there is 
a canonical local coordinate system of T*M such that the symplectic form oj 
can be written as 

n-l 

ijj = dpi A dq l . 
i=0 

Associated with oj there is defined on the dual tangent bundle T*M a volume 
2n-form 

S = oj A • • • A oj. 

Using canonical coordinates (q,p), the 2n-differential form can be written as: 
S(p, q) = dp A • • -dp n ~ l A dq° A • • -dq 71 ' 1 . 

Let us assume the existence of a nowhere zero vector field V on TM (therefore 
the Euler characteristic of TM must be different from zero). Then we can 
construct the (2n — l)-form 

Uq (p) = l v S, v e rT(T*M) 

iyS is a non-degenerate (2n — 1) differential form whose value on V is zero, 
since Lz^zS = for any vector Z. Let us chose a distribution of commuting 
vector fields, {X{, [Xi,Xj] = 0, i,j = 1, ...,2n — 1} locally supplementary 
to V such that {V,Xi,...,X2 n -i} is a local frame of TM. The distribution 
{Xi, X2n-i} is integrable and tyS is a volume form on the integral manifold 
S^. On the other hand is a fibered manifold: 

tt:S ± := [J S\ — > M. 
xeM 

Therefore we can define an averaging operation. 
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An interesting thing about this example is that we can only construct local 
averaging procedures. The overlapping of open sets where the averaging oper- 
ation is applied non-trivial and in general one needs more structures to define 
consistently the averaging procedure globally. 

4. Hermitian Vector Bundles [25]. The construction is similar to the one in the 
Finslerian case. The sub-manifolds T, x are defined as: 

:= {yeT x M\H(y,y) = l}, 

where H is the hermitian structure on M. Therefore 

S := □ S x . 

xeM 

Let us assume that the hermitian structure is of the form H = rj + iu, where 
77 is a Riemannian structure and to is a complex structure. To define the 
measure and the volume form we can use either the complex structure uj or 
the Riemannian metric i]. 

3.2.3 Average operator acting on sections 

The averaging operation can be extended to a family of operators acting on sec- 
tions of tensor bundles. This is especially important for the next section. Let 
7r*,7ri,7T2, 7r*T( p ' 9 )M and T^M be as before. Then let us consider the sections 
S e r(T(P><>)M) and tt*S £ r(ir*T( p ' q ^~M.) and the isomorphisms 

ir 2 \ v : r(vr*T( p '^M) — ► r(T^M), S v ^ S z , 

vr* : T{T^M) — »• r(vr*T( p '^M), S z ^ ir* v S z . 
Both isomorphisms are defined pointwise. 
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Definition 3.2.2 Consider the family of fiber preserving endomorphisms 

{A(W) : r(vr*TM) — »• T^*Tm), W G 7r -1 (U), U G M}. 
T/ie averaged operator of this family is the map 

< A >: T(TuM) — > T(TuM) 

suc/i i/iaf af eac/i pomt i S U if is given fry: 

«^>-s)w ; =^(/ N ,^w^*' s )w). 

u G vr" 1 ^), S G TTM, 
where (Ait* ■ S) (u) is the evaluation of the section A(tt* • S) at u. 
A similar definition holds if the operators act on cartesian products of T (7r*T( p ' 9 )M) . 

3.3 Averaged connection of a linear connection on 7r*TM 

We adopt a differential volume form f(x, y) u x (y) such that (dw x (y))\j$ x = 0. There- 
fore we denote w x (y)\isi x = dvol(x,y). Let us assume that a non-linear connection is 
defined on £, with S — > M a vector bundle. 

Definition 3.3.1 Let M be a n-dimensional smooth manifold, n(u) = x and con- 
sider a differentiable real function f G J-~(M). Then ir* f G .F(XJ) is defined by the 
condition 

<f = f{x). (3.3.1) 
The horizontal lift of the tangent vector X l -J^\ x G T X M is 

h : TTM — ► TTN (3.3.2) 
X ^ x ^ Xi ^ "eTr-^x). (3.3.3) 
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Proposition 3.3.2 Let M be a n- dimensional manifold and assume that N is en- 
dowed with a non-linear connection, u G tt^ 1 {x) C N, with x G M. Let us consider 
a linear connection V defined on the vector bundle 7r*TM — > N. Then a linear 
covariant derivative along X < V >x is defined on M, and is determined by the 
following conditions: 

1. MX G T-jM and Y G TTM, the covariant derivative ofY in the direction X, 
is given by the following averaging operations: 

< V >x Y ■=< TT 2 \uV hu(x) 7r* v y >«, G U u , (3.3.4) 

where U u is an open neighborhood of u G Tt~ l {x). 

2. For every smooth function f G JM the covariant derivative is given by the 
following average: 

< V > x f ■=< 7r 2 UV MX) </ > u , Vi> G U„. (3.3.5) 

Proof: it is shown in reference [22, section 4] or in the appendix. □ 

For the physical examples that we are interested, the notion of volume that we 
use is obtained by isometric embedding of the ambient Lorentzian structure n on 
the unit hyperboloid times a positive weight function /. The function f x := f(x, •) 
will be required later to be at least L 1 (S X ) and with compact support S x . This 
implies that the volume function is finite, 

vol(E x ) := / f(x,y)dvol(x,y) < oo. 

The manifold N x is oriented. In particular, the integration is performed in the unit 
tangent hyperboloid, which is 

:= {y G T,M, | V (y,y) = 1, y° > 0}. 
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Note that proposition (3.3.2) also holds in the more general case where the function 
/ is not bounded and does not have compact support, but all the relevant integrals 
(in particular the volume function and average of the connection coefficients) arc 
finite. 

Definition 3.3.3 ( Generalized Torsion) Let V be a linear connection on 7r*TM — > 
N, then the generalized torsion tensor acting on the vector fields X, Y € TM is 
defined as 

Tor(V) : rvr*TM x IVTM — ► IVTM 

(7r*X,7r*Y) — > Tor(V)(TT*X,7r*Y) = V h{x) ir*Y - V h(Y) 7r*X - tt*[X,Y]. (3.3.6) 

This tensor is similar to the usual torsion tensor Tor, 

Tor(V) : TTM x TTM — > TTM 
(X,Y) ^Tor(V)(X,Y) = V x Y -V Y X -[X,Y]. (3.3.7) 

Proposition 3.3.4 The averaged connection < V > has a torsion Tor(< V >) 
such that 

Tor(< V >) =< Tor(V) > . (3.3.8) 

Proof: It is shown in reference [22] and in the appendix of this thesis. □ 

Corollary 3.3.5 Let M be an n- dimensional manifold and V a linear connection 
on the bundle vr*TM — ► M with Tor(V) = 0. Then Tor(< V >) = 0. 

Proof: It is directly shown with the proof of proposition (3.3.4). □ 
If Tor((< V >)) = we say that the connection (< V >) is torsion free. 

Corollary 3.3.6 Let M be an n-dimensional manifold. If the connection V on 
7r*TM has the connection coefficients T l jk(x, y), then the averaged connection < 
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V > has the coefficients 

< V jk > (x) = vd ( Nx) T i jk (x,y)dvol(x,y). (3.3.9) 

Proof: Let {e,}, {7r*ej}, {h(ei)} be local frames for the sections of the vector bun- 
dles TM, 7r*TM and the horizontal bundle % respectively, such that the covariant 
derivative is defined through the relations: 

^h( ej )^* e k = T l jk Tr*ei, i,j,k = 0, 1,2, ...,n - 1. 

Then let us take the covariant derivative 

< V > ej e k = — L_( J^ 2 (V i{e ^e k )dvol x (y)) = — J^ 2 V jk ^e t dvol x 

The relation (3.3.8) follows from the definition of the connection coefficients of the 
averaged connection, 

< V > e . e k (x) =< T >* jk {x) ei {x) 

□ 
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Chapter 4 



Comparison of the Lorentz force 
equation and the averaged 
Lorentz force equation 

4.1 Introduction 

We start this chapter considering several aspects of two different, although related 
topics. The first one is the notion of (semi)-Randers space in the category of met- 
rics with indefinite signature. The second deals with a geometric interpretation of 
the Lorentz force equation and the associated averaged Lorentz equation. Both are 
related and the discussion of the first topic helps to understanding the second. 

The discussion of the above leads to a framework where the results on the averaged 
Lorentz connection can be formulated properly. We also introduce a metric structure 
in the space of connections on some pull-back bundles. Using this metric structure, 
it is possible to compare geodesies of the Lorentz connection and averaged Lorentz 
connection, which is the main result of this chapter. 

The notion of Randers space, which introduces a non-reversible space-time struc- 
ture, can be traced back to the original work by G. Randers [26]. The main moti- 
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vation was the investigation of a geometric structure encoding the time asymmetry. 
One of the results of that study was a unifying theory of gravitation and electrody- 
namics of charged point particles, encoding both physical interactions in an unifying 
space-time metric structure. 

However, one can consider that the non-degeneracy of the associated metric ten- 
sor g(x, y) was not discussed in detail, in particular when the issue of the gauge 
invariance associated with the electromagnetic potential is also considered. From a 
physical point of view, gauge invariance is a natural requirement. The combination 
of this requirement with the non-degeneracy criterion is non-trivial. Combining both 
has lead us to define Randers spaces in the context of pre-sheaf theory (this relation 
between Randers spaces and pre-sheave theory is only slightly treated in this thesis, 
since the details are still under construction). We hope that the identification of 
the appropriate formalism could provide a tool to formulate problems on Randers 
spaces in a consistent way. 

There are other difficulties associated with the signature of the tensor g(x,y). 
Indeed, while for positive definite Finsler metrics there is a satisfactory treatment 
(for instance [24, chapter 11}), for indefinite signatures, the theory of Finsler spaces 
is less universally accepted and several proposals are currently being used in the 
literature. 

There are two general formalisms for indefinite Finsler spaces (that we call semi- 
Finsler structures): Asanov's formalism [27] and Beem's formalism [28], [29]. We 
will argue why both treatments and the corresponding physical interpretations are 
unsatisfactory, in particular when we try to apply them to Randers-type spaces. 
We can see the main problem with Asanov's definition when one considers gauge 
invariance issues related to the structure of Randers-type metrics. Also when one 
considers the possibility of light-like geodesies. The major problem with Beem's 
formalism is that there is no natural definition of Randers-type metric in that for- 
malism. This is because Beem's formalism is based on homogeneous functions of 
degree two in the velocity variables y, while Randers-type functions are by definition 
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homogeneous of degree one in y. 

In section 4-3 we will provide a definition of semi-Randers space which is gauge 
invariant [14]. It has the advantage that all the notions involved are obtained di- 
rectly from the Lorentz force equation and that it is a gauge invariant definition. 
However, this definition does not correspond to a Finsler or Lagrange structure. 
Indeed, we show that even being possible, there are severe practical difficulties to 
find a Lagrangian definition of semi-Randers space which is at the same time gauge 
invariant and globally defined in the tangent space TM. This happens even in the 
absence of topological obstructions like the existence of monopoles for the 1-form 
A. Due to these difficulties we adopted a non-Lagrangian point of view in defining 
semi-Randers spaces. 

In section 4-4 we wm propose a geometric description of the dynamics of one 
charged point particle interacting with an external electromagnetic field. This in- 
terpretation is natural in the framework for Randers-type space discussed in section 
4-3. All relevant geometric data is extracted from the semi-Riemannian metric rj 
and from the Lorentz force equation, which in an arbitrary local coordinate system 
reads 

d 2 a l _ , doi da k N da k I .da da. , . , 

where a : I — > M is a solution curve on M, Tjj, are the coefficients of the 
Levi-Civita connection ''V of r], dA is the exterior derivative of the 1-form A and 
the parameter r is the proper-time of rj along the curve a. Then we interpret the 
equations (4.1.1) as the auto-parallel condition of a linear connection (in a convenient 
bundle). We called it the Lorentz connection. This interpretation does not make 
any additional assumption beyond the information already contained in the system 
(4.1.1) and the space-time metric rj, except for some additional constrains on the 
generalized torsion, necessary to determine the connection coefficients completely. 

In section 4-5 we obtain the averaged Lorentz connection associated with the 
Lorentz connection, applying the averaging method discussed in chapter 3. 
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In section 4-6 we will compare the solutions of the original system (4.1.1) and those 
of the auto-parallel curves of the averaged Lorentz connection. The result is that for 
the same initial conditions, in the ultra-relativistic limit and for narrow one-particle 
probability distribution functions, the solutions of both differential equations remain 
similar, even after a long time evolution, since there is a competition between time 
evolution and other factors. Therefore the original Lorentz force equation can be 
approximated by the averaged Lorentz force equation. 

Remark. The natural object extracted from equation (4.1.1) is what we call 
almost connection (see the appendix for a formalization of the notion). In spite of 
this subtlety, we use through almost all the thesis the name connection (strictly 
speaking projective connection), since most of the calculations that we perform are 
also suitable for almost-connections (or projective almo st- connections) . 

4.1.1 On the physical interpretation of the formalism 

In the next sections we present a formal theory. However, the way the results 
are constructed depends on the physical problems which motivated them. The 
main problem was to model the dynamics of a bunch of particles in an accelerator 
machine, under the action of an external an electromagnetic field. In this chapter 
we will consider the point particle dynamics point of view, which is related with the 
system of differential equations (4.1.1). 

There are some hypotheses in the results that we consider which are related with 
the main problem, although it is not explicitly mentioned: 

1. Bound conditions. For instance, all the parameter and variables that will 
appear in our results are assumed bounded in compact domains K of the 
space-time M. The reason for this assumption is that we are trying to model 
systems like a bunch of particles in accelerator machines. The evolution of 
a bunch starts at a given instant such that t = with the injection and 
separation process of different bunches and the final time t = T, where the 
bunch reaches the target. Every coordinate and parameter of the model is 
bounded in K C M. The external electromagnetic fields are also bounded. 
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2. Ultra-relativistic regime. This is because the type of systems that we are 
describing are ultra-relativistic. We will define an energy function E(x) which 
resembles the energy function used in accelerator physics. The ultra-relativistic 
regime happens when E(x(t)) >> 1, where the mass of the specie of particle 
composing the bunch is set equal to 1 and c = 1. 

3. Narrow distributions. This is one of the characteristics of the bunches in an 
accelerator machine. The narrowness of the distribution function is defined 
through the diameter a of the distribution function in the velocity space. The 
narrowness condition means that this diameter is small compared to the rest 
mass of the charged particles, a << 1. 

4. Adiabatic evolution. It is true that the change in energy is very slow compared 
to the energy itself, once the ultra-relativistic regimen has been reached. This 
is expressed by the condition dl °jj. E << 1. 

The following reasons show why we have adopted the system of differential equa- 
tions (4.1.1) as starting point for our geometrization of the electrodynamics of point 
particles are: 

1. It seems that there is not a satisfactory and simple geometrization metric 
formalism for the interaction of a charged particle with an interacting external 
electromagnetic field (this is the main conclusion of sections 4.1-4.3). 

2. The system of differential equations (4.1.1) is simple, contains all the sym- 
metries that we are interested in and describes all the phenomenology of the 
dynamics of the charged point classical particles. 

3. There exists an standard theory of geometric differential equations and its 
associated non-linear connections. 

4. This geometric theory of differential equations provides the framework to apply 
the geometric averaging procedure described in chapter 3. 
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4.2 Criticism of the notion of semi-Randers space as 
space-time structure 



4.2.1 Randers spaces as space-time structures 

Before moving to the more specific problem of denning semi-Randers spaces, let us 
discuss the notion of semi-Finsler structure. Let M be a C°° n-dimensional manifold, 
TM its tangent bundle manifold with TM D N and with projection ir : N — > M, 
the restriction to N of the canonical projection n : TM — > M. Therefore, (N, tt) 
is a sub-bundle of TM. 

Let us consider the following two standard definitions of semi-Finsler structures 
currently being used in the literature: 

1. Asanov's definition [27], 

Definition 4.2.1 A semi-Finsler structure F defined on the n-dimensional 
manifold M is a positive, real function F : N — >]0,oo[ such that: 

(a) It is smooth in N ; 

(b) It is positive homogeneous of degree 1 in y, F(x, Xy) = XF(x, y), VA > 0, 

(c) The vertical Hessian matrix 

, , ld 2 F 2 (x,y) , t s 

is non- degenerate on N. 

gij(x, y) is the matrix of the fundamental tensor. The set N x is the admissible 
set of tangent vectors at x; the disjoint union N = N x is the admissible 

set of vectors over M. 

In the particular case when the manifold is 4-dimensional and (gij) has signa- 
ture (+, — , — , — ), the pair (M, F) is a Finslerian space-time. 

2. Beem's definition [28], [29] 
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Definition 4.2.2 A semi-Finsler structure defined on the n- dimensional man- 
ifold M is a real function L : TM — > R such that 

(a) It is smooth in the slit tangent bundle N := TM \ {0} 

(b) It is positive homogeneous of degree 2 in y, L(x, Ay) = A 2 L(x, y), VA > 0, 

(c) The Hessian matrix 

a (x v) - l d2L ^ (4 2 2) 

is non- degenerate on N. 

In the particular case when the manifold is 4-dimensional and gij has signature 
(+, — , — , — ), the pair (M, L) is a Finslerian space-time. 

4.2.2 Comparison of Asanov's and Beem's formalism 

Some differences between the above definitions are highlighted below: 

1. In Beem's framework there is a geometric definition of light-like vectors and it 
is possible to derive Finslerian geodesies, including light-like geodesies, from 
a variational principle [30]. By construction, in Asanov's formalism it is not 
possible to do that in an invariant way, because light-like vectors are excluded 
in the formalism from the beginning, since nothing is said about how to extend 
the function F 2 from N to TM. 

2. Let 6(M) be the set of all piecewise smooth paths a : I — > M. In Asanov's 
framework, given a parameterized path a : I — > M, I C R on the semi-Finsler 
manifold (M, F) such that & E N for all t El, the length functional acting on 
a is given by the following expression: 

£ A : 9(M) — ► R 

a(t)^S A (a):= F(a(t),a(t))dt, I = [t mm ,t max ]. (4.2.3) 

Due to the homogeneous condition of the Finsler function F, £a is a re- 
parametrization invariant functional. On the other hand, if we consider Beem's 
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definition, the energy functional is given by the following expression [30]: 



£ B : 0(M) — »• R 

f'tjnax 

a(t) — > f B (<7) := / L(a(t),a(t))dt, I=[UU]' ( 4 - 2 - 4 ) 

Formulated in this way, Beem's energy functional is not re-parametrization 
invariant, because the fundamental function L is homogeneous of degree two 
in y. 

3. A third difference emerges when we consider the category of Randers-type 
spaces: 

Definition 4.2.3 (semi-Randers Space as semi- Finsler Space) 

In Asanov's framework, a semi-Randers space is characterized by a semi- 
Finsler function of the form: 

F(x, y) = ^(x)yV' + A(x)y\ (4.2.5) 

where rjij(x)dx l <g>dx J is a semi-Riemannian metric defined on M and A(x, y) := 
Ai(x)y l is the result of the action of the 1-form A(x) = Ai(x)dx l on y G TM. 

In the positive definite case and when rj is a Riemannian metric, the require- 
ment that gij is non-degenerate implies that the 1-form (A\, A n ) is bounded 
by rj: 

AiAjrfi < 1, rf k r} kj = &j. 

The indefinite case is quite different, since there is not a natural Riemannian 
metric that induces a norm in the space of homomorphisms. Therefore, the 
criterion for non-degeneracy is not trivial and further structure is required. 

Secondly, for both the positive definite and indefinite metric, only the variation 
of the length functional (4.2.3) (for fixed initial and final point variations) 
is invariant under the gauge transformation A i— > A + d\ (not directly the 
integrand itself); the Finsler function (4.2.5) is not gauge invariant as well. 
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Even in the case that we could define a metric and norm, transforming the 
1-form A by a gauge transformation can change the norm and therefore the 
hessian (gij) can become degenerate. 

On the other hand, the notion of Randers space in Beem's formalism is even 
more problematic. In this case there is not a formulation of semi-Randers 
spaces (because eq. (4.2.5) is positive homogeneous of degree one in y). This 
suggests that a proper formulation of the notion of semi-Randers space requires 
going beyond metric structures. 

4. It is interesting to have a definition of semi-Randers structure capable of tak- 
ing light-like trajectories for charged particles into account. As we have seen, 
Asanov's treatment is not able to consider light-like vectors. On the other 
hand, Beem's formalism is not capable of considering this problem, since there 
is not a known Randers-type structure in Beem's formalism. However, the 
asymptotic expansion of the ultra-relativistic charged cold fluid model pre- 
sented in [6] is an example where those light-like trajectories appear naturally. 
In that model, the leading order contribution to the mean velocity field of the 
charged cold fluid is a light-like velocity vector field, interacting with the exter- 
nal electromagnetic field; perturbative corrections change the velocity vector 
field to a time-like vector field. 

The above observations make it reasonable to introduce a non-metric interpre- 
tation for semi-Randers spaces. The option that we have adopted has been to 
formalize a geometric structure from the geometric and physical data that we have: 
the Lorentz force equation and the Lorentzian metric 77. This will lead us to solve 
some of the problems mentioned before. It provides a rigorous framework to discuss 
further developments. 
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4.3 Non-Lagrangian notion of semi-Randers space 

4.3.1 Non-Lagrangian notion of semi-Randers space 

Let us assume the existence of a smooth semi-Ricmannian structure r] on the man- 
ifold M. This implies that the function 

7] : TM x TM — > R 

(X, Y) ' V r ]ij (x)X i Y 3 , X,Y £ T X M 

is smooth in the variables x, X 1 , Y 3 . Since we will use the square root \/t]ij( x ) X l Y 3 , 
we also require that \Jrnj(x) X l X 3 ' is smooth in ''N := UxgTM {X e T X M, %(x) X l X 3 > 
0}. The null-cone is ^NC := UxeM iv e T^M | r](y,y) = 0}. We propose a notion 
of semi-Randers space based on the following 

Definition 4.3.1 A semi-Randers space consists of a triplet (M,r/, F), where M 
is a space-time manifold, rj is a semi-Riemannian metric continuous on TM and 
smooth on TM \ and a 2-form F G /\ 2 M such that dF = 0. 

F is in the second de Rham cohomology group H 2 (M). Due to Poincare's lemma, 
there is a locally smooth 1-form A such that dA = F. Any pair of locally smooth 
1-forms A and A such that A = A + dX, with A a locally smooth real function 
defined on the given open neighborhood, are equivalent and produce under exterior 
derivative the same cohomology class [F] 6 H 2 (M) that contains the element F: 
d(d\ + A) = dA = F. Note that we are speaking of locally smooth 1-forms A and 
of globally smooth 2-forms F. Therefore, instead of giving F, one can consider the 
equivalence class of 1-forms A, 

[A] := {A = A+dX, dA = F in the intersection of the opens sets where A and A are defined}, 

with A a locally smooth 1-form defined on the open set U C M, A a locally smooth 
1-form defined on the open set U C M and A a locally smooth function defined on 
UflU. Then if two locally smooth forms 1-forms ^A and U A, representatives of [A], 
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are defined on and U XJ respectively, one has that that d( ^ A — v A) = 0. For each 
point of the open neighborhood ^XJ = M Un "U, there is a locally smooth function 
defined in an open neighborhood of x XJ(x) C ^ U U such that (^A — u A) = d( tMU X) 
(a consequence of the Poincare lemma). 

Based on these arguments, we give an alternative definition of semi-Randers space: 

Definition 4.3.2 A semi-Randers space consists of a triplet (M, rj, [A]), where M 
is a space-time manifold, n is a semi-Riemannian metric continuous on TM and 
smooth on TM\ and the class of locally smooth 1- forms A is defined such that 

dA = F for any Ae [A]. 

Proposition 4.3.3 These definitions of semi-Randers space are equivalent. 

Proof. We proved already one of the directions of the equivalence. To show the 
other direction, one needs to construct locally 1-forms which produce the required 
2-form F under exterior differentiation. This is achieve by the Poincare lemma in a 
star-shaped domain [32, pg 155-156]. The formula for the 1-form A is 



We will adopt definition 4-3.2, since it has the advantage that it allows a discussion 
of some local issues related with the inverse variational problem of the Lorentz force 
equation. Essentially, this is the reason that even if the topology of M is trivial, the 
1-forms A are only locally smooth. 

As we have learned from the discussion above, a proper treatment of semi-Randers 
spaces combined with gauge invariance requires consideration locally smooth poten- 
tials. There are also locally smooth functions and localcompatibility conditions. 
These kind of structures are formalized by the notion of pre-sheaf structure (and 
the related notion of sheaf structure) [32, 33]. 




I n—1 



□ 
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Let us denote the set of locally smooth functions over M by /\f loc M. This is a 
pre-sheaf structure. The pre-sheave of locally smooth functions on open sets of M is 
denoted by Ti oc {M). Given a Lorentz semi-Randers structure (M,r/, [A]), for each 
of the representatives A G [A] G /\j oc M, there is on M a function Fa defined by the 
following expression: 



The following properties follow easily from definition (4.3.1) and from the definition 



Proposition 4.3.4 Let (M, rj, [A]) be a semi-Randers space with rj a semi-Riemannian 
metric, A G [A] and Fa given by equation (4.3.1). Then 

1. On the null cone NC^ := {y G T X M. \ r]ij(x) y l yi = 0} Fa is of class C° , for 
r] G C° andAe AL M - 

2. The subset where rj x (y,y) / is an open subset o/T x M and Fa is smooth on 
T X M \ NCr, for 7] smooth and A G Azoe M - 

3. The function Fa is positive homogeneous of degree 1 in y. 

Remark. There is no constraint on the non-degeneracy of the fundamental tensor 
gij. Therefore, it is not required that any representative A of [A] be bounded by 1, 
as of is the case for positive definite Randers spaces [24, chapter 11]. 

4.3.2 Variational principle on semi-Randers spaces 

Let *0(M) be the set of piecewise smooth curves on M with time-like tangent vector 
field. The functional acting on a is given by the integral: 



F A (x,y) = < 



y / r]ij(x)y i y j + Ai{x)y l for r] ij (x)y t y^ > 0, 



(4.3.1) 



V~ r lij( x )y i y j + Ai(x)y l for rnj{x)y % y 3 < 0. 



F A - 



£f a : 6(M) 



R 




(4.3.2) 
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with r the proper time associated with n along the curve a. The functional is 
gauge invariant up to a constant: if we choose another representative A = A + dX, 
then £p A {cr) = £f a (c) + constant, the constant coming from the boundary terms 
of the integral. Therefore, the variation of the functional is well defined on a given 
semi- Flanders space (M, n, [A]), for fixed initial and final points variations. 

In order to guarantee the construction of the first variation formula and the exis- 
tence and uniqueness of the corresponding solution it is necessary that the vertical 
Hessian gij be non-degenerate [31]. However, given a representative A G [A], one 
can not guarantee that the Hessian of Fa is non-degenerate. Due to the possibility 
of doing gauge transformations in the representative A(x) \-> A(x) + d\(x) we have 

Proposition 4.3.5 Let (M, n, [A]) be a semi-Randers space. Assume that the im- 
age of the curve a on the manifold M is a compact subset. Then 

1 . There is a representative A G [A] such that the Hessian of the functional Fa 
is non- degenerate. 

2. The functional (4.3.2) is well defined on the Randers space (M, n, [A]), except 
for a constant depending on the representative A G [A]. 

3. If the geodesic curves are parameterized by the proper time associated with the 
Lorentzian metric n, the Euler- Lagrange equation of the functional Fa is the 
Lorentz force equation. 

Proof: There are several steps in the proof: 

1. Using the gauge invariance of £f{o~) up to a constant, we can obtain locally 
an element A G [A] such that \AiAjn l ^\ < 1 in an open neighborhood in the 
following way. Consider that we start with a 1-form A which is not bounded 
by 1. The 1-form A{x) = A{x) + d\(x) is also a representative of [A]. The 
requirement that the hessian of A is non-degenerate is, using a generalization 
of the condition [24, pg 289] 

< 2 Mtfy* + V 7 Hj( x )y i y j (v ij Ai(x)Aj(x)) 
V 7 lij(x)y i y j + A i (x)y i 
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This condition is obtained in [24] relating the determinant of the metric r)ij 
and the determinant of gij, the fundamental tensor of a Randers-type metric 
structure. Although the authors are considering positive definite metrics, the 
result is valid for arbitrary signatures of the metric 77^, for y G N x 

On the unit tangent hyperboloid ~E X this condition reads 



< 



2 + 



A i (x)y i + rfi(x)Ai{x)Aj(x) 



1 + A i (x)y i 

Let us assume (if it is negative, the treatment is similar) that 

«. ( ,,, ) . 2 + M±^lW > a 

1 + Ai(x)y % 

To write down this condition, one needs that 1+ A\(x)y l / 0; the region 
where this does not hold is the intersection of the hyperplane 

P. :={ye T-jM I 1+^ = 0} 

with the unit hyperboloid £ x . The intersection is such that P x n S x C {y G 
Sx I F(x,y) = 0}. 

Let us write in detail the above condition of positiveness for the potential 
Ai = Ai(x) +diX(x): 

2 , , „ (Mx) + diX(x))f + rfi (Mx) + d t X(x)) (Aj(x) + djX(x)) 
e(x,y) = 2 + - 



1+ (Ai(x) + d i X(x))y i 



If y G S x , then for (3 7^ 1, (iy is not in S x . The situation is different if 
y G NC X . Then j3y G NC X , even if /3 / 1. Now we make the approximation 
E x — 7- NC X in the asymptotic limit y° — > 00. In this approximation, one 
can perform limits in the expression for e 2 (x, y). In particular, one can consider 
e 2 (x,/3y) for large j3: 

(Mx) + djXjx))^ + rg (Mx) + ftA(s)) (A,(x) + djXjx)) 
hm e (x,py) = 2+ lim 7 — — — — — — - 

/3^oo v ,Hyj /9^oo 1+ (^(x) +d i A(x))/3y t 
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- 2+ lim (M^)+diX(x))W 3 

That this limit is well defined for large enough y° can be seen in the following 
way. Let us assume that 1 + (Ai(x) + diX(x))y % = 0. Then we consider the 
expression 1 + (^4j(x) + 9jA(x))/3y* for (3 » 1. It is impossible that the 
second expression is zero except if j3 = 1 and 1 + (^4j(x) + <9jA(x))/3y* = 0. 

One should prove that e 2 (x,y) is positive for any y G 5V This is achieved 
because e 2 (x,y) is invariant. Therefore, we can change to a coordinate system 
where y° is arbitrary and the value of the bound does not change. 

2. If the curve a is such that its image <r(I) C M is covered by several open sets, 
for instance M U and Fa{c) is evaluated using both representatives. In 
principle, if we consider two representatives ^A and "A, there is a contribu- 
tion coming from boundary terms coming form the evaluation of dX on those 
points in the intersection M U n "U. This contribution do not appear in the 
first variation of the functional. Therefore, the first variation of (4.3.2) exists 
and does not depend on the representative A, for fixed initial and final point 
variations of a. The corresponding extremal curves exist and they are unique. 
They correspond to the solutions of the Lorentz force equation. 

3. If the image <r(I) C M is a compact set, the number of sets that we need 
is finite. Then we obtain a globally defined section of /\] oc M. 

4. Once a locally differentiable 1-form A with the required properties is obtained 
globally over the variation Var(a) of a : I — > M, one can follow the standard 
proof of the deduction of the Lorentz force equation from the variation of a 
functional [26], [2, pg 47-52]. The variation Var{a) must be constructed such 
that the 1-form A is globally defined on Var(<j). □ 

Remarks 

1. It is important to notice that the non-degeneracy of the metric g, L j does not 
necessarily imply that it has the same signature as the semi-Riemannian metric 



61 



Tj. Further investigations are required to determine the criteria for conservation 
of signature. 

2. If the curve a is parameterized with respect to a parameter such that the Fins- 
lerian arc-length F(a, a) is constant along the geodesic, the geodesic equations 
have a complicated form (see for instance [24, pg. 296]) and they are not in- 
variant under arbitrary gauge transformations of the 1-form A — > A+dX, X G 

^oc(M). 

3. As we have mentioned before, the fact that we are speaking of local data natu- 
ral leads us to consider notions from pre-sheaf and sheaf theory [32, chapter 6, 
33, chapter 2} as the basic ingredient in the definition of Randers spaces. Sheaf 
theory is a theory that allows us to treat local objects, for example germs of 
locally smooth functions or locally smooth differential forms [33, chapter 2]. 
Hence, we think this is a natural framework to study the geometry and varia- 
tional properties of semi-Randers spaces. 

With definition (4.3.2) at hand, the problem of how to introduce the gauge sym- 
metry in a Randers geometry is solved. The price to pay is 

1. The class [A] and the Riemannian metric rj are unrelated geometric objects. 
This is in contradiction with the spirit of Randers spaces as a space-time 
asymmetric structure. 

2. One needs to consider locally smooth 1-forms A instead of globally defined 
1-forms. This happens even if the topology of M is trivial or the cohomology 
class of A is trivial. 

We have seen that given a 1-form A G [A] one has to work hard to find another 
representative A £ [A] such that Douglas's theorem [31] holds. Douglas's theorem 
states the condition under which a system of ordinary differential equations can be 
interpreted as the Euler-Lagrange equation coming from a Lagrangian. One of the 
requirements is that the vertical hessian of the Lagrangian must be non-degenerate, 
det(gij) / 0. 
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On the other hand, light-like trajectories cannot be considered in this formalism 
for semi-Randers spaces in a natural way. Beem's formalism allows the treatment of 
light-like geodesies as extremal curves of an energy functional for some examples of 
indefinite Finsler space-times. However, it is not known if a semi-Randers function 
exists in Beem's formalism. 

In conclusion, the advantage of the definition given here over Asanov's definition of 
semi-Randers space is that it is consistent with gauge invariance. Nevertheless, our 
formalism is still not completely satisfactory, since it can not be used for llight-like 
trajectories. 



4.4 Geometric formulation of the Lorentz force equa- 
tion 

4.4.1 The non-linear connection associated with the Lorentz force 
equation 

Let us consider a semi-Randers space (M, 77, [A]) and the bundle N — > M, N := 
UxeM {y G TrM, rj(y, y) > 0} C TM. Then the following diagram 



TM 



N^^M. 



where e is the following natural embedding 



e : N — > TM 



(x,y) ^ (x,y), x G M, y G N x 

We will consider the differential map dir : TN — > TM. Recall that the vertical 
bundle is defined as the kernel V := ker(dir); at each point one has ker(dir\ u ) := V u , 
with u £ N. 
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The system of second order differential equations (4.1.1) determines a special type 
of vector field on N called spray. It is well known that a spray defines an Ehresmann 
connection on TN[34]. Let us denote by ij(Z,Y) := r\ij{x)Z % YK 

Definition 4.4.1 Let (M, rj, [A]) be a semi-Randers space. For each tangent vector 
y €E T X M with rj(y, y) > ; the following functions are well defined, 

L r jk {x,y) = Tjt^) + - i j F^(x)y m Vmk + F* k (x)y m Vmj )+ 

+ F i m (x) — V (r> jk - , l , VisVkiy s y l ), (4.4.1) 

y '2^y~y) yhk v(y,y) h 1 1 y ! 

T T*jfc(x), (i,j,k = 0, 1, 2, ...,n) are the connection coefficients of the Levi-Civita 
connection T 'V in a local frame, Fjj := diAj — djAi and F l j = rf k F k j, for any 
representative A £ [A]. 

The unit hyperboloid sub-bundle S acquires an induced connection, whose connec- 
tion coefficients are 

L r J jfc (x,y)| s = T« jk (x)+-(F* j (x)y m r ?mfc +F l ^y^^+F* m { X )^{n jk -n ]S ri kl y s y l ). 

The structure of these functions is clear: v T l j k (x) are the connection coefficients of 
the Lorentzian metric 77; the other two terms are tensorial. Indeed, one can define 
the following expressions: 

Vj k [x,y) = -= (F i j (x)y m 7 lmk + F t k (x)y m 7 lmj ), 

2\/r){y,y) 

y m /i \ 

T^x,y) = F^(x)^=={ Vjk - —r, jam rt). 

Therefore, 

1 jk — 1 jk + -* jk + L jk- 

With the functions L l - k and Tj fc , one can construct the following maps: 



L u : T U N x T U N — ► T n N 
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(X,Y)^ U jk (x,y)XiY k -?-. 

Recall that the section -£j is the horizontal lift of the section <9j. 
The second operator that we define is 

T u : T M N x T U N — > T U N 

(X,Y)^ T i jk {x,y)XiY h -^ i . 

u = (x, y) and X, Y are arbitrary tangent vectors X, Y € T U N. This can be 
generalized to homomorphisms acting on vector sections. 
Note the following elementary property 

T U (Y, Y) = 0, Vy G N x , y = z/^-, « = (x, y). 

However, T u (-, y) / in general. 

4.4.2 The Koszul connection L V on TN associated with the Lorentz 
force equation 

Let {eo, e n _i} be a local basis for the sections of the frame bundle associated 
with the tangent bundle TTM and let us assume that each ei is a time-like tangent 
vector at the point x G M (therefore, the metric cannot be diagonal in this basis, 
since r] is a semi-Riemannian metric). Then {7r*eo, vr*e n _i} is a local frame for 
the fiber tt^ 1 C 7r*TM, u G N, {/io, /i n -i} is the local frame of the horizontal 
distribution 1i u C T U N obtained by the horizontal lift hi = h(ei) and {vq, ...,v n -i} 
is a local frame for the vertical distribution V u C T U N. 

Given the set of functions { L T l j k , i,j,k = 0, ...,n — 1} and the non-linear con- 
nection associated with the system of differential equations, there is an associated 
Koszul connection on TN [35] (a Koszul connection is a linear connection defined 
through the corresponding covariant derivative). 

Proposition 4.4.2 Let (M, 77, [A]) be a semi-Randers space and M a n- dimensional 
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manifold. There is defined a covariant derivative L D on TN determined by the 
following conditions: 



1. For each X 6 T~L U and Z G TV. 



L 



D X Z = X kL T i jk {x,y)Zi h t , X = X 1 h t 



Z 



Z l h 



with {hi} a local frame for the horizontal distribution, u = (x,y) and v an 
arbitrary point of an open set OcN containing u. 

2. The covariant derivative of arbitrary sections Z G TTN along vertical direction 
is zero: 



3. The covariant derivative D is symmetric, i.e, has zero horizontal torsion: 

l DjjV — L DyU — [U, V] = 0, VU,V eH. 

4. For all X eV. u and Z £ TV, L D X Z = 0. 

Proof: Let us consider the Finsler geodesic equation associated with the semi- 
Randers space Fa = \fr}ij{x)y l y^ + Ai(x)y : ' but parameterized using the arc-length 
of the Lorentzian metric n: 



where F = dA is the exterior differential of the 1-form A and w(X, Z) = r]ij(x)X l Z 3 . 
From these equations, we can read the value of the semi-spray coefficients, which 
are 



L D v Z = 0, VVGV, ZGTTN. 




dT 2 



dx j dx k 



dr dr 



ff(x, v ) = "F jt (x)yif + r, , i(x)(dA) jk (x)^/ri)h¥)y t . 



Taking the first and second derivatives with respect to y, we obtain 



-^-G\x,y) 
2dv 3 



V^y'+rf 1 ^) (dA) lm (x) 



l 



Vjs (x)y s y m +rf l (x) (dA) y (x) y/r](y, y) . 
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l^W Gl{x ' y) = vr kj{x) - rf l {x){dA)im{ x ) 2Hy ] y)f ,2 u.W v s v m ^) y p + 

+r] il (x)(dA) lm (x) -= =- Vjk(x)y m + rf\x){dA) lk {x) - -= m s y s + 

2(V?7(y,y)) 2( v / r ? (y,y)) 

+rf l (x) {dA)ij (x) Vks (x)y s ■ 

2Wv(y,y)) 

One can check that 3^^G i (x,y) = L r jk (x,y). 

From the structure of these coefficients one can check (following a standard pro- 
cedure, for instance in [34]) that there is a splitting of the tangent vector spaces 
T U N for each u G N. In addition, one can check that the connection coefficients are 
symmetric, L Y % j k = L T l k j. The fact that the covariant derivative along the vertical 
directions is zero is an additional hypothesis used to make the covariant derivative 
unique. □ 

Remark 1. In this interpretation of the Lorentz force equation, the role of the 
function Fa is not fundamental, since all the information is obtained directly from 
the differential equations (4.1.1). 

Remark 2. We have imposed the restriction that the covariant derivatives along 
vertical directions are zero. However, there can be covariant derivatives which are 
compatible with the equations (4.1.1) but which have non-zero vertical covariant 
derivatives. 

Proposition 4.4.3 The following properties hold: 

1. The Lorentz connection L D is invariant under gauge transformations A — > 
A + dX of the locally smooth 1-form A(x) = Ai(x)dx l . 

2. Given a point x G M, L D admits a normal coordinate system centered at x 
which coincides with the normal coordinate system associated with W centered 
at x iffF(x) = 0. 

Proof: 



1. The first property is a consequence of the fact that all the geometric objects 
appearing in the connection coefficients L T l j k are gauge invariant. 
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2. If L D is affine, for any point x G M there is a coordinate system where the 
connection coefficients are zero L Y % jk{x, y) = 0. This implies that L T % j k y^y k = 
at the point x G M. Due to the decomposition L T l j k = T l j k + T % j k + L l j k 
this is equivalent to 

= (T* jfc + TV + U jk )yiy k , G S x . 

Since the transversality condition holds {T l j k yiy k = 0), one obtains 

0=(T < jfc + L i jk )y^y\ \/y G S x . 

Assume that there is a normal coordinate system centered at x for L V and that 
this coordinate system coincides with the normal coordinate system associated 
with W. Then at the point x, one has the relation 

i L T l jk y 3 y k ) = F*,y = 0, Vy G T,M. 
This last condition is strong enough to imply F = 0. 

□ 

Remark. In the positive definite case, it is well known that the requirement that 
the Chern connection of a Randers space lives on the manifold M is that the 1-form 
A must be parallel in the sense that V VA = 0, [24, chapter 11]. This is a stronger 
condition than F = dA = 0. The parallel condition indicates that the structure 
(M, F) is a generalization of the Berwald structure in Finsler geometry: a Berwald 
space is a Finsler space where the connection coefficients live on M; the closeness 
condition indicates that the space is Douglas [24, pg 304]; Douglas's spaces are such 
that they have the same geodesies as the underlying Riemannian metric rj. One 
can easily move the proofs to the Lorentzian and indefinite category, if one adopt 
Asanov's framework. 

Corollary 4.4.4 Let (M, rj, [A]) be a semi- Randers space. Then the Lorentz force 
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equation can be written as 

L D-x = 0, 

where x : I — > M is a time-like curve parameterized with respect to the proper time 
of the Lorentzian metric i], x is the horizontal lift on N and L V is the non-linear 
connection determined by the system of differential equations (4.1.1). 

Proof: A solution of the auto-parallel condition of the Lorentz connection defines 
a curve on N given by (x,y)(r) = (x(t),x(t)). Projecting this curve into M by it, 
one obtains a curve x(t) which is a solution of the Lorentz force equation. □ 

4.4.3 The Lorentz connection on the pull-back bundle 7r*TM 

We introduce the third framework, which will be directly used later to define the 
averaged connection. Given the non-linear connection L D on TN — > N, there is a 
natural linear connection on the pull-back bundle 7r*TM — > N that we denote by 
L V characterized by the following: 

Proposition 4.4.5 The linear connection L V on the pull-back bundle 7r*TM — > N 
is determined by the following structure equations, 

1. L V on 7r*TM — > N is a symmetric connection, 

L V^vr*y - L Vyir*X - tt*[X, Y] = 0, (4.4.2) 

where X, Y G TTM, X, Y G TTN are horizontal lifts of X,Y G TTM to 
TTN, with r](X,X) > and r){Y,Y) > 0. 

2. The covariant derivative along vertical directions of sections o/7r*TM are zero, 

L V V] ir*e k = 0, (j,k = 0,...,n-l). (4.4.3) 

3. The covariant derivative along horizontal directions is given by the formula 

L V hj 7r*e k = L T i jk (x,y)7r*e t , (i, j, k = 0, n - 1). (4.4.4) 



69 



4- By definition the covariant derivative of a function f € .F(N) is given by 

L V^f :=X(f), VleT M N. (4.4.5) 

Proof: One can check by direct computation that the above relations define a 
covariant derivative on 7r*TM and that they are self-consistent. A general covariant 
derivative can be expressed in terms of the connection 1-forms as 

u l j := L V jk dx k + L T jk 5y k , 

where we have used a local frame of 1-forms {dx°, dx n ~ x , Sy°, Sy 11 ^ 1 }. Since 
the covariant derivative of sections on 7r*TM along vertical directions is zero, one 
obtains 

L T jk 5y k = Q L T jk = Q 
at each point (x, y) € N. Since the torsion tensor is zero, one has 

L-pz L-pi 

1 jk kj- 

Therefore, we have to provide the rule for deriving sections along horizontal direc- 
tions. Since the coefficients given by formula (4.4.1) are symmetric, this rule is 
consistent with the torsion-free condition. Finally, we require that the covariant 
derivative to be a local operator. This is satisfied by (4.4.5), which guaranties that 
L V satisfies the Leibnitz rule. □ 

Corollary 4.4.6 Let M, N, 7r*TM and L V be as before. Then the auto-parallel 
curves of the linear Lorentz connection L V are in one to one correspondence with 
the solutions of the Lorentz force equation, 

L V.7r*± = L D ± x = 0, x = %. 
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Proof: If in some coordinate system the Lorentz connection L V has the connection 
coefficients L T l jk, the auto-parallel condition is 

dx(T),y (d 2 x % (T) L , dx(r) dx 1 (r) dx k : (r) 



dr 



0= U^V^TT*-^) = -^F+ Lr ^^) 



dr J V dr 2 dr dr dr 



_/dV(r) / fl 1 ^(r)dx-(r) dx»»(r) 

"I dr 2 + I 1 ^ (rfA,im 2W ^ ) ^l )) 3/2 J?JS dr dr Vkp dr + 

x 1 dx m (T) a,,.. 1 dx s {r) 

+rf\dA) lm -i] jk / ; +rf(tL4)z fc - Vjs —±^ + 

I n ii (dA)l 1 ^(r)x ^(r)dx fc (r) x 

Using rj( dx ^ , ) = 1, the above expression simplifies to 

0=[n( )j = 

/d 2 x l (T) ,7/,^ 1 dx s (r) ,7,,. N 1 dx s (T)\dx j (T)dx k (T)\ 

This is the Lorentz force equation (4.1.1). □ 

We need to translate from the non- linear connection L V on TN to the linear 
connection on 7r*TM because we will consider the averaged connection, which was 
defined in [22] and in chapter 3. 



4.5 The averaged Lorentz connection 

Given the linear connection L V on the bundle 7r*TM — > 5] we can obtain an 
associated averaged connection using the theory described in section 3.3. 

Usually the measure is given as f(x,y)dvol(x,y). The function f(x,y) must be 
gauge invariant and such that the low order moments are finite. The volume form 
is 

dvol(x,y) = y/^detrj -i dy 1 A • • -dy™ -1 , y° = y°(x°, x n ~\ y\ y"' 1 ). 
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Although the dimension of the manifold £ is 2n — 1, we will use the extrinsic 
formalism explained in chapter 2. This in particular means that all Latin indices 
run from to n — 1, if nothing else is stated. Then one can prove the following 

Proposition 4.5.1 The averaged connection of the Lorentz connection L V on the 
pull-back bundle 7r*TM — > X is an affine, symmetric connection on M. The 
connection coefficients are given by the formula 

< L V jk >= T jk + (F* j < l -y m > Vmk + F* k < l -y m > rj mj )+ 

+ F i m ±(< y m > Vj k-v js vki< y m v s y l > )• (4.5.1) 
Each of the integrations is equal to the y-integration along the fiber, 

vol{T, x )= j f(x,y)dvol(x,y), < y % >:= — t4tt / V % f(x,y) dvol{x,y), 
JS X VOL(2j x ) Jv x 

< y m y s y l >:= 1 f y m y s y l f(x, y) dvol{x, y). 

Proof: Equation (4.5.1) follows easily from the definition of the averaged connection 
by linearity. We only need to prove that < y % > and the other moments are given 
by the corresponding integrals and the identity operator Id : 7r*TM — > 7r*TM, 

(x,y) i y (x,y): 

< V 1 >= \ / TT 2 Idy t TT*(y t )f(x,y)dvol{x,y) = 1 / y l f(x,y)dvol(x,y) 

V01{2j x ) J^ x V0L{2j x ) J^ x 

and similarly for other moments. Note that since y G H x r)(y,y) = 1, the factors 
\A?(y> v) do not appear in the connection coefficients. □ 
Remarks. 

1. If we consider the bundle 7r*TM — > N, the coefficients of the averaged 
Lorentz connection are 
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< L T >* kj = T* kj (x) -rf l {x){dA) lm {x) < 2My 1 y))3/2 7 l3 s(x)y s y m Vk P (x)y p > + 

+r?\x)(dA) lm {x) < * =r Vjk(x)y m > +rf\x) < (dA) lk {x)—-^==-ri js y s > + 

+rf l {x){dA) lj (x) < -= Vks{x)y s > . 

In this case, the averaged connection of higher moments of the distribution 
function. 

2. Formula (4.5.1) holds in any local natural coordinate system. 

The following proposition enumerates basic properties of the averaged Lorentz 
connection. The proof is straightforward; one needs to check the properties of the 
coefficients given by the formula (4.5.1): 

Proposition 4.5.2 Let (M, rj, [A]) be a semi-Randers space, f : S — > R non- 
negative with compact support in each T, x and < L V > be the averaged Lorentz 
connection. Then 

1. < L V > is an affine, symmetric connection on M. Therefore, for any point 
x 6 M, there is a normal coordinate system such that < L T l j k > (x) = 0. 

2. < L V > is determined by the first, second and third moments of the distribu- 
tion function f(x, y). 

Remark. While the first property is a general property of the averaged connection, 
the second one is a specific property of the averaged Lorentz connection, which 
results because we are considering trajectories whose velocity vectors y in the unit 
hyperboloid S. Also note that < L v > does not preserves the norm r\ij{x)y % yi for 
arbitrary y G T^M: < L v > z {r]ij{x)y i y j ) / for an arbitrary Z £ T X M. 



73 



4.6 Comparison between the geodesies of X V and 

< L V > 

4.6.1 Basic geometry in the space of connections 

Let us consider the Lorentzian manifold (M,r/) with signature (+, — , — ) and M 
n-dimensional. Let us choose on M a time-like vector field U normalized such that 
rj(U,U) = 1. Then one can define the Riemannian metric fju [36] 

fju(X, Y) := - V (X, Y) + 2 V (X, U) V (Y, U). (4.6.1) 

fju determines a Riemannian metric on the vector space T^M that we also denote 
by fju and that in local coordinates can be expressed as fju = Vij(x) dy l dyK Note 
that the metric fju on M is in local coordinates fju = fjij (x) dx % <8> dx^ . We hope that 
the meaning of the symbols of the type fju is clear from the context. 

The pair (T X M, fju) is a Riemannian manifold. The Riemannian metric fju in- 
duces a distance function d^ v on the manifold T X M., 

d nu : T X M x T,M — > R 

(X, Y) ' y inf{ f 1 Jfju(a,&)dT, a : I — > T X M, «r(0) = X, a(l) = Y}. 
Jo 

We will say that the time-like vector field U defines a field of observer [60, pg 45]. 

We assume that f(x,y) has compact support on each unit hyperboloid S x . Then 
the diameter of the distribution f x = f(x, •) : S x — > R; (x, y) h-> f(x, y) is a x := 
sup{dfj(yi,y2) |yi,2/2 G supp(f x )}. We define the parameter a := sup{a x , x G K}, 
with K C M a compact domain of the space-time M. We will restrict always our 
considerations to K. Note that a (and other parameters) can depend on K. 

Let us fix the coordinate system (x, y) on N. Let us denote the space of linear 
connections on TN by Vn- This space is a finite dimensional manifold whose points 
are coordinated by the set of functions {r* jk(x, y), i,j,k = 0, 1, ...,n — 1}. We will 
introduce a distance function on Vn- First, we recall the definition of the norm of 
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an operator. 



Definition 4.6.1 Given a linear operator A x : T X M. — > T X M in the finite dimen- 
sional normed linear space (T X M, || • \\fj z ), its operator norm is defined by 



llz (x) := sup { llA{ f lfjz (x), y G T,M \ {0} }. 



The norm || • is constructed from the Lorentzian metric rj using a particular local 
time-like vector Z through the definition (4.6.1). We will later specify a vector field 
Z which will be of special interest for our purposes. 

Proposition 4.6.2 On the space Vn, there is a distance function. The distance 
between two points 1 V, 2 V G Vn is given by 

L \fr) Z (X,X) 

X G TTN, X V, 2 V G V N } (4.6.2) 

for a Riemannian metric (4-6-1) constructed using the local time-like vector field Z. 

Proof: The function (4.6.2) is symmetric and non-negative. The distance between 
two arbitrary connections is zero iff 



Vvzi^xX - 2 V X I, l V x X - 2 V X X) = 

for all X G TTS. This happens iff 1 V X X = 2 V X X for any X G TTS. The triangle 
inequality also holds, from the triangle inequality for f). This reads 

, m f ^yzi^xX-^VxX, l v x x- 3V X X) } 

L V i]z{X, X) > 

r \/%( ! Vxi - 2 v x i, x v x x - 2 v x J) 

< sup <^ . ^ + 



y/vz(X,X) 

*v x x, ^ _ 
vV(.v..Y: ' ^ 
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r 2 Vii - 3 v x x, 2 v x i - 3 v x i) -i 

+sup < — = !• < 



^d^Cv, 2 v) + ^ z ( 2 v, 3 v). 

□ 

Recall that for an arbitrary 1-form oj we denote by a;" := r/ _1 (w,-) the vector 
obtained by duality, using the Lorentzian metric 77. Similarly, given a vector field X 
over M, one can define the dual one form X b := n(X, •) ; lx^ is the inner product of 
the vector X with the form u>. Also recall that the difference between two connections 
is a tensor. In this sense, one can consider the difference between connections on 
the pull-back bundle 7r*TM given by the Lorentz connection L X7 and the pull-back 
connection of the averaged connection ir* < L V >. Then it makes sense to take the 
difference L V — tt* < L V >, which is a tensor along ir. We use the hat- notation to 
indicate integrated variables. For instance, < y > means an integration operation, 
where the variable integrated is y. Let us consider a local frame {eo, e n _i} on 
TM. 

Proposition 4.6.3 Let f(x,y) be the one-particle probability distribution function 
such that each function f x has compact and connected support supp(f x ) C S x . Then 

( L V y y-n* < L V > y y)(x) = -(i 5 F)^x)ii y (5))(x,y))+0 2 (5\y))(x,y)+0 3 (5 3 (y))(x,y), 

(4.6.3) 

where 6(x, y) =< y > (x) — y does not depend on the 2-form F. The tensors Oi(x, y) 
are given by the following expressions: 

2 (S 2 (y))(x,y) = ir m ( < y m > {x)5 s {x,y)5 l {x,y)+ < y m >< 5 s (x,y)6 l (x,y) > 

+ 2<y l >(x)< 5 m (x,y)5 s (x,y) > ) %j n lk y i y k vr*^, (4.6.4) 
3 (5 3 (y))(x,y) = ± F l m < 5 m (x,y)5 s (x,y)5 l (x,y) > n^tf y k ix* e t . (4.6.5) 
Proof: From the expressions for the connection coefficients, 

n* < L V > y y- L V y y = \(F^(x)(< y m > {x)-y m )r hnk +F i k (x)(< y m > ( x )-y m ) 11mj ) + 
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,(*)((< y m > (x) - y m ) Vjk - r, jsVkl (< y m y s y l > (x) - y m y s y l )))y j y k ir*e t , 



since r]ijy t y 1 = 1. Here {e^, i = 0, ...,n — 1} is an arbitrary frame unless otherwise 
is indicated. The difference between the two connections can be expressed in terms 
of the following tensors: 

5 rn {x, y) =< y m > (x) - y rn , S msl (x, y)y s y t =< y m y s y l > {x)r] sjmk yiy k - y m 

and is given by the following expression: 

(V< L V > y y- L V y y)(x) = ±(F* J (x)(5 m (x,y))i lmk + F\(x)(5 m (x,y))i lmi + 

+F* m (x)(5 m (x,y) Vjk - r ljs r lkl d ms \x,y)))u j y k 7r*e i = 
= (F j y^6 k (x,y)y k + F i m (x)(6 m (x,y) - %m <5^'(x, y)))^*^, 
where yj = r/ jk y k . 

In the above subtractions the second contribution is of the same order in 5{y) as 
the first one. To show this, recall from the definitions 

S mal (x,y)y a y l =< y m \x)y s (x)y l (x) > W kV ] V k ~ ■ 

Then we can use the following relations: 

f =< y s > (x) - 5 s {x,y). 

One substitutes this relation and taking into account that < S k (x, y) >= 0, one gets 

5 ms \x,y)y syi = ( K <y m >{x)<y s > (x) < y l > (x) + < y m > (x) < 5 s {x,y)5 l {x,y) > + 

+ <y l >{x) < 5 s (x,y)5 m {x,y) > + < f > (x) < 5 m (x, y)5 l (x, y) > - 
- < 5 m (x,y)5 s (x,y)5\x,y) > )r, sj r, lk y j y k -y m . 
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Now we use a similar relation to go further in the calculation. Let us write 

y s =< y s > (x) - 5 s {x,y). 
We introduce these expressions in the calculation of S msl y s yi: 

5 msl (x,y)y syi =( K <y m >{x) (y l + 5 l (x,y))(y s +5 s (x,y))+ < y m > (x) < 5 s (x,y)5 l (x,y) 

+ (y l + S l (x,y))<5 s (x,y)5 m (x,y)> +(y s + 5 s (x, y)) < S m (x, y)5 l (x, y) > - 
- < 5 m (x,y)5 s (x,y)5 l (x,y) > )vsjVlkV i y k - y m . 

Using again the fact that y l i]i k y k = 1 we get (again using < y m >= y m + 5 m (x,y) 
to recombine the first and last term), 

5 msl (x,y)y syi = 5 m (x,y) + 2 < y m > (x) 6 s (y) Vsj y^ + 20?(5 2 ) + 20 3 m (,5 3 ), 

the tensors O2 and O3 are given by the formulae 

Ol(5 2 (y))(x,y) = ±F\ l (<y m > (x)5 s (x,y)5 l (x,y)+ < y m >< 5 s (x,y)5 l (x,y) > 

+2<y l > (x) < 5 m (x,y)5 s (x,y) > ^VsjVlky j y k , 
Ol(5 3 (y))(x,y) = \? i m< 5 m (x 1 y)5 s (x,y)5 l (x,y) > Wk y 3 y k . 
The transversal contribution to the difference between the connections is given by: 

If* m (x) (S m (x, y) Vjk - W ki5 msj (x, y) y V) = ±F j m {x) (5 m (x, y) - 5 m (x, y) 

-2<y m > (x)6>(x,y)r, 8j yi) - (0?(5 2 ) - 3 m (5 3 ))F* m = 
= -F* m (x)( <y m > {x)5 s {x,y)n S3 yi - \o?{5 2 ) - \o?{8*j) = 
= -F* m (x)( <y m > (x)6 s (x,y)y s - \o?{5 2 ) - \oT(5 3 )) = 



> + 
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The longitudinal contribution is 

i(F* j(x)(< y m > (x) -y^Vrnk+F' k (x)(< y m > (x) -y^ mj )y^y k = F%-(x) y^5 k (x, y)y k . 

Adding together the longitudinal and transversal contributions and taking into ac- 
count the formula 5 m (x,y) =< y m > (x) — y m we get the following expression: 

(r^)(5 k (x,y)y k ) - (F m (x)( < y m > (x)6 s (x,y)r ]sj y^ + \o 2 {5 2 ) + \o^)) = 
= -{F* m 5 m (x,y)) (5 k (x,y)y k )- {\o?{5 2 ) + \o^{5*)) F* m . 

□ 

Let us consider a frame {e^, i = 0, n — 1} such that fjz is diagonal at the point 
x G M in this frame. After calculating the distance of the connections using the 
formula (4.6.3), the leading term in S is quadratic. Also recall that in section 1.3 we 
have fixed our units of energy and momentum in such a way that they are given by 
dimensionless numbers. 

Remark. Given a norm fjz on T^M, one can define an induced distance on 
7r^ 1 (ii), !i£ 7r~ 1 (x). This norm is just defined as dfj z (£, C) '■= d^ z 7Ti(£)), £, £ £ 

Proposition 4.6.4 Xei (M, rj, [A]) and L V be as before and assume that f x has 
compact and connected support for each fixed x G M and i/iai a := supja^, x G 
M} << 1. Then the following holds, 

dfj z (TT* < L V>, L V)(x) < ||F||^(x)C(x)a 2 + 2C 2 (x)a 2 (l + a) + Cf (x)a 3 (l + a), 

(4.6.6) 

with C(x), C-2{x), Cs(x) being functions depending only on x with value of the order 
of unity. 

Proof: Let {7r*eo, 7r*e n _i} be a local orthonormal frame for the induced fiber 
metric on the pull-back bundle from the Riemannian metric fj. From equation 
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(4.6.3) one obtains 



L V y y-7r* < L V > y y\\ nz = \\F t J (x)P(x,y)5 k (x,y)y k 7T*e t +02(5 2 )(x,y)+0 3 (5 3 )(x,y)\\ 
< \\F i j P(x,y)8 k (x,y)y k e i \\ fjz + \\0 2 (5 2 )(x, y)\\ nz + \\0 3 (S 3 )(x, 



Each of these three terms can be bounded. 

Recall that we are using a local frame such that ||ej||jj z = 1. Then we can bound 
the first term in the following way: 



For a fixed x the support of the distribution function f(x,y) is compact and 
connected, thus one can write the decomposition < y > (x) = e{x) + z{x) with the 
property z{x) G supp(f x ). In the case that r\ is the Minkowski metric one can check 
that ||e(x)||f; z < a(x) < a by geometric inspection. This bound of e(x) follows from 
the shape of the unit hyperboloid and is proved in the following way. First, note 
that the domain "E x := {y G T X M. \ rj(y, y) > 1 y° > } is a convex set with respect 
to fj. Indeed, we note that 9(S X ) = H x and that < y > (x) G ~E X . Secondly, each 
(£ x ,fjz x ) is a Riemannian manifold, with fjz x '■= (vz)ij(x)dy l ® dyK Therefore, we 
can use the standard definition of center of mass [48], in this case with a measure 
given by f(x, y, s) dvol(x, y). The function f(x, y, s) is such that 



s is the parameter of the line connecting y G supp(f x ) with < y >. Let us denote 
by supp(f x ) the convex hull of supp(f x ). By construction < y >G supp(f x ). One 
can check that < y > is the center of mass of the convex set supp(f x ) [48]. 

We have the following bound 



F i jS j (x,y)d k (x,y)y k e i \\ ijz < ||F|| fe • \\S(x,y)\\rj z ■ \5 k (x,y)y k 




\\5(x,y)\\ flz < \\ < y > (x) - y\\rj z < \\e(x) + z(x) - y\\fj z 



< 
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< IkWIhz + 11*0*0 - Vhz ^ a + a = 2a. 
For the third factor, one has the following bound (S(x, y) =< y > —y) 

\5 k (x, y)y k \ = | < y k > (x)y k -l\ = \ < y k > (x) (y k - < y k > (x) + < y k > (s)) -1| < 

< | < y k > {x){y k - <y k > {x))\ + \ < y k > (x) < y k > (x) - 1\. 
Using Cauchy-Schwarz inequality for f)z we obtain 

\S k {x,y)y k \ < \\<y k > {x)\\ nz \\(y k - < y k > {x))\\ nz +\ < y k > (x) < y k > (x)-l\ < 

< || < y k > (x)\\ij z a+\ < f > (x) <y k > (x)-l| < + \\e\\^ z a +{^ I + \\e(x)\\ nz 

< vTTa a + (a/1 + a - 1) < (1 + a)a + (1 + a - 1) = 2a + a 2 . 
Therefore, one obtains 

\\(F l ^(x,y))(5 k (x,y)y k )e l \\^ z < \\F\\ nz (x)C{x)a 2 + 0(a 4 ). 

The function C(x) in equation (4.6.6) is bounded by the constant 4 in the coordinate 
frame determined by the vector field Z. This bound is universal, independent of the 
Lorentzian metric rj, the vector field Z and it has a geometric origin. 

Using homogeneity properties on the variable y, one can see that the following 
relations hold: 

\\0 2 (S 2 )\\r lz < C 2 (x)a 2 (l + B 2 (x) a), (4.6.7) 

and 

||0 3 (5 3 )|| % < Cl(x)a 3 (l + B 3 (x)a). (4.6.8) 

The functions Ci(x) depend on the particular shape of the support of the distribution 
function / and on the curvature of the metric r]. Using geometric arguments (and in 
particular compactness and connectedness of the supp(f (x)))m K, one can bound 
these functions in terms of a in a similar way as we did for C(x). The constants 
are of order 1 because this was the case for C(x) and there are no new divergence 
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factors in the functions Bi(x) and Cj(x). 



□ 



Corollary 4.6.5 Let (M, 7], [A]) be a (semi)-Randers space. Let us consider a com- 
pact domain K C UxeM su PP(fx) compact, with L V and < L V > as before. Then 
there is a global bound: 

dfj z ( L V,ir* < L V >)(x) < C\\F\\ nz a 2 +2C$a 2 (l+B 2 a) +Cla 3 (l+B 3 a), Vx G M. 

where the constants C, C2, C3, B2, -B3 are bounded by a constant of order 1. 

Proof: It follows from proposition (4.6.4) and compactness of the domain K C 
LLeM su PP(fx) that we are considering. □ 

4.6.2 Comparison between the geodesies of L V and n* < L V > 

There are several ways of defining the energy of a bunch of particles. We have chosen 
one which will be useful for our comparison results. We define the energy function 
E of a distribution / to be the real function 

E : M — > R 

x i-> E(x) := inf{y°, y G supp(f x )}, (4.6.9) 

where y° is the 0-component of a tangent vector of a possible trajectory of a charged 
point particle, measured in the laboratory coordinate frame. The name energy for 
this function is deserved because of the choice of the units that we have adopted, 
even if energy is a function on the co-tangent bundle T*M instead of the bundle 
TM, where the velocities are defined. 

Let us restrict our attention to the case that the Lorentzian metric is the Minkowski 
metric in dimension n. We can define 9 2 (t) = y 2 (t)— < y > 2 (t) and 6 2 {t) =< y > 2 
(t) — y 2 ^). Here y{t) is the spatial component of the velocity tangent vector field 
along a solution of the Lorentz force equation and y(t) is the spatial component 
of the tangent vector field along a solution of the averaged Lorentz force equation, 
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with both solutions having the same initial conditions. The spatial components 
are defined respect the observer Z = J^. We will call this observer the laboratory- 
frame. Since we are fixing this frame, the corresponding Riemannian metric || • ||^ z 
is denoted simply by || • ||^, simplifying the notation. The maximal values of these 
quantities on the compact domain K of the space-time manifold M are denoted by 
9 2 and P. 

Theorem 4.6.6 Let (M, n, [A]) be a semi-Randers space and r\ the Minkowski met- 
ric. Let us assume that 

1. The auto-parallel curves of unit velocity of the connections L V and < L V > 

are defined for time t, the time coordinate measured in the laboratory frame 

7 - 9_ 
z - at- 

2. The ultra-relativistic limit holds: E(x) » 1 for all x € K. 

3. The distribution function is narrow in the sense that a « 1 for all x G K in 
the laboratory frame. 

4- The following inequality holds, 

\e\t) -£ 2 (t)|«i, 

5. The support of the distribution function f is invariant under the flow of the 
Lorentz force equation. 

6. The change in the energy function is adiabatic: ^logE « 1. 

Then for the same arbitrary initial condition (x(0),x(0)), the solutions of the equa- 
tions 

L \7 ± x = 0, < L V >~ x = 

are such that 

\\x(t)-x(t)\\ nz < 2(C(x(t))||F|| f)z (x(t))+C|( a; (t))(l+i? 2 (x(t))a))a 2 ^- 2 (x)t 2 +0(a 4 ), 

(4.6.10) 
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where the functions C{x(t)), Ci(x(t)) and Bi(x(t)) are bounded by constants of order 
1. 

Proof: At the instant t, we calculate the distance measured in the laboratory 
frame between x(t) and x(t), solutions of the geodesic equations of the connections 
L V and < L V > respectively. Both geodesies have the same initial conditions 
(x(0), x(0)). Let us start writing the general expression of the solution of the Lorentz 
force equation for those initial conditions: 

x\t) = x^O) + £ ds(x l (0) + J* dlx l {l)y (4.6.11) 

Since the initial conditions for both geodesies are the same, the equivalent relation 
for the geodesies of the averaged connection is 

x*(t) = ^(0) + J ds(x i (0) + J dlx\l)y (4.6.12) 

We estimate the distance between both solutions at the instant t. Since we know 
the distance between the connections L V and tt* < L V >, it is possible to give a 
natural bound for the distance between the solutions. The main tool that we use is 
the smoothness theorem on the dependence of solutions of differential equations on 
the external parameters (for instance [21, Appendix 1] or [40, chapter 1] or in the 
appendix of this thesis) . 

Let us consider the family of connections depending on the distance between the 
two connections 

Uax = e^( L V,< L V>) 

given by the convex sum: 

f V := -^—(Uax -0 L V + T^—&* < L V >, £ £ [0,^]. (4.6.13) 

Qmax Qmax 

For £ = one has — ^V, while for £ — £max 

one has the averaged connection. 
Using the result of the smoothness of the solutions of the differential equations, 
one can expand the solution ^x 1 of the geodesic equation for the connection with 



84 



parameter £. The second derivative with respect to the coordinate time t reads 

& = °x*+ (<J ; v ; ) ; o-i • Oit 2 l (4.6.14) 

We need to bound the derivative (<9g ^x 1 )\^ = q. The first think is that if x (the tangent 
velocity vector to a Lorentz geodesic) is not on the support of the distribution /, 
this derivative can be done arbitrary large. From the formula (4.6.13) and (4.6.14), 
one obtains: 

(% ^)| 5=0 = 7— — ■ ( L V X ± - vr* < L V >i x)\ 

smax 

where x(t) is the solution of L V±x = with the given initial conditions. This is 
because the support of the distribution function / is invariant under the flow of the 
Lorentz force equation. We are dividing by the distance £ max > thus the derivative is 
such that its norm is bounded by 1, 

ll(^)k=olh =-^—-\\( L V ± x- vr* < L V > x x)^ < 1, 

Kmax 

because of the definition of i ma x and the formula (4.6.14). Note that for writing 
this condition, it is essential that the support of the distribution / must be invariant 
under the flow defined by the Lorentz force equation: if this is not the case, the 
parameter ^ m ax can not be defined and the difference on the covariant derivatives 
cannot be bounded. On the other hand the relations between proper times and 
coordinate time in the laboratory frame are 

dr = 7~ 1 dt, df = ^~ l dt. 

This implies the following relation between derivatives, 

d d d d 

dt ^ dr' dt ^ df 
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Using the hypotheses \9 2 {t) — 9 2 (t)\ « 1 and j^logE « 1, one obtains the 
following relation, 



where by the adiabatic hypothesis, the time derivatives of the energy function have 
been dropped out; the factor 2 comes from the bound of the term which contains 



Corollary 4.6.7 Let (M, fj) be as before and such that there is a global bound for 
||F||^(a;) < ||F||fy < oo, the energy function is bounded from below by a constant E 
and the curves x(t) and x(t) are compact. Then there are some constants C{ such 
that Ci(x) < Ci and the following relation holds 



with F the maximal value of F(x) attached along the compact curves x(t) and x(t). 



1. In the above result the 2- form F is physically interpreted as the Faraday form. 

2. Global bounds occurs in two possible scenarios: 

(a) When manifold M is compact. In this case, one has to consider space- 
like boundaries, since it is well-known that if the space-time is compact 
without boundaries, there exists closed time- like curves ([23, pg 58]) and 
this violates causality. 

(b) The trajectories that we consider are compact in space and bounded in 
time between an initial time t = and a final time t = T. In this case, 
one can define an effective compact space-time manifold with a spurious 
boundary and apply the first case, with the exclusion of closed time-like 
curves, which do not occur in physical situations. 




the derivative of the energy. 



□ 



x{t)~ x(t)\\fj< 2(C\\F\\ fj +Cl(l + a))a 2 E~ 2 t 2 +0(a 4 ), 



(4.6.15) 



Remarks 
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3. We have assumed in our calculations that the external field F does not depend 
on the energy of the beam of particles. However, this is not necessarily the 
case in some situations (as in betatron accelerator machines [9,10]). 

4. There are effects which could reduce the beam size (adiabatic damping and 
Landau Damping [9]). If this happens, there is a strong reduction of the size 
of the dispersion in energy and momenta of the beam. This implies that one 
can describe this as an effective exponent E~ 2+/3 with /3 < 0. Then our results 
are safe under these kind of effects. 

5. That the curves x(t) and x(t) have compact image has a physical interpreta- 
tion: all the trajectories start in a source region and finish in a target region. 

6. Although we have chosen Z = to be the observer corresponding to the 
laboratory frame, the same kind of calculations can be performed in any frame. 

Theorem 4.6.8 Under the same hypothesis as in theorem 4-6-6, the difference be- 
tween the tangent vectors is given by 



i{f)-x{f)\\ f) < (K(x)\\F\\fj(x) +Kl(l + D 2 (x)a))a 2 )E- 1 t + 0(a 4 ). (4.6.16) 



with K(x), K2{x) and D2{x) functions bounded by constants bounded by respective 
constants of order 1 . 

Proof: The proof of this theorem is similar to the proof of Theorem (4.6.6), although 
based on the following formula for the tangent velocity field along a curve: 



Corollary 4.6.9 Under the same hypothesis than in Corollary 4-6.7, there are some 
constants of order 1, K , K2 and D2, such that K(x) < K, K(x)2 < K2, D2(x) < D2 




(4.6.17) 



□ 
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and the following relation holds, 

\\x(t) - x(t)\\rj < {KWFW^x) + K\{\ + D 2 {x)a))a 2 ) E^ 1 1 + C(a 4 ). (4.6.18) 

4.7 Discussion 

4.7.1 Structural stability of the approximation 

L V > 7T* < L V > 

In the proof of proposition (4.6.3) there is a cancelation of the leading orders of 
transversal and longitudinal contributions when we calculate the difference between 
< L V > and L V. The contribution coming from the transversal terms was: 

1^ m (x) (5 m (x, y) Vjk - n js n kl 5 ms i (x, y) y^y k ) = \^ rn {x) (5 m (x, y) - 5 m {x, y) 

-2<y m > (x)6 s (x,y)ri sj yi) - 0?(5 2 ) - 0™(5 3 ) Y i m (x) = 
= -F\ m (x)( <y m > (x)5 s (x,y) Vsj yi) - 0^{5 2 ) - Of{P) F* m (x). 
The contribution coming from the longitudinal terms was 

^(F%(x)(< y m > (x) -y-) r?mjfc+ F i fe (x)(< y m > (x) -y m )v mj ) y J y k = F* j(x) y'5 k (x, y)y k . 

The reason for this cancelation of the first order term in 5 is based on the formal 
structure of the connection L V. This structure has a two-fold origin: 

1. The definition we have adopted for a general non- linear connection and 

2. The structure of the Lorentz force equation. 

The cancelation is independent of the details of the distribution, even if a is not 
small. It is independent of the value of F. 

The cancelation of the linear terms can be written in the following way: 
L V y y -7T* < L V> y y= 0(a 2 (x)). 
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Let us consider the connection 



V = L V - T, T jk = F* m (x) ^== U k - —L-njsriuySy 1 ) . 

The connection V is such that its auto-parallel curves coincide with the solutions 
of the Lorentz force equation. It is also gauge invariant. Therefore, it is also a 
good candidate for a geometrization of the Lorentz force equation. However, if we 
calculate the analogous difference V y y — n* < V > y y we get that in general it is 
linear in 5, which implies 

V y y - it* < V > y y = 0(a(x)). 

Finally, if we consider the covariant derivative X V associated with a spray vector 
field x e TTN, we obtain 

x V y y -7T* < X V >= 0(a(x)). 

We can write this relations as 

lim LV y y -«*< L V>yy = o, (4.7.1) 
a^o a 

lim %V-**<V>vV + o, (4.7.2) 

a-K) a 

Urn W »» -*'<^>,V o. (4.7.3) 
q^o a 

This fact can be stated in the following way. Let us denote a 1-parameter family of 
linear connections by V(q) such that both V(a) and V(a + h) are defined. Then 
one can define the derivative operator 

V(a + h) y y- V(a) y y 
h^O h 
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which is formally the differential (in the way we define, it is the Gateaux differential 
[52, chapter 5]) of the operator 

V(a,y) :I— >T X M 
a ^ L V(a) y y. 

for a fixed y G supp(f x ) and such that a, a + h G I. The notation makes sense such 
that for each y there is a given operator V y (a)y. Then equations (4.7.1), (4.7.2) 
and (4.7.3) can be stated in terms of derivatives. In particular, since V(q = 0, y) = 
L V y y, one has that equation (4.7.1) is equivalent to the statement that L V y y is a 
critical value of V(a,y). Therefore one can write 

2 j2 

k*< L V y y>= 7r*<V(0,y)> +^-_| Q=o V(a,y) +0(a 3 ). 

When the support of the distribution is invariant under the flow of the Lorentz 
equation, for a = one has that L V = vr* < L V >. Therefore, 

2 j2 

< L V y y>= L V y y+^-^\ a=0 V(a,y) + 0(a 3 ). (4.7.4) 

2 ,j2 i 

for any y G S. For arbitrary distribution functions ^- ^| Q , =0 ^( a 5 2/) 7^ 0- 

On the other hand, equation (4.7.2) can be rewritten, following the same steps 

vr* < V y y >= V y y + a -^| a=Q V(a, y) + 0(a 2 ) (4.7.5) 

for any j/eS. For arbitrary distribution functions ^| Q=Q V(a,y) / 0. Also note 
that = L V y y = V y y for any y G ~E, which coincides with the Lorentz force 
equation. 

In the case of a connection X V obtained from an arbitrary spray %, the relation is 
vr* < *V y y>= *V y y+ «^-| Q=0 x V(«,?/) +0(a 2 ). (4.7.6) 
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Apart from the relevance for calculational purposes, it is interesting to know if 
there is some reason why (4.7.4) holds for the Lorentz connection L V obtained from 
the Berwald connection associated to the spray L \- If one changes the way we 
obtain the connection from the spray, or changes the spray (for instance using the 
connection V instead of L V), one obtains conditions of the type (4.7.5). If one 
performs a similar calculation for a general connection, a differential equation like 
(4.7.6) is obtained. This suggests that there are two factors in obtaining the relation 
(4.7.4): 

1. The choice of the non-linear connection as the Berwald- type connection applied 
to the Lorentz spray L \- 

2. The particular structure of the Lorentz equation. This is apparent in the 
calculations in the proof of proposition (4.6.3). 

The conclusion is that the non-linear Berwald connection obtained from the Lorentz 
force equation is structurally stable with respect to the parameter a. As compared 
with other connections, this property happens (maybe only) for the Berwald-type 
connection. This notion of stability is similar to the one presented in [61, chapter 
3}. 

4.7.2 On the hypotheses on which the approximation 

L V — > n* < L V > is based 

We would like to discuss some reflections and interpretations of the hypotheses of 
theorem (4.6.6) and subsequent results. Some of the hypotheses are not essential to 
perform the approximation, but are very useful in the calculations and in writing 
the asymptotic expressions. Therefore, we can differentiate between fundamental 
hypotheses, which are 

1. The auto-parallel curves of unit velocity of the connections L V and < L V > are 
defined for the time t, which is the coordinate time measured in the laboratory 
frame. This is an hypothesis on existence, since nothing can be done if the 
curves are not defined for the parameter that we are speaking about. 
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2. The support of the distribution function / is invariant under the flow of the 
Lorentz force equation. If this hypothesis does not hold, the relation between 
the averaged connection and the Lorentz connection is arbitrary. Therefore, 
this is a fundamental hypothesis. Note that one does not need to assume 
that / is a solution of the Vlasov equation in the sense that the condition is 
weaker. This implies that our result holds for alternative kinetic models for 
the distribution function /. 

3. The dynamics occurs in the ultra-relativistic limit, E{x) >> 1 for all x G M. 
This is a hypothesis which in principle is not fundamental for the calculation 
but it is fundamental for the approximation L V — > tt* < L V > to be good. 
If this hypothesis does not hold, the difference between the solutions of the 
two differential equations will not be as small as in the ultra-relativistic limit 
and the expressions will be more involved. 

4. The distribution function is narrow, a << 1 for all x £ M. This is useful 
to interpret the formulas as asymptotic series in a. Note that for current 
accelerators, this condition holds. 

Apart from discussed above, there are other hypotheses which are not fundamental 
to the results, although they are helpful in the calculations 

1. The following inequality holds 

\e 2 -e 2 \ < i. 

This hypothesis is only used in the proof to simplify some expressions. There- 
fore, it is not fundamental. The approximation L V — > tt* < L V > can be 
good even if \9 2 — 9 2 \ <C 1 is non true, but the asymptotic expressions will be 
more involved. 

2. The adiabatic hypothesis dlo ^ t E << 1 is not necessary for the approximation 
L V — > tt* < L V > to be good. However, this hypothesis simplifies the 
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calculations and the final expression. Note that this is a condition which is 
satisfied in actual accelerator machines. 

There are limitations on the validity of the results that we have obtained in this 
chapter. 

1. The 2-form F is interpreted physically as the Faraday form of an external elec- 
tromagnetic field. We have not commented on the dependence on the strength 
of the electromagnetic field. It is clear that for any finite a, if the electromag- 
netic field is too strong, the approximation will not be good. However, for 
delta distribution functions with invariant support by the flow of the Lorentz 
vector field L x, the approximation is always valid. 

2. The main results of this chapter (theorems (4.6.6) and (4.6.8) and corollaries 
(4.6.7) and (4.6.9)) are not Lorentz covariant. However, as we have said before, 
given an arbitrary observer defined by a time-like vector field Z, it is possible 
to obtain similar results. 

4.7.3 Range of applicability of the approximation 

Let us discuss the limits of applicability of the averaged dynamics and in particular of 
the formula (4.6.10). We will make the assumption that during the time of evolution 
the 7 factors are increasing. Therefore, equation (4.6.10) takes the form 

\\x(t) - x{t)\\f, < Ca 2 ^ 2 (t )||F||^ 2 . 

Note that we are assuming E(t) > E(to). 

Let us assume a natural maximal spatial distance L max between points on supp(f). 
For instance, in an accelerator machine, L max can be related with the diameter of 
the pipe: L max must be smaller than it. Assume that the averaged model is a good 
approximation if the distance ||ic l (i) — x l (i)||^ is less than L max . This can start to 
happen after a time evolution such that the difference \\x l ( 1 t max ) — x l ( 1 t max )\\ f j = 
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Lmax- The characteristic time where the averaged model loses validity is 

There is a second constraint coming from the spread in velocities. Let us consider 
the relations (4.6.16) and (4.6.17). For any distribution function, one can expect 
that the approximation L X7 — >< L V > to be valid until \\x(t) — x(t)\\jj is of order a. 
This is because the unit hyperboloid is strictly convex subset of the tangent space. 
One writes the condition 

a=(K\\F\\ nz {x) +Kl{l + D 2 a))a 2 )E'\t )h max = K\\F\\ nz a 2 E~ l (t a ) H max , 

The maximal time calculated in this way is 

2 t^=K.?M. 1 ± r . (4.7.8) 

The asymptotic behavior of the time where the approximation is broken is universal. 
We have to point out that the approximation can be broken before t max . However, 

Proposition 4.7.1 The following consequences are true 

1. lim£_ ! . 00 t max = oo, if all the other parameters are finite, 

2. lim^^o t m ax = oo, if all the other parameters are finite, 

3. lim|| F ||__>o t max = oo, if all the other parameters are finite. 

The existence of two maximal times up to where the approximation is valid provides 
a definition of maximal length L max : it is the length for which 



) — ■teJ = K - 



E(t ) 



All the quantities appearing are known, therefore one has 



L ma x = A - (4.7.9) 



|F||, 
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with A a constant which does not depend on x, F, a or E. 

If we compare this definition with L max := diam(ir(supp(f))), we have a criterion 
for a definition of weak 2-form: 

Definition 4.7.2 The 2-form F is said to be weak iff L max < L max . 

Proposition 4.7.3 If the 2-form F is not weak, then the approximation L V — >< 
L V > is not valid. 

Proof: If Lmax > L max the system cannot reach the boundary d(ir(supp(f)) before 
the approximation breaks down because the 2-form F produces undesirable effects 
on the dynamics of the system. □ 

Finally, we have the following corollary: 

Corollary 4.7.4 Under the same hypothesis of theorem 4-6-6, the following hypoth- 
esis hold: 

1 . In the limit a — > the Lorentz force equation and the averaged Lorentz force 
equation coincide. 

2. In the limit E — > oo the Lorentz force equation and the averaged Lorentz force 
equation coincide. 
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Chapter 5 



Charged cold fluid model from 
the Vlasov model 

5.1 Introduction 

Despite limitations concerning the mathematical description of the discrete nature of 
the particles comprising a plasma, modeling the dynamics of relativistic non-neutral 
plasmas and charged particle beams by fluid models is common place. The relative 
simplicity of these models, compared with the corresponding kinetic models makes 
them appealing. 

We propose in this chapter another justification for the use of fluid models in 
beam dynamics. We will concentrate on the charged cold fluid model. However, 
we should notice that the same philosophy is also applicable to more sophisticated 
models. 

In high intensity beam accelerator machines, each bunch of a beam contains a 
large number of identical particles contained in a small phase-space region. In such 
conditions, a number of the order 10 9 — 10 11 charged particles move together under 
the action of both external and internal electromagnetic fields. Often in modern 
applications, such bunches of particles move ultra-relativistically. 

One is interested in modeling these physical systems in such a way that: 
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1. The model for a bunch of particles must be simple, in order to be useful in 
numerical simulations of beam dynamics and for analytical treatment, 

2. It allows for stability analysis and a qualitative understanding of the dynamical 
behavior of the system. Three dimensional numerical simulations can be also 
desirable. 

The standard approach has been to use fluid models as an approximation to 
a kinetic model. These derivations of fluid models from kinetic models can be 
found for instance in [8, 37-39] and references therein. They are based on some 
assumptions, usually in the form of equations of state for fluids or assumptions on 
the higher moments of the distribution function /. These constraints are necessary 
in order to close the hierarchy of moments of the distribution function and to have 
a sufficient number of differential relations to determine the remaining moments. 
This is a general feature of all the derivations of fluid models from kinetic theory: a 
truncation scheme is required for the fluid model to be predictive. 

We present in this chapter a new justification of the charged cold fluid model from 
the framework of kinetic theory. The novelty of the new approach is that it uses 
natural hypotheses suitable for particle accelerator machines and exploits only the 
mathematical structure of the classical electrodynamics of charged point particles 
interacting with external electromagnetic fields. We estimate the covariant derivative 
of the mean velocity calculated with the one-particle distribution function. This is 
given as an asymptotic formula in terms of the time of the evolution, diameter of the 
distribution and the energy of the beam. The charged cold fluid model is described 
by only one dynamical variable, the normalized mean velocity field. The variance 
and the heat flow tensor are not necessarily zero, but are finite and given externally. 
In our treatment both the fluid energy tensor and the flux tensor are assumed to be 
given. Our aim is not to give an equation for the mean velocity field, but to evaluate 
how much certain differential expressions (formally equivalent to the charged cold 
fluid model equations) differ from zero. Then one can stipulate the validity of the 
model from the estimates of the corresponding differential expressions. On the 
other hand, in the models presented for instance in [8, 37-39], the variance and 
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the covariant heat flow are dynamical variables and a system of partial differential 
equations is used to determine the dynamics of these fields. We think the analysis 
performed in this chapter, can be extended to other fluids models in future work. 

The method used in this thesis to obtain these results is the following: 

1. One considers the results from chapter 4, which compares the solutions of the 
Lorentz connection with the solutions of the averaged Lorentz force equation. 
In this sense, we are under the same hypotheses as in theorems 4-6.6 and 4-6.8. 
We will use the bounds of the differences of the corresponding geodesies. 

2. It happens that under the same assumptions as used for the particle dynamics, 
the corresponding solutions of the Vlasov equation / and the averaged Vlasov 
equation / are similar. This result is based on the comparison results of the 
point particle dynamics. 

3. Each of the distributions / and / determine a different mean velocity field. 
One can prove under the same hypotheses that these mean velocity vector 
fields are similar. This means that the difference between them is controlled 
by powers of small parameters. 

4. Finally, we show that the auto-parallel condition of the velocity field of the 
averaged Vlasov equation associated with the averaged dynamics is controlled 
by the diameter of the distribution /. Together with the above point, this 
result allows us to provide estimates for the auto-parallel condition of the 
mean velocity field of the solution of the Vlasov equation. 

Therefore, the methods presented here and the usual derivations of the fluid mod- 
els contained in [8, 37-39] are different. The standard approaches assume an asymp- 
totic expansion of the differential equations for the moments, in terms of a perturba- 
tion parameter which is similar to the diameter a of the distribution function. These 
asymptotic expansions are particularized at low orders as a truncation scheme in the 
hierarchy of moments. Those references discuss systems of partial differential equa- 
tions which are self-contained and consistent with physical constraints and with the 
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asymptotic expansions. On the other hand, our approach is based on the structure 
of the Lorentz force equation of a charged point particle, which lies at the basis 
of the kinetic models. Written in a geometric way, the Lorentz force equation is 
replaced by the averaged Lorentz force equation. The key point is that the averaged 
Lorentz connection admits normal coordinates, which simplifies in a fundamental 
way our calculations. Then under some regularity assumptions on the distribution 
function, we can place bounds on the differential expression of interest. 

We have assumed that the distribution functions are smooth (at least of class C 1 ) 
in the coordinates x % . Although we do not currently have a proof that we can extend 
our results to bigger functional spaces for the distribution functions, since the main 
results are written in terms of Sobolev norms, it is conjectured that they can be 
extended to Sobolev spaces. Indeed our proofs suggest that we require smoothness 
on the x coordinates and the existence of weak derivatives on the y coordinates. 

5.2 Comparison of the solutions of the Vlasov and av- 
eraged Vlasov equations 

In this section we estimate the difference between the solutions of the Liouville 
equations associated with the averaged Lorentz connection and the original Lorentz 
connection. The Liouville equation associated to the Lorentz force equation is called 
Vlasov's equation in kinetic theory. In a similar way, we call to the Liouville equation 
associated with the averaged Lorentz force equation the averaged Vlasov equation. 

5.2.1 Examples of Liouville equations 

Given a non-linear connection characterized by the second order vector field % G 
TTN, the associated Liouville equation is %(/) = 0. In the following two examples 
presented here, which are related with our purposes. 

1. From the coefficients of the Lorentz connection L T l jk(x, y) one can recover the 
spray coefficients L G l (x, y), using the homogeneous properties on y of L G 1 (x, y) 
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and Euler's theorem of positive homogeneous functions. In particular, the 
spray coefficients are 

L (?(x,y)= L r i jk (x,y)yiy k = (V jk + —t==(F i j (x)y™r ]rnk + 

v m 1 \ 

+F\(x)y m Vmj ) +F l m {x)—^==(ri jk - -^—^ Vjs 7 lkl y s y l )jy^y k = 

= hr jk + -j^== (F* j(x)y m Vmk + F* k (x)y m Vmj ))y>y k = 

= ^r, fc yV + AAK^y)F%(xV. 

Then one can define the vector field L x- 

L x{x,y)=y 1 ^-- (^r% fe (x)yV + ^y)F%(xV) A 

2. A similar procedure applies to the averaged Lorentz Vlasov vector field. In 
this case, however, we do not have the simplifications given above. Therefore, 
the spray coefficients are 

< L G i > (x,y) = L r jk (x,y)yJy k = ("I* jk (x)+ < - J- (F l j(x)y m r] mk + 

v ^Vv{y,y) 

+F l k (x)y m rj mj ) +F t m (x) 2 ^l__ = {rj jk - -^j^Y) > )v 3 y h - 
If y G S, semi-spray coefficients can be simplified to 

< L G" > (x,y)| s = (V jfe (x)+ < 1(F* j (x)y m r hnk + F i k (x)y m r hnj )+ 

+F J m (x)— (77j fe - vjsr] k iy s y l ) >)y 3 y k = 

= + ^(F l ,(^) < > + F* k (x) <y m > Vmj )+ 

+F i m {x){i ljk < — > -r] jsm < —y s y l >)J y j y k . 
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The averaged Vlasov vector field can be written in a similar way as before, 

< L X > Is = y'~ ( v T i jk y j y k + \{? l 3 {x) <y m > r^+F* k (x) < y m > r, mj )+ 

+F t m (x)(r ]jk < — > -vjsVki < -yi/V >)) y J y k q-- 

5.2.2 Comparison between the solutions of the Vlasov equation and 
the averaged Vlasov equation 

In the following (M,ry) is Minkowski space, since we will use theorems (4-6.6) and 
theorem (4.6.8). The Riemaniann metric f)z is determined by the vector Z = J^. 
Since we will use the Euclidean metric associated with Z = Jj, we simplify the 
notation and employ f) in place of f)z- Z = ^ will be the observer that we call 
laboratory frame. We will restrict our attention to a compact domain K C M. 

Proposition 5.2.1 Let f and f be solutions of the Vlasov equation L x{f ) = and 
the averaged Vlasov equation < L \ > (/) = 0, where L \ and < L \ > are the spray 
vector fields obtained from the non-linear connections L V and < L V >. Let us 
assume the same hypotheses as those in theorem 4-6.6. Then for the solutions of the 
Vlasov and averaged Vlasov 's equation with the same initial conditions, one has the 
relation 

\f(t, x(t), x(t)) - f(t, x(t),x(t))\ < (C(x)||F||^f (x)(l + B 2 {x)a))a 2 E~ 2 1 2 + 

+ (K(x)\\F\\fi(x) iff (1 + D 2 (x)a))a 2 E~ x t (5.2.1) 

for some functions C(x(t)) K(x(t)) along the geodesic of the Lorentz connection. 

Proof: / and / are solutions of the corresponding Vlasov and averaged Vlasov 
equations respectively. Therefore, / and / are constant along the corresponding 
auto-parallel curves; x(t) and x(t) are the projections on the space-time manifold 
M of the integral curves of the vector fields L \ and < L \ >■ Here t is the time- 
parameter in the laboratory frame determined by the vector field 4. Then the 
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Vlasov and averaged Vlasov equation can be written as 

L xf = j t f(x(t),x(t)) = 0, < L x > f = j t Hm,m) = o- 

For the same initial conditions, the geodesic curves corresponding to the connections 
L V and < L V > are nearby curves at the instant t in the way described by theorem 
4.6.6. 

Let us introduce the family of interpolating connections, 

L V e := (1 - e) L V+ e < L V >, ee[0,l]. 

Each of them has an associated spray vector field L \e- Therefore, let us consider 
f e (x,y) to be the solution of the following Liouville equation L Xefe = for some 
given initial conditions. Since the dependence on (e, x, y) of the vector field L Xe is C 1 , 
the solutions of the Liouville equation are Lipschitz with respect to the parameter 
e. We can see this fact in the following way. The Liouville equation can be written 
as 

L X J e = &j t f(x e (t),y £ (t))=0, 

where (x e (t),y e (t)) is an integral curve of the vector field L \e restricted to the unit 
hyperboloid bundle and such that it is parameterized by the coordinate time t. Then 
one can use standard results from the theory of ordinary differential equations to 
study the smoothness properties of the solutions of the above equation. In particular, 
the connection coefficients for the interpolating connection are, 

( L r e )v = (i-e) L rv + e< L v Jk >. 

From the formula (4.5.1) for the coefficients L T l jk{x,y) one can check that ( L T e ) 1 jk 
are smooth functions in an open set of time-like vectors y and the parameter e. From 
here it follows the Lipschitz condition for f € in e. 

We will give an upper bound for the difference \f(t,x(t),x(t)) — f(t,x(t),x(t))\. 
Note that in this expression the point where both / and / are evaluated are (t, x(t),x(t)). 
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In order to achieve this, standard results on the smoothness of the solution of differ- 
ential equations are used (see chapter 1 of [40]). In particular we use that for each 
(e, x(s),y(s)), there is an open neighborhood U f of [0,1] x supp(f) T - K containing 
(e,x(t),y(t)) such that the solutions of the differential equations are Lipschitz in Uj. 
Therefore, using the Lipschitz condition, one obtains the bound 

\f(t, x e (t),x t (t)) - f(t, xi(t),xi(t))\ < c i(e, x(t),x(t))S((e, x(t), £(*))) + 

+c 2 {e,x(t),5(t)) \\x e (t) - Xi(t)\\fj + c 3 (e,x(t),x(t)) \\x € (t) - xz(t)\\fj. 

Cj(e, x(t),x(t)) are constants which depend on the open neighborhood XJf, <5((e, x(t), x(t))) 
is the diameter on the e component where we are applying the Lipschitz condition. 

One can always choose a refinement of an open cover of [0, 1] x s«pp(/)|tk such 
that both the Lipschitz condition, theorem (4.6.6) and theorem (4.6.8) can be applied 
simultaneously. Since [0, 1] is compact, we can consider a finite open covering of [0, 1] 
for each instant t. Then using the above local bound in each of the open sets U e , 
one obtains the global bound 

\f(t,x(t),x(t)) - f(t,x(t),x(t))\ < Cl +c 2 \\x(t) - x{t)\\ n + c z \\x{t) - Z{t)\\ n . 

The constants q are finite (by definition of Liptschitz and by compactness of the 
interval [0,1]). The functions / and / are constant along the respective geodesies. 
Therefore, 

\f(t,x(t),x(t)) - f(t,x(t),i(t))\ = |/(0,x(0),x(0)) - /(0,£(0),£(0))|. 

Let us assume the same initial conditions x(0) = x(0) and x(0) = x(0) for the 
geodesies of the Lorentz connection. Since the difference \ f(t, x(t),x(t))—f(t, x(t),x(t))\ 
is a smooth function of \\x(t) — x(t)\\fj and \\x(t) — x(t)\\fj, one obtains 

< ci < ifi||z(t) " X{t)\\fj + Ki||x(t) - X(t)\\rj 
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for some constants K^. Then we have, 

\f(t,x(t),x(t))- f(t,x(t),x(t))\ < 

< \f(t, x(t),x(t)) - f(t, x(t), x(t))\ + \ f(t, x(t), x(t)) - f(t, x(t),x(t))l 

The first term is bounded by c±, which is bounded by Ki\\x(t) — x(t)\\fj + Ki\\x(t) — 
x(t)\\fj. The second term can be developed in Taylor series in the differences \\x(t) — 
x(t)\\fj and \\x(t) — x(t)\\fj, since / is smooth. Therefore, 

\f(t, x(t),x(t))-f(t, x(t),x(t))\ < (C(x(t))\\F\\ n (x(t)) Ci(l+B 2 (x(t))a))a 2 E' 2 t 2 + 
+ (K{x(t))\\F\\ fj (x{t)) K 2 (l + D 2 (x{t))a))a 2 E' 1 1. 

□ 

5.3 The charged cold fluid model from the averaged 
Vlasov model 

In the following results (M, rj) is Minkowski space, since we will use theorem (4-6.6) 
and theorem (4-6.8). The Riemannian metric fjz is determined by the vector Z = ^. 

There is another local observer related with the vector field < y > . Since the norm 
is not continuous on M, one needs a local smoothing procedure. Given a subset of 
the paracompact manifold, S x , we can take the induced bump function from the 
bump functions defined on S x . Using these bump functions, we can smooth vector 
fields [51, pg 25]. 
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5.3.1 Comparison of the Vlasov model with the averaged Vlasov 
model 

Definition 5.3.1 Given a semi-Randers space (M,r/, [A]), the averaged Vlasov model 
is defined by the dynamical variables f determined by 

< L X>f = 0, (5.3.1) 

where < L x > is the Liouville vector field of the averaged Lorentz dynamics as- 
sociated with the external electromagnetic field F. The dynamical variable f(x,y) 
defines the following 

V ■= yf(x,y)dvol(x,y), vol(S x ) := / dvol(x,y)f(x,y). (5.3.2) 

Since we will use the results of chapter 4, the distribution function / is at least of 
type C 1 in the x-coordinates and Lipschitz on the y-coordinates. Since the support 
/ G fx is compact, several Sobolev norms are defined [41, chapter 3\. We will write 
our results in terms of those norms. 

Proposition 5.3.2 Let < L x > f{ x iV) = an d L xf( x iV) = be such that the 
domain of definition of the vector field < L x> is an open sub-manifold o/S. Then 
one can reduce supp(f x ) — > supp(f x ) for all x 6 M. 

Proof: Let us consider the product of the functions f(x,y)g(x,y), where < L x > 
f = 0, the function g(x, y) is a bump function adapted to the support in K of the 
vector field < L x >■ Since both the support of < L x > and supp(f x ) are sub-sets 
of the paracompact manifold X, this function exists [32]. Therefore, we select the 
function g(x, y) such that 

9x(y) = 0, (x, y) e supp{f x ) \ U (f x ), g x {y) = 0, (x, y) G dsupp{f x ) 

where U{f x ) D supp(f x ) and all the derivatives are zero on dsupp(f x ). Then one 
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can perform the following calculation: 



< L X> (/<?) = 9 < L X > f + / < L V > g = 0. 

We can always restrict the solutions of < L \ > f = in such a way that formally 
supp(f x ) = supp(f x ) and points 2 and 3 are proved. □ 

Using equation (5.2.1), it follows that the error induced by the substitution / — > 
f is of order a 2 . Hence, in the following calculations, when it is useful, we can use 
dvol(x,y) as a measure and substitute supp(f x ) by supp(f x ) and / by /. 

Let F be a closed differential 2-form defining the Liouville vector field L \- Let us 
consider the Sobolev spaces (W 1 - 1 ^), || • ||i,i) and (W 0,2 ^), || ■ ||o,2> [41]. Recall 
that the space of smooth functions is denoted by F{T> X ) (an introduction to the 
notions of Sobolev spaces can be found in appendix 5 or [40, chapter 3]). 

Recall that we have denoted 6(x,y) =< y > (x) — y. For the next results we 
will restrict to the Minkowski space-time (M,r/). Let us denote T(supp(f x )) the 
characteristic function of supp(f x ). 

Theorem 5.3.3 Let M be an n- dimensional space-time manifold, K C M a com- 
pact domain and < L \ > the vector field associated with the averaged Lorentz force 
equation. Assume that: 

1. The distribution function is such that f x , djf(x,-) G T(E X ) C W 1 ' (S x ), 

2. The function 5(x, •), (djS)(x, •) G W ' 2 ^). 

Then 

i 

\\ < L V > v n*)k < "f^f^ (E \\9olo 9 (S k x )\\o,2) ■\\f x \\ 1 , 1 -a*+0(a% 

vol(supp(f x )) 

(5.3.3) 

where S x (-) := S(x, •) and 

V\x) =< f > f (x) := - r — 1 —— r f dvol(x,y)f(x,y)y i . 

1 h x f(x,y)dvol(x,y) 
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The volumes are 



vol(supp(f x )) := vol(E x ); vol E {^ x ) ■= / T(supp(f x )) ■ dvol(x,y); 
the derivative in equation (5.3.3) refers to the local frame such that the vector U := 



<y> 



(U o ,0). 



y/v(<y>,<y>) 

Proof: Because the averaged Lorentz connection is an affine connection on M, given 
a point x 6 M, there is a coordinate system where the connection coefficients are 
zero at that point, < L T > l jk{x) = 0. Therefore, for any given point x € M 
one can choose a normal coordinate system such that the averaged Vlasov condition 
holds, 

y j d j f(x,y)\ x = 0. (5.3.4) 

Using this normal coordinate system, one can get a simplified expression for the 
covariant derivative of V along the integral curve of V: 



using a coordinate frame {-^k, k = 0, ...,n — 1}. Note that this expression is not a 
partial differential equation because it only holds at the point x. 

From the relation (5.3.5) we obtain that 

< L V >y V(x) = V%V k = —3— f dvol(x,y)yif(x,yy 

,/V n / dvol(x,y)y k f(x,y)). 
It is the right hand of this equation that we shall estimate, 

,:» (/ dvol(x,y)y J f(x,y)- 
vol{supp(f x )) Ws, 

l^ dvol ^y)y k ~^y))) = 
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a u ^ iif \ I k x dvol(x,y)f(x,y)y k 
dvol(x,y)y J J{x,y) { 



vol(E x )\Jv x ' ' vol 2 (£ x ) 

■djij dvol(x,y)f(x,y)))^+———(f dvol{x,y) y 3 f{x,y) ■ 

■dj{ f dvol{x,y)y k f{x,y))) = 

= ~ ^2( S ) ( dvol(x,y)y :} f(x,y)d j (J^ dvol(x,y) f(x,y) ) < y k > ) + 

■dj([ dvol{x,y)y k f\x,y))). 

Shifting the variable of integration —y k + < y k >= —5 k (x,y), one obtains the fol- 
lowing for the above expression 

< L V > v V(x) = - —pzj^ ^ { j dvol(x,y)y j f(x,y)- 



■dj( f dvolyf(x,y)) <y k >) 



+ 



■dj ( I dvol(x,y) (< y k > +5 k (x,y))f(x,y)))d k = 

= vol 2(p ) ( j?, dvol ( x ^v)y 3 h x ^v)dj{j^ dvol(x,y)S k (x,y) f(x,y))^jd k . 

Since y l dif(x,y) = at x and since f x is a smooth function of y, we can Taylor 
expand the integrand, obtaining: 



■dj ^oZ(x,y)<5 fe (x,y) (/>,y) + J| - </) 
There is a coordinate system such that < L r > l j k {x) = at the point x. This is 
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reflected in the averaged Vlasov equation, which has the form yi djf(x,y) = at 
one given point x 6 M. Then we get the following for the above expression 



( < L D >y V(X] 



1 



vol 2 C£ 



dvol(x,y)y 3 f(x,y) ■ 



■a,-(jf dvol(x,y)6 k (x,y)^(y l -y 1 ))). 

{y l — y l ) an d S k (x,y) are bounded by the diameter a(x) (remember that in taking 
the moments we can substitute the pair (f x , supp(f x )) by (f x , supp(f x )) if we desire, 
since by proposition (5.2.1) the difference between the two distributions functions 
is small and because by proposition (5.3.2) we can replace the supports as well). 
Therefore, 



< L V >yV 



dvol{x,y)y 3 f(x,y) 



v vol 2 (£ x ) 

■dj ( / dvol(x, v) 5 k (x, y) % (y l - y l )d k ) 
J-E x oy 

< ^ dvol(x,y)y 3 f(x,y)d j ^ i 

dvol(x,y)5 k (x,y) (y l - y l )d k ^j 



< 



vol 2 CE a 



+ 



voP(Z x ) 



df\ 

dvol{x,y)y 3 f(x,y) — x j 



< 



dvol{x,y) dj5 k (x,y) (y l - y l )d k ^j 



< 



vol 2 (X x ) 



df \ 

dvol(x,y)y : > f(x,y)dj ) 



dy l 



vol 2 C£ x ) 



dvol(x,y)\\5 k (x,y) (y l - y l )d k ) + 

r) / 

df \ 

dvol(x,y)y J f(x,y) — % J 



dvol(x,y)d j \\6 k {x,y) (y l - y l )d k \\ _ ) . 
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One can find a bound for each of these integrals. For instance, using the Hoelder 
inequality for integrals in an arbitrary space X [13] 



X(j) dfi 



< 



(| x iA^r) 1/p (| x i0^i 



, M V9 1 1 

, - + - = 1, 1 < p,q < oo. 

p q 



We will use this inequality several times for the case p = q = 2, obtaining 



dvol(x,y) (y l -y l )5 k (x,y)d k ) 



< 



< 



< 



( / dvol(x,y)\(y l -y l )\\\ 5 k {x,y)d k \\.) 

< n dvol(x,y)\(y l -y l )\ 2 f ■ ( f dvol(x, y) \\ S k (x, y)d k f. ) \ 

Note that the index I is contracted with a factor J^j-. Therefore y l J-j/j- is Lorentz 
invariant and it can be computed in any inertial system, in particular in the labo- 
ratory frame. If we do this computation on this frame, we can continue with the 
above bound in the following way: 

( f dvol(x,y)\(y l -y l )\ 2 ) 1 • ( f dvol{x, y) \\ 5 k (x, y)d k \\ 2 _ ) 1 . 

< voll(Y, x ) « • ( / dvol(x, y) \\ 5 k (x, y)d k 2 . 

In order to bound the second factor we use the following argument (that we used 
already in section 4-6), 

1 3 

\\5{x,y)\\fj < || < y > (x) -y\\fj < \\e + y-y\\fj < ||e||^ + ||y-y||^ < -a +a = -a. 
y is in the support of the distribution /. Therefore, a bound on the integral is 



( f dvol(x,y)(y l -y l )6 k (x,y)d k )\\_ < | • vol E (Z x ) ■ a 2 . 



no 



Similarly one obtains the following bound: 



< 



dvol(x,y)\(y l -y l )\ 2 y 



< 



dvol(x,y)d J (5 k (x,y))(y l -y l )d k 

■ ( f dvoiix^Wd^ix^mwif 

< voll(V x )a- ( ! dvol(x,y)\\d j (5 k (x,y))d k \\iy. 
Then because of the definition of the corresponding Sobolev norm || • || 0,2 : 



dvol(x,y) WdjVMMtf-MdkW-) < voll( S upp(f x ))a- \\djS%, 2 
voll^a ^ \\6%log(6 k x )\\ ,2 < volj,^) a 2 ■ £ \\d 3 log(5 k x )h t2 . 



Similarly, 



dvol(x,y)y 3 f(x,y)dj —j^j < ^ ( / ^oZfx, 2/) |y 5 /(»> 2/) | 2 ) ' 



dvol ( x >y)\ d j^k 



df(x,y) ,2 \ i 



)' 



The second factor is equal to the Sobolev norm H^/a; ||i,i - The first factor is bounded 
in the following way: 

dvol(x,y) |y J f(x,y)\ 2 ^j 2 <( / dvol(x, y) \y ] | 2 f(x, y) ) 2 = 

= toI(S,) •(< |^| 2 >)i 
Therefore, we get the bound: 

n— 1 

( / d w /(x,y)^/(x,y)a,|4)| <«o/(S :c )(^(<|^| 2 >)5)||^/ :c || 1 , 1 . 
In a local frame where the vector field U = (U°, 0), this contraction can be re-written 
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as 



(< (vT >) 5 • \\doM\i,i = ll(< \vT >) 5 • fhfxhi = ll(< (vT >) 5 ■ 3>/*lki = 

= ||(< (y ) 2 • (d f x ) 2 >)t\\ ltl = ||(< (y* • a,/,) 2 >)*||i,i. 

The last expression is covariant. Using normal coordinates associated with the affine 
connection < L V > we obtain ||(< (y J • djf x ) 2 >)^||i,i = 0. 

Finally, we can bound the following integral 



df 

dvol(x,y)y J f(x,y) — { 



< 



df 

dvol(x,y)\ y J f(x,y) J < 



<£(/ dvol{ X ,y)Wf{x,y))\ 2 ) h -{l dvol{x,y)\^pfi\^) h 

7 '_q -'Sx y 



As in the previous integral, we get 



5/\ 



dvol{x,y)y J f{x,y)^ 1 j\<vol('E x )- (< |^'| 2 >)5 • \\f x \\ hl 



dy 



Using these bounds, we obtain the following relation: 



< L D >y V(x)\\ n < 



( dvol(x, y) y J f(x, v) d j^i) \ vol E {supp{f x )) 



•a 2 + 



n-l 



+ 



^l(£,)-a 2 -(E \\djlog(6 k x )\\o,2) < 
k 

< .}._ . -Vol(Ji x )- 



vol 2 C£ x ) 



n-l 



> 3 



3=0 



ii,i • ^4(^)-« 2 -(E ii^M^)iio, 2 )) = 
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v -§Hf • ( E« w\ 2 »* (E m^im) • n/.iii,! • « 2 - 



j=0 fc 

In a local frame where the vector field £7 =< Y > has components (C/°,0), the 
following relation holds: 



E(< ly'l 2 >)' (E ll^iM^) llo,2 = « IvT >)* (E II* ||o, 2 = 

= Ell < ^°) 2>5 /o ^)ll0,2. 

fc 

Note that the normal coordinate system (that we are using) coincides with the 
adapted coordinate system, associated with the vector field U = (U°,0). In this 
coordinate system, there is a bound < y° >< 1 + a, where a is the diameter 
measured in the co-moving frame. It is of order 1 or smaller than 1. Therefore, 

E(< W\ 2 >)* (E Wog{S*) ||o, 2 < (< |y°| 2 >)* • (E HftM*£) Ita ■ (l + fi). 

j=0 k k 

Then we have the following result: 

1 

LV7 ^ T7/--MI / ^K^V 



< 



□ 



Corollary 5.3.4 For compact domains K C M and under the same hypotheses as 
in theorem 4-3, the following relation holds: 

|| < L V > v V(x)\\a < n ■ C(K) ■ a 2 + 0(a 3 ), 

for some constant C(K). 

Proof: Take the constant C(K) to be 



r vol 2 (Yl 1 1 
C(K) = max x6K { (£ \\d log(5 k x )\\ , 2 ) ■ \\f x \\ ltl }. 
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□ 

These expressions are asymptotic formulas if 1 >> a. 
Remarks 

1. In the preceding results the ultra-relativistic limit (E >> 1) was not essential. 
However, the series in power of the energy has asymptotic meaning if E >> 1. 

2. There are several notions of normal coordinates, since we have several affine 
connections: "?V, W and < L V >. However, we have only used the normal 
coordinates associated with < L V >. 



5.3.2 Bound on the auto-parallel condition of the unitary mean 
vector field of the averaged Vlasov model 

Let us consider the normalized mean velocity vector field: 

„ _ V 
U ~ T](V, V) 1 / 2 ' 

Since < L V > does not preserve the Minkowski metric rf, the covariant derivative 
of u in the direction of u using the Lorentz connection L D is 

< L V> il u=^j- ) < L V>yV + ±(v-(log(r l (V,V))))v. (5.3.6) 

The first term is bounded by theorem 5.3.3, since r](V, V) > 1. The total derivative 
of r)(V, V) along a trajectory of V is 

Cy( V (V,V))=V-( V (V,V)) = 

= 2r ? (< L V >y V,v) + (< L V >yv)(V,V). 

We have proved that the first term is of order a 2 . Using normal coordinates for 
< L V >, one can compute the second term: 

( < L V > v v ) (V, V) = V (V, V)F jm < 5 rn (x, y) 5 s (x, y) S l (x, y) > V^V s V t . 
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We can estimate these contributions 



Proposition 5.3.5 Under the same assumptions as in theorem 5.3.3, the following 
relation holds: 



< L V>uU< V °lf^j (E H<9oM^)llo, 2 ) ■||/«||i,i-« 2 +0{c?). (5.3.7) 



Proof: The first term of the right hand of the equation 5.3.6 is bounded by theorem 
5.3.3. The second term is bounded using Hoelder's inequality for integrals [42, pg 
62] 



/ dvol{ Z ) h{z)- f m {z)\ <II (/ dvol(z)\f k (z)\ Pk ) Pk , £ft = 1, 1 < Pk < 
^ X k=l JX k 

In particular one can apply this inequality to the third order moment 

< 5 m (x,y) 5 s (x,y) 5 l (x,y) >: 



oo. 



< 5 m (x,y)5 s (x,y)5 l (x,y) > 



< 



dvol(x, y) f(x, y) 5 m (x, y) 5 s (x, y) 5 l (x, y) 

S p 3 j(/ a _*»«».»)i/(x,,)ff"(»,,)i') 1 - 

•( I dvol(x,y)\f(x,y)5 l (x,y)\ 3 y. 

The distribution function is positive on supp(f). Also one can choose a distribu- 
tion function such that f x < 1. By proposition (5.3.1), one can substitute in the 
integrations f x — > f x , which implies that 



<6 m (x,y)S s (x,y)6 l (x,y) > =0(a 3 ). 
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Since the norm n(< y >,< y >) > 1, one gets a third degree monomial term in a 
for the covariant derivative < L D >u u. □ 

Corollary 5.3.6 Under the same assumptions as in theorem 5.3.3 in a compact 
domain K C M, one obtains 

< L V >a u < •C(K) • a 2 + 0(a 3 ), (5.3.8) 

for a convenient constant C(K). 



5.3.3 Bound on the auto-parallel condition of the mean velocity 
field of the Vlasov model 

Let us consider a local Lorentz congruence, which is a set of auto-parallel curves of 
the Lorentz connection L V, for a set of initial conditions at each (to,x), x G M to 
where M <0 > M is a 3-dimensional spatial sub-manifold. One can consider in 
a similar way the congruence associated with the averaged Lorentz connection for 
the same initial conditions. Note that, while the Lorentz connection preserves the 
Lorentz norm r)(x, x) of the tangent vectors of the geodesies, this is not the case for 
the averaged Lorentz connection. 

Theorem 5.3.7 Let F be a closed 2-form and L V the associated non-linear Lorentz 
connection. Under the same assumptions as in theorem 5.3.3 the solutions of the 
Lorentz force equation v V±x = (tiF)f can be approximated by the integral curves 
of the normalized mean velocity vector field u = v ( x ) ^ = j ^ e distribution 

y/T)(V(x),V(x)) 

function f(x,y), where f(x,y) is a solution of the associated Vlasov equation L \f = 
0. The difference is controlled by polynomial functions at least of order 2 in a, 

|| L V u n||(x) < a 2 (x)a 2 +0{a 3 ) (5.3.9) 

where the function a 2 {x) is a bounded function of x. 

Proof: We repeat an argument that we have used before. By proposition 5.3.1, 
both distribution functions / and /, solutions of L \f = and < L \ > f, are such 
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that 



\f(t, x(t), x(t)) - fit, x(t),x(t))\ < ((7(x)||F^C 2 2 (x)(l + B 2 {x)a))a 2 E~ 2 t 2 + 

+ (K(x)\\F\\ fj (x) K 2 (l + D 2 (x)a))a 2 E' 1 1. 

Therefore, the corresponding mean velocity fields are nearby as well, because of the 
linearity of the averaging operation and because of the above relation. Then their 
corresponding integral curves and the associated local congruences are also similar. 
By corollary 5.3.6, for narrow distributions, the normalized mean field u associated 
with / is such that 

|| < L V >u u(x)\\fj < a 2 a 2 + 0{a 3 ). 

for some function a 2 (x). Remember that one can interpolate smoothly between the 
connections L X7 and < >. Therefore, locally, one can interpolate smoothly 
between their integral curves. Also, because of the smoothness of the solutions of 
the geodesic equations with respect to the parameter of interpolation, there is a 
function a 2 in a small open neighborhood of M such that 

L X7 u u < a 2 (x)a 2 +0(a 3 ). 

□ 

5.4 Discussion 

Theorem 5. 3. 7 shows when the charged cold fluid model is a good approximation to 
the Vlasov equation in the description of the dynamics of a collection of particles 
interacting with an external electromagnetic field, in the ultra-relativistic regime. It 
is interesting that we have obtained this result without using additional hypotheses 
on the higher moments of the distribution function, except that the distribution is 
narrow and smooth enough for our calculations (we need some smoothness conditions 
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in order to use Taylor expansion for the function f(x, y) in the velocity coordinates. 
Indeed, it seems that one only requires weak differentiability in y). 

One of the important hypothesis on which the calculation is physically relevant 
is the requirement that the diameter of the distribution a must be small in the 
laboratory frame. Also, the ultra-relativistic regime E >> 1 is useful in order to 
have good estimates and holds in current particle accelerators. 

There are some technical issues that we would like to mention briefly: 

1. We have assumed that the distribution functions / and / are at least C 1 in x. 
However, let us consider the Dirac delta distribution with support invariant 
by the flow of the Lorentz force, 

f(x,y) = V(x)6(y-V(x)). (5.4.1) 

Since the width of the distribution is zero, a = 0. One can use this distribution 
as a solution of the Vlasov equation in 2-dimensional space-times, for a proper 
value of the function VP. This example and the fact that the bounds found in 
section 5.3 are formulated using Sobolev norms suggest the possibility of gen- 
eralize the results to bigger function spaces. The results to use in this case are 
Sobolev embedding theorems [41, 57, 58]. However, we are not investigating 
this question in this thesis. 

2. The same method can be applied to other fluid equations. Depending on the 
specific bounds and parameters, one can decide which model is better in each 
particular situation. 

5.4.1 On the validity of the truncation schemes in fluid models 

Given a kinetic model, usually the Maxwell- Vlasov system of differential equa- 
tions, one defines a fluid model in terms of the low moments (typically, first, 
second, third and fourth moments) of the distribution function / [8, 37-39]. 
All the higher moments are set equal to zero. Usually, the typical reasoning 
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is that with the low moments one can write down models that are consistent 
and explain a reasonable number of phenomena in plasma physics. 

We can argue that the reason why these models work is that in some situations 
the underlying Vlasov model can be substituted by the averaged Vlasov model. 
Then the Vlasov model depends only on the first, second and third moments 
of the distribution function /. Therefore, as soon as the the hypothesis of 
a given fluid model are compatible with the hypothesis of the approximation 
Maxwell- Vlasov model — > averaged Maxwell- Vlasov model, fluid models whose 
dynamical fields can be written in terms of the first, second and third moment, 
are equivalent to the underlying averaged Vlasov model. The equivalence must 
be understood in an approximated way, since there is an approximation in this 
argument. 

The variables that one considers in fluid models are the mean velocity field 
(2.1.6), the covariant kinetic energy-momentum tensor (2.1.7) and the covari- 
ant energy-momentum flux tensor (2.1.8). Therefore, one can propose the 
following 

Definition 5.4.1 Two kinetic models are equivalent if their corresponding 
mean velocity field, covariant kinetic energy-momentum tensor and covariant 
energy-momentum flux tensor are the same. 

Definition 5.4.2 Two fluid models are the same if their corresponding mean 
velocity field, covariant kinetic energy-momentum tensor and covariant energy- 
momentum flux tensor are the same. 

We propose the following conjecture in the form of a theorem 

Theorem 5.4.3 If two kinetic models are equivalent, the corresponding fluid 
models are the same. If two fluid models are the same, the underlying kinetic 
models are equivalent up to the order of approximation of the kinetic model by 
the averaged kinetic model. 
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The first implication is trivial. The second implication is true for the Vlasov 
model, as we have proved in this chapter. 

Therefore, when one works with a kinetic model, there is an underlying equiv- 
alence class of fluid models. We can call this an universal class. The elements 
of an universal class are, by construction, fluid models of ultra-relativistic 
narrow distributions. Then it is a useful idea to consider for each class the 
simplest model possible. In practice, the simplest model will dismiss higher 
order moments. 
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Chapter 6 



The Jacobi equation of the 
averaged Lorentz connection 
and applications in beam 
dynamics 

6.1 Introduction 

6.1.1 The Jacobi equation of an affine connection on M 

Let M be an n-dimensional manifold. Given an affine connection V on the tangent 
bundle ir : TM — > M, the auto-parallel curves of V are the solutions c : I — > M 
of the system of differential equations 



V T T = 0, 




where t is an affine parameter of V. The curvature tensor R of the affine connection 
V is the tensor field defined by the expression 

R : TTM x TTM x TTM — > TTM 
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(X, Y, Z) i ^ R(X, Y, Z) = V X V Y Z - V Y V X Z - V [X ,Y]Z, VA, Y, Z G TTM. 

(6.1.1) 

The curvature tensor has an associated family of curvature endomorphisms {R x (X, Y), X, Y G 
rTM, x G M} defined by 

#(A, Y) : rTM — ► rTM 

i? x (X,y)Z = R(X,Y,Z)(x), Z G rTM, G M. 

A vector field J along the parameterized geodesic c : I — > M of an afflne, torsion- 
free connection is a Jacobi field if it satisfies the Jacobi equation 

dc 

V X V X J-R(X,J)X = 0, X = —. (6.1.2) 

This equation can be re- written using a local frame {ej, i = 0, 1,2, ...,n — 1}. It 
corresponds to the system of the second order differential equations 

n2 ji J yj J vm 

^_ R . j , m(c(t))J ^^ = , (6.1.3) 

where is the covariant derivative along the reference geodesic and the curvature 
tensor is given by the expression 

R jkm = d m Y d}fT j m + (r jaT rrn r jmY r fc) . 

For the next fundamental results one can consult [35, section 10.1; 54, § 14]. It 
is well known that for an affine connection there are 2n linear independent Jacobi 
fields along any given central geodesic c. This fact is a consequence of the existence 
and uniqueness of the solutions to second order differential equations and their 
smoothness properties on the initial values. A Jacobi field is completely determined 
by the value of J(0) and ^(0). 
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A smooth map 



C : (—A, A) x I — > M t— >-c s (t) C M 

is a geodesic variation of c(t) if each c s (t) : (—A, A) — > M is an affine parameterized 
geodesic curve for all s 6 I. The variation vector field of a geodesic variation is the 
vector field along the curve c defined by dC (■§;), with dC : T((-A, A) x I) — ► TM 
the differential of the smooth function C. The vector field [dC(J^), ^|r(c(£))] along 
the curve c(t) vanish. The variation vector field acts on an arbitrary smooth function 
as a derivation 

The geometric interpretation of a Jacobi field is obtained through the following [35, 
54]: 

Proposition 6.1.1 Let M be a manifold equipped with an affine torsion- free con- 
nection V. Then each Jacobi field J is a variation vector field of a geodesic variation 
C. Conversely, any variation field of a geodesic variation C defines a Jacobi field. 

6.1.2 Jacobi equation for linear connections defined on the pull- 
back bundle vr*TM — > N 

Let M be a smooth n-dimensional manifold, N TM a sub-bundle of the tangent 
bundle TM and let us consider an Ehresmann connection defined on TN. 

In order to formulate a Jacobi equation for linear connections on 7r*TM — > N, 
we mimic the standard derivation of the Jacobi equation for affine connections [35] . 
The bending term is determined by the /i/i-curvature endomorphisms: 

R{h{X),h{Y)) : 7r*TM — »• vr*TM 

C ^ R(h(X), h(Y))( = {V h{x) V h{Y) - V h{Y) V h{x) - V M[Xl Y]))C, 

VC e TTM, 
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where the horizontal lift h, h(X l di) = X l {x)-^ was introduced in section 3.3. 

Similar to the case of affine connections, one obtains the Jacobi equation for linear 
connections on 7r*TM — > N. A basic fact is that we require that the Jacobi vector 
field commutes with the vector field X(t) along the curve c(t), which means that 

[J,X]\ c{t) =0. 

Then the torsion-free condition along the curve c is 

Vft(j)vr*A" = V A( x)vr* J, 

and the curvature endomorphism along the curve c is such that 

R(h(X), h(J)) C = {V h(x) V h(J) - V h(J) V h(x) ) C, C € IVTM. 

We can compute the second covariant derivatives along the curve c, 

V M x)V M x)7r*J = V h{x) V h{J) n*X = V h{J) V h{x) TT*X - R(h(X), h(J)) ir*X. 

(6.1.4) 

Let us assume that X G TTM is a auto-parallel respect to V, which means V/j(x)vr*X 
0. One obtains the following second order differential equation for J(t) 

V M x)V MX)7 r* J - R(h(X), h(J)) tt*X = 0. (6.1.5) 



Definition 6.1.2 A field J(t) along the curve c : I — > M satisfying equation (6.1.5) 
is a Jacobi field o/V. The corresponding vector field iT2( J (t)) along the curve is the 
associated vector field. 

We can identify the vector field 7T2(J(t)) with J(t). 

Definition 6.1.3 Let V be a linear connection on the bundle 7r*TM. An auto- 
parallel variation of the auto-parallel curve c : I — > M is a map C : (— A, A) x I — > 
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M such that for each value of the parameter s, the curve c s (t) := C(s,t) is an 
auto-parallel curve, V^(c s (t))7r*c s (t) = 0. 

Proposition 6.1.4 Let V a linear connection on 7r*TM. The variation field along 
c(t) of a variation C(t,s) is a Jacobi field ofV. 

Proof: It is clear from the deduction of the Jacobi equation for V. □ 

Remark. For covariant derivatives such that they are zero along the vertical 
direction, the covariant derivatives does not depend on the particular lift of X G 
T^M to T U N. Therefore for those covariant derivatives the expression that one 
obtains in this case for the Jacobi equation is defined as before and is independent 
of the vertical component of the lift that we are using. 

6.1.3 Physical Interpretation of the Jacobi equation 

For an affine connection, the Jacobi field represents the deviation vector of a given 
trajectory from the reference geodesic. However, if the reference trajectory is ob- 
servable, the assumption that the reference trajectory coincides with the central 
geodesic provides physical meaning to the Jacobi field. Let us consider a geodesic 
variation C(s,t). Each of the geodesies c s (t) = C(s,t) corresponds to a possible 
trajectory for a charged point particle. Then the Jacobi field corresponds to the 
deviation variable from a particular trajectory follow by a particle to the reference 
trajectory. 

The above property holds at least for affine connections. Therefore, let us fix a 
central geodesic c(t) and consider the set of all geodesic variations of the central 
geodesic c(t). Hence there is a relation between the set of geodesic variations (which 
is equivalent to the set of Jacobi fields along c(t) by proposition 6.1.1) and the set 
of all the trajectories allowable by the dynamics and the topology of the space-time 
manifold M. 

If the Jacobi vector field is given by J(t), one identifies the components of J(t) 
with the relative coordinates of a given geodesies respect the central affine geodesic 
c(t), J k (t) = u k (t). 
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Assuming this interpretation, the Jacobi Equation (for both affine and non-affine 
connections) is a second order Riccati equation 

-^+R{t)u = 0. 

This type of equation appears when one considers small deviations from a solution of 
another differential equation. One example is Hill's equation in celestial mechanics 
[55]. 

However, in the case where the connection is affine, the form of the Riccati equa- 
tion is the same, but the endomorphism R(t) along the curve c(t) is simpler, since the 
curvature endomorphism depends only on the point c(t) and not on the derivative 
^jp- as in the general case of a non-linear connection. 

6.2 Jacobi equation of the averaged lorentz connection 

In this section M is a n-dimensional manifold. The averaged Lorentz connection 
< L V > is an affine connection on M. 

Using local coordinates, a Jacobi field along the central geodesic can be written 
as J = &(s)dj. The reference trajectory will be X(t), which will be assumed to 
be a geodesic of the averaged Lorentz connection. Then the Jacobi equation for an 
affine connection on the tangent bundle TM — > M can be expressed as 

+ 2T ifc (X( T ))__ + £ dl V jk ( X )—— - o. 

This equation is called the geodesic deviation equation. The central geodesic is 
denoted by X{t) and a neighborhood geodesic is given by x(r) = £(r) +X(t). The 
parameter r is the proper-time of the central geodesic measured with the metric rj. 

The averaged Lorentz connection < L V > is an affine connection on the tangent 
bundle TM. Therefore, we can apply the standard Jacobi equation to the averaged 
Lorentz connection. Given an arbitrary semi-Randers space (M, 77, [A]), in a lo- 
cal natural coordinate system, the averaged Lorentz connection has the connection 



126 



coefficients 

< L F jk >= jk + i(F* <y m > i] mk + ¥ l k < y m > Vmj )+ 

+F i m (< y m > Vjk - VjsVkl < y m y s y l > ), 

with F = dA, with A being a representative of [A], A € [A]. The tangent vector 
y are on the unit hyperboloid S x . 

The Jacobi equation of the averaged Lorentz connection is 

+F* m (< y m >v jk -VjsVki < y m y s y l > ))+2^(^(F* < y m > r^+F* fe < y m > 7? mj )+ 
+F i m (< y m >r/ jfc - ^ S 77 fc , < y m y s y l > )) ^ ^ + 

From the form of the system of differential equations (6.2.1) we conclude that: 

1. There is a term representing the inertial acceleration: 

AM'T V + { Vr^(^ + 2^f). (6.2,, 

The inertial acceleration Ax is universal, in the sense that it is independent of 
the particle mass. 

2. We can not say that Ax is independent of the electromagnetic field, since ^p- 
can depend implicitly on the electromagnetic field when defining the reference 
trajectory. The typical example is the reference orbit of a betatron [9, chapter 
3}. 

These properties of the inertial acceleration Ax will help us to write it without doing 
explicit calculations, using establish formulae for elementary cases. 
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6.2.1 The Jacobi equation of the Lorentz connection versus the 
Jacobi equation of the averaged Lorentz connection 

We have shown that the deviation equation from a given reference trajectory defines 
a Jacobi equation. However the non-linearity of the Lorentz force equation creates 
difficulties in view of the applicability of the above interpretation: 

1. The evaluation of the covariant derivatives with respect to the original Lorentz 
connection requires a reference vector. Due to the dependence of the connec- 
tion coefficients F l jk(x,y) on the direction y, there must be assigned a partic- 
ular point in the tangent space yo G T^M at which the connection coefficients 
should be evaluated. This implies a specification of the direction where the 
connection coefficients are evaluated. Attaching a physical significance to the 
choice of the vector yo implies the selection of a particular model, which re- 
quires additional justification. This difficulty is resolved using the averaged 
Lorentz connection, which is an affine connection and whose connection coef- 
ficients do not depend on the reference vector yo- 

2. It was proved in [22] (although in the category of Finsler spaces and for con- 
nections which covariant derivative vanish along vertical directions) that the 
averaged curvature of the original connection is the curvature of the averaged 
connection. This result can be extended to arbitrary linear connections on 
7r*TM with vanishing covariant derivative in the vertical directions. Hence, 
we can apply this result to the Lorentz connection. In the corresponding av- 
eraged Jacobi equation appears the averaged curvature, which is the same as 
the curvature of the averaged connection. Therefore we can work with the Ja- 
cobi equation for the averaged connection, as an attempt to give an averaged 
description of the dynamics of a beam of particles where the bending term is 
the averaged force. 

Motivated by the above reasons, in this chapter we replace the Jacobi equation of 
the Lorentz connection by the Jacobi equation of the averaged Lorentz connection as 
a description of beam dynamics of bunches of particles in accelerators. We consider 
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systems in the ultra-relativistic regime and we also assume that the distribution 
function / is narrow in velocity space. We show how the Jacobi equation for the 
averaged connection < L V > provides a geometrical formulation of the transversal 
(in the case of dipole and quadrupole fields) and longitudinal beam dynamics (when 
the external fields are linearizable in the relative coordinates). We then provide a 
method to introduce corrections to the averaged Lorentz dynamics caused by the 
composed nature of the bunch of particles. These corrections are expressed in terms 
of known or observable quantities. 

6.3 Transversal beam dynamics from the Jacobi equa- 
tion of the averaged connection 

Let us assume that (M, 77) is Minkowski space with M being 4-dimensional. There is 
a global coordinate system denoted by (r, x 1 , x 2 , x 3 ). r is the proper time considered 
from a given initial point of the reference trajectory, the coordinate x 2 is given by 
the Euclidean length of the path of the reference trajectory measured from the initial 
position in the reference frame defined by the vector field jl, which corresponds to 
the laboratory frame, x 1 is the horizontal coordinate and x 3 the vertical coordinate 
respect to the central geodesic. The longitudinal direction at each instant r is given 
by the vector By definition (x^x 3 ) are the transverse coordinates, while x 2 is 
the longitudinal coordinate. 

6.3.1 Relation between the transverse dynamics and the Jacobi 
equation 

We choose the laboratory reference frame for our calculations. Under the transverse 
dynamics, the difference ^ — will be constant. The external electromagnetic 
fields are static magnetic fields in Minkowski space. In the subsequent calculations 
we will only consider the lower order terms in the degree a + b + c of the monomials 
£ a (^) 6 e c appearing in the expressions, with £ = x{t) — X(t), e =< y > — 

Firstly, we linearize the equations with respect to the degree defined by the vector 
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fields along the central geodesic £ and the powers of the difference e appearing on 
each term Recall that the transverse component, which are the terms in the geodesic 
equation proportional to the tensor T l jk = F* m (< y m >r/jk — rjj s r]ki < y m y s y l > ), 
are neglected systematically because they are of higher order in the degree (a + fe+c) 
than the longitudinal component. The longitudinal component is proportional to 

V jk = -(F* j < y m > V mk + F* k < V m > Vmj) 

After linearizing, the system of differential equations are 

d 2 C d£idX k (l.: _ , ,\ 



+2£'$ (1 (F* < y m > r? mfe + F l fc < y m > Vmj ) 



dXi dX k 
dr dr 



+ (T, fc + ^r'r, fc )(^^ + 2^^)=0, ij,k,m = 0,...,n-l. 

Since e is small, we can replace < y > — > obtaining the following differential 
equation: 

d 2 C dgdX h fl,! dX m ^ dX m .\ 

Remark. This equation is the geodesic deviation equation of the Lorentz force 
equation. Therefore, the difference between the averaged Jacobi equation and the 
deviation equation associated with the Lorentz force is of higher order in the degree 
(a + b + c). At leading order both equations coincide. 

The condition of transversal dynamics is 

dr dr 

where O 2 indicates first order in the general degree (a + b+c). This is a more general 
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condition than what is usually stated as the transverse dynamics in accelerator 
physics: the magnetic field is perpendicular to the velocity field ^2L of the particles 
in the beam [9]. 

Due to this condition, we suppress the respective term in the differential equations, 
getting the equations 

d 2 e # dx k / 1 dx™ x , / 1 dx m x 

Therefore, at first order, the differential equations are 

dH l & dX> , dXi dX k , dXi dX k dX? di k _ 

In the transverse dynamics, one assumes by construction that ^rF* j = 0, if there 
is not dispersion. Hence, the differential equations in this regime are 

As we mentioned before, the last term of this equation corresponds to the inertial 
term. In a similar way as in the circular motion, the inertial terms is already well 
known [9,10]; the components of the inertial acceleration will be assumed to be 

(VT 2 .rf«,p2 v dX' J dX k dXJde* _ 
(T jk + awr ifc )(__ + 2 — — )-0, 

(T ^+^^( ——+2—— ) = (— ) ^-(^) - = (^) -(-7)" 

In the last expression we consider £ to be small in relation to the curvature radius 
of the central geodesic p. Note that p is not necessarily constant. However, we are 
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assuming a planar trajectory. In the usual formalism, the last term corresponds to 
the relative centripetal force between two particles following close trajectories. 



6.3.2 Examples of transverse linear dynamics 

We study some examples of transverse dynamics using the linearized version of the 
averaged Jacobi equation. 

1. Motion in a normal magnetic dipole 

As we said before, the reference frame is the laboratory frame. The reference 
trajectory is a solution of the averaged Lorentz force equation. In this case 
the electromagnetic field is given by the expression 



where 60 is the dipole strength. Since the magnetic field is constant, ^ l d{F l j = 
0. Therefore, the equations of motion for the transverse degrees of freedom 



/ \ 



F = 



6 



-6 



\ J 



(e 1 ^ 3 ) are 



dr 2 




d 2 e 



dT 2 



=0. 



Changing the parameter of the curve from r 



x , one has that: 



d 2 e 




d 2 e 



= 0. 



dp 



dp 



These are the standard equations for a normal dipole. 



2. Motion in a skew magnetic dipole 
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The electromagnetic field is given by the expression 



( \ 

-bo 

b 

\ J 



In this case, the Jacobi equations are 



d 2 e ( dx 2 i e, 

dr 2 { dr ) p [ p } 



2f3 



o, 



dX 



dr 2 



0. 



Following the same procedure as before we end with the equations for the 
deviation equation in a skew magnetic field: 

d 2 e e d 2 e 

\ = 0, — = o. 



dl 2 



dl 2 



3. Motion in a normal quadrupole field combined with a dipole 

In this case the electromagnetic field has the form 



F(x) 






y o 



\ 

o b - he o 

-&o + b^ 3 

o -he o j 



The Jacobi equation reduces to t 



dr 2 dr 



^fs+(— r-(-M=o, 



e^F 



3,=0. 



dr 2 dr 

Let us consider the respective contributions ^p-^'c^F 1 j and ^^'c^F 3 j. Using 
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Euler's theorem on homogenous functions one gets the relations: 



dXi. 

~dT~ 



Then the differential equations are 



0, ^ + fh = o. 



Using the Euclidean length as a parameter of the curve, one obtains 
d 2 £} i dB 3 rdX^i 1 d?£? ^ 2 ,i ^dB 1 



dl 2 s d^ 1 v dl ' p 2 



In the second differential equation, the second term is zero, since (bj, b^, b|j) = 
(0, 0, b[j). Then we obtain the following differential equations for the transverse 
motion: 



2 C 3 



d 2 t 



dT 2 



+e *i = o. 



These are the equations of the linear transverse dynamics in quadrupoles com- 
bined with magnetic dipole fields using the proper time parameter r. If we 
use the Euclidean length /, the equations are 



dp 



dp 



(6.3.2) 



which are the standard equations in transverse dynamics in accelerator physics 
[9, 10]. 

4. Motion in a normal dipole combined with a 45 degrees quadrupole 

In this case the electromagnetic field is 



F = 








\ 

b + b^ 3 

-b - he o he 

o -he o / 
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Following the same procedure as before, we get the Jacobi equations 



d 2 e 



+ e h + 4 = o, 



d 2 e 



-eb =o. 



(6.3.3) 



dl 2 



dl 2 



The above examples show how the linear transverse dynamics can be obtained from 
the Jacobi equation of the averaged Lorentz connection. 

It is not possible to use this formulation for higher multipole magnetic fields 
because the linear approximation breaks down. One possibility to incorporate higher 
modes is to consider the generalized Jacobi equation [59], which is a non-linear 
geodesic deviation equation. 

We remark again that with the approximation < y > — > the deviation equa- 
tion of the averaged connection and the Lorentz connection coincide. If one considers 
higher order effects, one can obtain differences between the Jacobi equations of the 
Lorentz connections and averaged connection. 

6.4 Calculation of the averaged off-set effect between 
the reference trajectory and the central geodesic of 
the averaged Lorentz connection 

In this section we calculate the averaged difference between the reference trajectory 
and the solutions of the averaged connection. 

6.4.1 Calculation of the dispersion function in beam dynamics 

We follow the formalism developed in [9] for the treatment of linear perturbations 
and dispersion. However we will maintain the proper time r as the parameter of the 
curves, in contrast with the usual treatment, which uses the Euclidean length along 
the reference trajectory. In the following, primes indicate derivatives with respect 
to the proper time. 

It follows from section 6. 3 that the transverse dynamics is determined by equations 
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of the form 

u" +K(t)u = (6.4.1) 
The general solution is of the form 

u(r) = C(r)u + S(t)u' , u'(t) = C'(r)u + S'(r)u' , 

with initial conditions 

C(0) = 1, C"(0) = 0; 5(0) = 0, S'(0) = 1, 

for arbitrary initial values uq and u' . The functions C(r) and S(t) satisfy 

C"(r) + #(t)S(t) = 0, 5"(r) + ^(r)5(r) = 0. 

However, small perturbations can change the dynamics. The perturbed equation 
has the form: 

u"(t) + K{t)u{t) = P (t) (6.4.2) 
A particular solution for (6.4.2) is 

P(t)= [ T p(f)G(T,?)df, (6.4.3) 
J o 

where G(t,t) is the Green function associated to the differential equation (6.4.1). 
One can prove that in the absence of dissipative forces (that is, which do not depend 
on the velocity of the particle) , the Green function of the differential equation is given 
by the following combination: 

G(t, f) = S(t)C(t) - C(t)S(t). (6.4.4) 

Therefore, the general solution for the equation (6.4.2) is 

u(t) = a C(r) + b S(t) + P(r). (6.4.5) 
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This solution breaks down if there are synchrotron radiation or other dissipative 
effects. 

We will use standard notation of beam dynamics. If all the particles in a bunch 
do not have the same energy, one obtains for the transverse degrees of freedom the 
following differential equation [9, pg 109], [10]: 

u" +K{t)u =-(r)A u , A = ^, 5p= V($Pi) 2 + (W + fe) 2 - 
Po Po 

We need to assign a value to dp. One natural value is the maximal value of 
x e M}. Since M is non-compact, we restrict to a compact domain 
K C M. This definition does not depend on the particular trajectory of each parti- 
cle. The general solution for u is linear in the perturbation, therefore 

u(r) = a u C(r)+b u S{T) +AD(t) := a u C(r)+b u S(r) +Off u (r), P(r) = AD(r). 

where a and b depend on the initial values. 

6.4.2 Calculation of the off-set due to the deviation e- 7 

In this subsection (M, 77) is the 4-dimensional Minkowski space-time. We have shown 
in the previous section that in first order of approximation with respect to the 
degree (a + b) of the monomials £ a • (e™) 6 and its derivatives, when we take the 
approximation < y > — > the differential equation for the transverse motion is 
the Jacobi equation of the averaged connection. Therefore we can consider the terms 
on e k in the averaged Jacobi equation as a perturbation and apply the method of 
the Green function. 

From the definition of the off-set function for the transverse degrees of freedom, 
we obtain 

Ofti\r) = u^{r) - a 1)3 C x >\t) - 6 1>3 ^(r), 
the super-index refers to the transverse components x 2 and x 3 in the laboratory 
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frame defined previously. Using the corresponding Green function we obtain 

Offt 3 (r)= [ T p^{T)G(T,f)df. 
Jo 

The perturbation p{r) is in this case defined by all the terms of the averaged Jacobi 
equation which are not contained in the linearized equation respect to the degree 
(a + 6 + c). Therefore let us re- write the Jacobi equation of the averaged connection. 
Using e k =< y k > we get 

Ofti\r) = jT df2 d ^(f) . (F^f) e k (f) + F 1,3 fc (f ) e 3 (f) ) + 

+^^(F 1 ' 3 m (f) ( < y m > (f)r, jk - < y m y a y k > (f) W ) + 



df df 

+ (F 1 ' 3 m (f) ( < y m > (f) Vjk -< y m y a y k > (r) Wfe ))). (6.4.6) 

1 3 

This is an integro-differential equation for Off u ' as we can show. In the integrand 
of equation (6.4.6) we can make the substitution 

de< 3 



df 



-(f) — > (a 1 ,3C" 1 ' 3 (f) + 6 1>3 5' 1 ' 3 (f) + 0// u (f) / ). 



We can also consider the derivatives in the longitudinal and temporal direction using 
this notation, with a convenient choice of the coefficients ao,2 and 60,2 • Then we can 
write 

JL( f ) _^ ( 0j c«{f) + bj S»(f) + Offt{f)'), j = 0, 1, 2, 3. 

In this expression repeated indices are not summed! For the transverse degrees of 
freedom, the unperturbed solutions are the same as before [10], 

u 1 - 3 ^) = u C(t) + u' S(t) C"{t) + K(t)C = 0, S"(t) + K(t)S = 0. 

For the longitudinal j = 2 and temporal j = degrees of freedom, one gets the 
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following relations by comparison with the Jacobi equation, 

^(f) -> (a 2 c\f) + b 2 s'Hf) + off u (fy), 



df 



^(f) — ► (oo C'°(f) + 6 S'°(r) + 0//°(f )') 



Let us consider the regime where 0//° = Offl = 0, Vu. Then the off-set function 
is 

0//^' 3 (t) = ^ df (X;2(a i ^(f)+6 J -^(f)+O/^(f)0(f)~ (-(F 1 ' 3 J (f) efe (f) + 
+F 1 ' 3 fe (f)e j (f)) +— - — - (F 1 ' 3 m (f) ( < y m > (f) Vjk -< y m y a y k > (f) VjaVlk ) + 



df df 

+ {a l C l {f)+b l S l {f)+Off l u {f))d l (F 1,3 m ( < y m > (f) Vjk - < y m y a y l > (f)VjaVlk)))- 

1 3 

This is an integro-differential equation for Off u ' that we formally can solve itera- 
tively. In the Born approximation one puts Off l u (f) = in the integrand: 



ft . j vk . 1 

Off\\r) = j( df (2(^(5)+^^ 

+^^{ Fl ' 3 m(f) {<y m > (f)r, jk - < y m y a y k > (f) W ) + 
+ {a l C\f) + b l S l {f))d l (F 1 - 3 m (f) ( < y m > (f)r, jk - < y m y a y l > (f) Wfc ))). 
Therefore, we get the expression 

Offl'Hr) = jT df (2^ • ^(l(pi.s, (f)cjfc ( f ) +F^ fc (f) Cj -(f)) 

+^^( Fl ' 3 - ( < V m > (f)r, jk -< y m y a y k > (f) Wfc ) + 
+ ?d l (F 1 ' 3 m (f) ( < y m > (f)7ft fc - < y m y a y l > (f) Wfc ))). (6.4.7) 
This expression depends on the particular solution u. A way to eliminate this 



+ 
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dependence is to take the following average 

f T / dX^ /1 

<Offt 3 >(r) = ^ df(2^(f)-^ r (-(F 1 - 3 j (f) efe (f) + F 1 ' 3 fc (f) ej (f)) 

+^^( Fl ' 3 m(f) {<y m > (f)r, jk - < y m y a y k > (f)r ]jaVlk ) + 
+ <C l > (m (F 1 ' 3 m (f) ( < y m > (f) Vjk -< y m y a y l > (f) W i*)))- (6.4.8) 
The averaged off-set is therefore an observable quantity. It is determined by: 

1. The reference trajectory X(t), which is a geodesic of the averaged connection 
and as we have discussed before, it is known theoretically. 

2. The tangent velocity field ^ along the reference trajectory. This is known 
theoretically. 

3. The external electromagnetic field F 1,3 m (x), 

4. The value of the vector field e k (r) = < y k (r) > 

5. The first, second and third moments of the distribution function f(x(r),p{r)) 
along the reference trajectory. 

Finally, in the case that the perturbation does not change significatively along the 
trajectory, we obtain that the term containing derivatives are neglected. Therefore, 

< Of ft 3 > (r) = jf df< (2^(f)— • (F^rMf) + F 1 ' 3 fe (f ) e, (f )) > + 

rl Xi d X^ / \ \ 

+ V V( Fl,3m(f) ( <ym> ^ k - < y mya y k > ( 6 - 4 - 9 ) 

In the case of a delta function distribution we have< Offt 3 > (r) = 0. This 
means that the averaged off-set effect is a collective effect. 
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6.5 Longitudinal beam dynamics and corrections from 
the Jacobi equation of the averaged connection 



Let (M, 77) be the Minkowski space-time and consider an inertial coordinate system 
denned by the vector field Z = J^, that corresponds to the laboratory frame. The 
interaction of an ultra-relativistic bunch of particles with an external longitudinal 
electric field is described by the Faraday tensor 



F = 



( E 2 (x) A 



-E 2 (x) 

/ 



For narrow distributions one obtains the following condition, 



dr dr 



This relation can be seen as follows. For the linear dynamics £ = (£, 0, — £, 0) in 
the laboratory frame. Using the ultra-relativistic limit = (1 + E,0, E, 0), with 



E » 1. 

Then the averaged Jacobi equation for the limit e J 
regime are 



in the ultra-relativistic 



d 2 C n d£, j dX k /l dX m 
+2-? — (-(F l 



dr 2 dr dr 

dXi dX k 



dr 



+ ( T jfc + £ 0, T j7c ) ( , , + 2-7--T-) - 0. 



dr dr v ' dr dT dr dr 

For the above longitudinal electric field, the equations of motion are 

d 2 £° d£, 2 dX k m d^dX 2 ^ 

+ 4 —E 2 <y m > Vmk + -i—^E 2 <y m > Vmk + 



dr 2 dr dr 
'dX k dX 2 



dr dr 
dX 2 dX> 



„i . /dA™ dA" „ dA" dX^ m \ 

H l di[ ^ ^ E 2 <y m > Vmk + ——E 2 <y m > r lmj j + 
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d 2 i 2 dZ dX* rn d£ k dX° 

-J-o- + - A 1—^2 < y > Vmk ~ -1 1—^2 < y > Vmk- 

dr z dr dr dr dr 

m dX Q dX'J m \ 

■^r + Cr ifc + fft'r ^x—-^-^— — )-o, 

rf 2 ^ 2 / nT , 2 ./o nT ,2 x,dXidX k dXidi k , 

In an inertial coordinate system, the inertial terms are zero. Therefore, the system 
of equations in the linear longitudinal dynamics in the ultra-relativistic regime is 

d 2 e° di 2 dx k ^ rn di k dx 2 ^ 

+^ l d l (^ r ^ r E 2 <y™> Vmk ) = 

d 2 i 2 di 2 dX k ^ rn d£ k dX°„ 
dr z dr dr dr dr 

-^ l di(^^E 2 <y m > Vmj ) = 

d 2 X 1 

0, 



dr 2 

<fx 2 

dr 2 

dX k 



0. 



If e =< y > —^f- ~ and since the distribution function has support on the unit 
hyperboloid, < y k > ps 1 + a. Using also the decoupling condition ps 0, 



we have that 



g-f*-tf*<^«)-0, (6.5.2) 
= 0, (6.5.3) 
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rl 2 Y 2 

l*-0. (6.5.4) 
The only non-trivial equation has the form 

-d^ + ^ E2 -^ dl{ ^ E2) -°- 

In the ultra-relativistic limit the velocity field = 7 (in units where the speed of 
light is equal to 1). Therefore, the equation above can be written as 

+ fr^ ~ 2 ^ l diE2{X + £ 2 ) = 0. (6.5.5) 

We perform the following approximation in equation (6.5.5): 

E 2 (X + e) = E 2 {X)+e^:E 2 . 

Due to the translational invariance of the partial derivatives, di = in the above 
expressions, by the chain rule. Then we have 

^ + ^e 2 - 2 1 e^ k mx +o- MX)) = 0. 

If E 2 (X + £) can be approximated linearly on £, using Euler's theorem of homoge- 
neous functions one gets the following expression: 

d 2 £ 2 df 2 

+ ^T^(t) - 2 1 (t)(E 2 ((X + 0- E 2 (X)) = 0. 

6.5.1 Examples 

1. Constant longitudinal electric field. 

In this case the equation of motion is 

dr 2 dr 
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A particular solution is 



£j 2 



2. Alternate longitudinal electric field. In this case, the electric field is 
E 2 (X 2 + £ 2 ) = E 2 (0)sin(w rf (X 2 + £ 2 )). 
The differential equation is 

^+^E 2 (0)sin(w rf {X 2 H 2 ))-2 1 E 2 (0)(sin(w rf (X 2 H 2 ))-sin(w rf X 2 )) = 0. 
We can expand this equation in £, since £ is small 

+ ^£ 2 (0)(sin(to r/ X 2 ) + cos(w rS X 2 )t 2 ) - 2~fE 2 (0)(cos{w rf X 2 )t 2 ) = 0. 



At first order in £ 2 we have the equivalent expression 
d 2 P 2 df 2 

+ ^E 2 (0)sin(w rf X 2 ) - 2 1 E 2 (V)(cos(w rf X 2 )e) = 0. 

We choose the initial phase such that sin(w r fX 2 ) ~ 0; therefore cos(w r fX 2 ) 
1 and the equation is 



^-2 7 £ 2 (0)e 2 = 0. (6.5.6) 

This equation is a linearized version of the ordinary linearized longitudinal 
dynamics and has exactly the same structure [11]. 



6.6 Conclusion 

We have seen that the linear beam dynamics in accelerator physics can be obtained 
from the Jacobi equation of the averaged Lorentz connection for electromagnetic 
fields F(£) linear in £. In particular we have proved that in the case where the 
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magnetic fields are linear on the deviation £, like in a dipole and quadrupole magnetic 
fields, the transverse dynamics can be interpreted as the dynamics of the Jacobi 
equation. A similar conclusion follows for the so-called longitudinal dynamics. 

The advantages of this derivation are that it involves only observable quantities. 
Another advantage is that the averaged dynamics is linked through the distribution 
function to the collective behavior of the system. 

The off-set effect calculated in section 6.4 can be of relevance in the control 
and diagnostics of the beam parameters, since it is a direct observable quantity and 
because it is related with the collective description of the bunch of particles. 
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Chapter 7 



Conclusions 

7.1 General conclusions 

This thesis describes the foundations of the averaged Lorentz force equation and its 
applications in the mathematical modeling of ultra-relativistic bunches of charged 
particles. Through the averaged Lorentz dynamics, a new theoretical justification of 
the use of fluid models in beam dynamics has been obtained. We have seen that for 
relativistic dynamics and for narrow probability distribution functions (in the sense 
that the diameter a of the distribution function obtained using the Euclidean metric 
in the laboratory frame is very small compared with the mass of the particles at rest), 
it is justified to substitute the original kinetic model based on the Vlasov equation 
by an averaged charged cold fluid model. One can control this approximation in 
terms of the energy of the bunch, the diameter of the distribution a and the time of 
the evolution of the bunch, all these variables measured in the laboratory frame. 

Our method does not provide a system of differential equations for the fluid mod- 
els. Instead it provides estimates of the differential operations which occur in the 
definition of the charged cold fluid model. Given a particular physical situation one 
can decide whether the given model is satisfactory or is a bad approximation using 
those bounds. 

The averaged model has been an essential tool in obtaining those results. The 
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reason is that the averaged Lorentz force equation is simpler than the Lorentz force 
equation. The existence of normal coordinate systems associated to the averaged 
connection < L V > has been crucial for the calculations performed in chapter 5. 



There are some advantages using the averaged Lorentz equation in the description 
of the dynamics of a bunch of particles instead of the Lorentz force equation: 

1. At the classical level the electromagnetic field is measured by the effect on 
charged point particles. In case of the electromagnetic interaction of charged 
particles with an external electromagnetic field, the Lorentz force equation is 
the geodesic equation of a complicated non- linear almost- connection (for the 
notion of almost- connection see the appendix). The equation can be simpli- 
fied by considering the associated averaged Lorentz connection. The averaged 
Lorentz connection is simpler than the original one, since the averaged connec- 
tion is an affine connection on M. This property allows us to have important 
technical tools (in particular normal coordinates). 

2. Since the averaged Lorentz connection is simpler than the original Lorentz 
connection, one can use it to perform numerical simulations of the dynamics 
of a bunch of particles. The simplified model must allow a better numerical 
implementation in the simulation of the dynamics of a bunch containing a 
large number of charged particles. 

In a similar way, there are advantages using the averaged Vlasov equation instead 
of the original Vlasov equation: 

1. The calculations using the averaged Vlasov equation can be simplified using 
normal coordinates systems. This is basically because the underlying averaged 
Lorentz connection is an affine and torsion free connection. 

2. It involves only the low moments of the distribution function /. This depen- 
dence on low moments, together with the fact that / is an approximation of / 
in the regime when the dynamics is ultra-relativistic and the distributions are 
narrow, can explain why current fluid models need only to consider differential 
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equations for the first, second and third moments, while the fourth moments 
are neglected. That is, the justification for the truncation schemes in fluid 
models comes from the structure of the averaged Lorentz connection. 

Since the metric approach to geometrization contains intrinsic problems, we have 
considered an alternative geometric treatment of the Lorentz force. The framework 
has been the theory of non-linear connections associated with second order differen- 
tial equations. We have applied this theory to the Lorentz force equation (4.1.1). 

However, starting a dynamical model from a differential equation can be insuffi- 
cient for some purposes. For instance, one could not recover a canonical Hamiltonian 
formalism. From a geometric point of view, if one has only a differential equation, 
one can not speak of variational problems. 

The Lorentz force equation has been interpreted as the non-linear Berwald con- 
nection of a spray vector field L \- However, the above points imply that one has 
to look for a variational formulation. A consistent formulation has been described 
in section 4-3. This formulation is technically complicated and requires sheaves and 
pre-sheaves theory. It is a non-metric approach to semi-Randers spaces. 

7.1.1 Brief discussion of the main problem 

The main problem formulated in section 2.3 has been addressed: we have obtained 
a recipe such that we can decide when the charged cold fluid model is a good 
approximation or not to the underlying kinetic model, in the context of accelerator 
physics. The conclusion is that for the actual accelerator machines, the charged cold 
fluid model is a good approximation to the underlying Vlasov model and that it can 
be used in the description of the beam dynamics. 

As a byproduct we have obtained an averaged Lorentz force equation which is 
simpler than the original Lorentz force equation. We have proved that in physical 
situations, the Lorentz force equation can be substituted by the averaged Lorentz 
equation. 
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7.2 Generalizations and open problems 



Some open problems that the present thesis leaves for future investigation are the 
following: 

1. The geometric averaged method can be applied to other dynamical systems. 
As an example, let us consider the structure of the Lorentz connection L V. If 
viu^v) = 1) tlie- connection coefficients are polynomial in y up to third order. 
Therefore, it could be interesting to consider it in a similar way as an effective 
connection that has semi-spray coefficients of the form: 

(?{x,y) = a 1 jk (x)y j y k +a i jk i{x)y j y k y l +a i jk i m {x)y j y k y l y m +a i jk imn(x)y j y k y l y m y n . 

(7.2.1) 

This is the simplest generalization of the semi-spray coefficients of the connec- 
tion L V. 

There are some restrictions on these semi-sprays coefficients 

(a) The resulting geodesic equations must be gauge invariant. This means 
that there is an intrinsic gauge symmetry transformation and the tensors 
depend only of gauge invariant quantities. 

(b) The resulting equation must be Lorentz invariant. 

(c) The study of the basic dynamics of these connections compatible with 
the laws of the Electrodynamics (in particular with the Larmor law [3, 
pg 469]). A second order dynamics of a charged point particle must be a 
particular class of these dynamics, since it is well known and experimen- 
tally checked. 

(d) It is well known that if the back-reaction force is taken into account, then 
the Lorentz-Dirac equation follows from a balance equation [1-4]. One 
can investigate if a generalization of the type (7.2.1) can accommodate 
a second order equation which remains second order, after considering 
back-reaction. 
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(e) It could be interesting to clarify whose semi-sprays of the type (7.2.1) 
or possible generalizations of the Lorentz force coming from non-linear 
electrodynamics [53]. 

2. The method used in this thesis to justify the use of the charged cold fluid model 
is applicable to other fluid models, like the warm fluid model. Indeed one can 
discriminate which model is better in some particular application, depending 
on the diameter of the given distribution function and the energy of the beam. 

3. One can use the averaged Lorentz equation as a model in numerical simula- 
tions. Since the structure of the averaged equation is simpler than the original 
equation, it could be convenient to use it in numerical simulations in beam 
dynamics. 

4. We have assumed some technical hypotheses. Although these assumptions are 
well defined and hold for the physical examples that we have in mind, it could 
be interesting to reduce the number of assumptions, obtaining more general 
results. 

Apart from these points, directed to the core of this thesis, there are several points 
which could deserve more attention: 

1. The notion of almost- connection. As is explained in the appendix, it is a 
natural generalization of the notion of connection. We think that it is non- 
trivial, since the non-extensibility of the covariant derivatives and parallel 
transport is a difficult property to prove. 

2. The notion of semi-Randers space [14, work in progress]. This a basic notion 
that we have need to discuss but which is of interest on its own. Also it 
is interesting to generalize the notion of semi-Rander space associated with 
non-abelian symmetries. 

3. The notions of structural stability introduced in section 6.4- It indicates a 
topological structure behind current fluid models which deserves investigation. 
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Appendix A 



Mathematical appendix 

A.l Proofs for Chapter 3 

The following proofs are adapted from reference [22]. 

Proof of proposition 3.3.2. The consistency of the equation (3.3.4) is proved 
in the following way, 

< V > X f =< 7T 2 | U V MX) </ > u = 

=< ir 2 \ufhi(X)(ir*f) > u =< TT 2 \u<(X(f)) > u =< Xf > u =Xf. 

The fourth equality holds because the definition of the horizontal local basis { -£p \ u , .... 
in terms of i = 0, ...,n- 1} and {^|-, i = 0, ...,n- 1}. 

We check the properties characterizing a linear covariant derivative associated 
with the averaged connection < V >: 

1. < V >x is a linear application acting on sections of TM: 

< V > x (Yi + Y 2 ) =< V > x il+ <V> X Y 2 



< V > x XY = XV X Y, 
V Y 1 ,Y 2 ,Y g rTM,A G R, X G T^M. 
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(A.l.l) 



For the first equation, the proof consists in the following calculation, 

< V > X (Yi + Y 2 ) =< -K2\iNh u {X)<{Yl + Y 2 ) > = < ^2\uV hu (x)<yi > + 

+ < 7T 2 | U V ' hu(X )K Y 2 >v = < V > X Yi+ < V > x Y 2 . 
For the second condition we have that 

< V > x (AY) =< n 2 \ u V Lu{x) ir* v (\Y) > v = A < n 2 \ u V Lu{X )K Y >= A < V > x Y. 

2. < V >x Y is a J'-linear with respect to X: 

< V > Xl +x 2 Y =< V > Xl Y+ < V >x 2 Y, 

<V> fX (Y) = f(x) < V >x Y, 
MY G TM,t) G 7r _1 (z), X,X u X 2 eT x M, f eT(M). (A.1.2) 
To prove the first equation is enough the following calculation: 

< V >xi+x 2 Y =< k 2 \ u (V hu ( Xl +x 2 ))KY > v = 

= < TT 2 \ U V hu ^ Xl )K Y >v + < TT 2 \ U V hu (x 2 )K Y > = 

= (<V> Xl Y) + (<V>x 2 Y). 

For the second condition the proof is similar. 

3. The Leibnitz rule holds: 

< V > x {fY) = df(X)Y+f <V> X Y, V Y G TTM, / G F(M), X G T,M, 

(A.1.3) 

where df(X) is the action of the 1-form df G A 1 M on X G T^M In order to 
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prove (3.3.7) we use the following property: 

7r*(/y) = Kf<Y, v y G TTM, / G j-(m). 

Then 

< V > x (/y) =< 7r 2 | u V^ (x)7 r:(/y) > u =< Tr 2 \ u Vh u (x)<(fKY > u = 

= < TT 2 \ u {V hu{X) (TT* v f))TT* v {Y) > u + < Tr 2 \ u {7r*J)V hu{x) Tr* v {Y) > u = 
=< 7r 2 | u (/t u (X)(7r*/))7r*(y) > u +f x < tt 2 | u V hu ( X )K( Y ) >«= 

= < (X x f)lT 2 \ u lT* u (Y) > u +f x < 7T 2 \ v V Lu(X) TT* v (Y) > u . 

For the first term we perform the following simplification, 

< (Xf)TT 2 \ u TT* u (Y) > u = (Xf) < TT 2 \ U TT* U (Y) > u = 
= (Xf)(< TT 2 \ U TT* U > U )Y = (Xf)(< I > U )Y. 

Returning to the above calculation, we obtain 

Vx(/y) = Vx(/)y + fV x Y = df(X)Y + fV x Y. 

The generalization of < L V > to higher order tensor bundles is as usual. 
Proof for Proposition 3.3.4 and corollary 3.3.5 

T<v>{X, Y) =< vr 2 | u V feu (x)7r^ > u Y- < tt 2 \ u X7 h u {Y)^w >u X - [X,Y] = 

= < K2\vSh u {X)K >u Y- < TT 2 \ U V hu( Y)K >u X 

-<ir 2 \u<[X,Y] >= 
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=< 7r 2 | u (V MX)7 r*y - V hu(Y) 7T*X - TT*[X,Y}) > u =< T(X,Y) > . 
On the other hand: 

T <V> (X,Y) =< 7r 2 | u V MX) 7r; > u Y- < 7r 2 | u V My) 7r; > u X - [X,Y] = 

=< K2\,Nh u (X)K >u Y- < 7r 2 | u V hji (y)7r* > u X- < ir 2 \uK [ x , Y ] >= 
=< 7r 2 | u (V MX) 7r*y - V hu{Y) 7r*X - tt*[X,Y]) > u = O.D 
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A. 2 Ordinary differential equations 



We state here a result on differential equations that we have used several times in the 
text. It was used in [21] to show the local existence and uniqueness of parameterized 
geodesies and similar results on the existence and uniqueness of solutions of ordinary 
differential equations, 

Theorem A. 2.1 Let fi(t,y,s) be a family of n functions defined in \ t \< 5 and 
(y,s) G D ; where D is an open set in R n+m . If fi(t,y,s) are continuous in t 
and differentiable of class C 1 in y, then there exists a unique family 4>(t,y,s) of n 
functions defined in | t |< 5 and (rj,s) G V, where < 5 < 5 and D is an open 
subset o/D such that 

1. (f>(t,7],s) is differentiable of class C 1 in t andrj. 

2. ^ = f i (t,<t>(t,ri,s),s). 

3. (f)(0,r),s) =r). 

If f(t,y, s) is differentiable of class C p , < p < oo, in y and s, then (p(t,rj,s) is 
differentiable of class C p+1 in t and of class C q in r] and s. 

The main use in the thesis was the following. Assume that the geodesic equations 
are of the form 

% = &{x,v,s); %=v\ i = i,...,n. 

We have applied the theorem to the case where s defines a smooth homotopy. 
Therefore, we have to apply to the case where s is 1-dimensional and / = (y l ,G l ), 
with i = 1, 2n. 

We have used this theorem to prove smoothness properties of the solutions of 
several differential equations. 
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A. 3 Basic notions of asymptotic analysis 



In this appendix we follow the notation of [43, chapter 1]. 

Definition A. 3.1 For two complex functions f(z) and g(z) we write f(z) = 0{g{z)) 
as z — ^ zq iff there is a positive constants K and c such that for < \z — zo\ < c one 
has that \ f\ < K\g\. 

Definition A. 3. 2 An infinite sequence of functions {& n (z), n = 1, 2... } is an asymp- 
totic sequence as z — > zq if 

lim — — =0, Vn G N. 

z^rzo ® n (z) 

Definition A. 3. 3 Given the asymptotic sequence {$ n (z), n = 1,2...} as z — > zq 

the expression 

N-l 

f( z ) = Yj a n®n{z) + 0{<5> N {z)) aSZ -»• Zq, 
n=l 

is said to be an asymptotic expansion of the function f(z) in terms of {Q n (z), n = 
1,2...}. 

Given an asymptotic sequence {& n (z), n = 1,2...}, if for a given function / the 
asymptotic expansion exists it is unique, with coefficients given by 

at=lim (/M-^°W 

Asymptotic expansions have the following elementary properties: 

1. The first order term in the expansion is the leading term, which provides the 
major contribution to the series. 

2. The asymptotic expansion of a function depends on the choice of the asymp- 
totic sequence. 

3. The asymptotic expansion as z — > zq is a linear operation respect to the 
function which is being expanded: / i-> (a 1; ajy^i) is linear in /, ViV > 1. 
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4. If the derivative of the function / has an asymptotic expansion, the asymptotic 
expansion of the derivative of a function is the derivative term by term of the 
expansion of /. 

5. The asymptotic expansion of the real integral of a function is the real integral 
of the expansion, integrated term by term. 

6. The product of asymptotic expansions is in general non-asymptotic. However, 
the product of asymptotic expansions in power series around the same point 
is also an asymptotic expansion. 

If the asymptotic expansion is convergent, there is a convergence region such that 
the approximation of a function by an asymptotic series is becoming more accurate 
when we consider more terms in the expansion. If the asymptotic series is a power 
series and divergent, the better accuracy attainable by an expansion consists on 
taking the expansion 

N 

such that N(z) is the last term that the magnitude of the terms |a n $„| is decreasing 
with n. After this term ajy^N, the rest of the terms are increasing. Usually, the 
error in the expansion in a power series is of the order of the first term neglected. 
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A. 4 Notion of almost-connection 



During our analysis of the Lorentz force equation, we have considered the associated 
non-linear connection on TN and the associated linear connections on 7r*TM — > "E. 
However, strictly speaking the system of differential equations (4.1.1) does not define 
a non-linear connection on TS. The reason is the following. Let us fix a point 
«eS. There exists a natural embedding e : S ^ TM. One can consider the sub- 
bundle e(S) <^-> TM and the covariant derivative L D acting on elements of TTe(S). 
This covariant derivative can be extended to derive sections of the extended bundle 
LLes T « N ' which is a sub-bundle of TM: 

L D : r( | | T U N) x T( [J T U N) — > r( □ T tt N), 

«es «es «es 

From the definition of L D, this operator cannot be extended in a smooth and natural 
way to be a covariant derivative on TM. This is because of the appearance of factors 
y/rj(y, y) in the connection coefficients: the function y/rj(y, y) is not defined in the 
whole TM. The non-extendibility of this function to TM is the origin of the problem 
to extend the covariant derivative along arbitrary directions in TM. 

The above fact suggests the existence of a mathematical object which is a gen- 
eralization of the ordinary notion of covariant derivative in the sense that allows 
covariant derivatives along outer directions to the manifold, although not obtained 
as a restriction of an ambient covariant derivative operator. 

Example. The Lorentz connection was obtained from the Lorentz force equation 
assuming that the torsion is zero. We obtained the following connection coefficients: 

L-pi V~pi _i_ I j i 

1 jk 1 jk t -L jk * ■ L/ jkj 

L V = TTT—^^mk + F^y" 1 ^), 
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V / 1 \ 

T jk = Y l m ( rj jk - — ^jsVkiy'y 1 ) ■ 

This defines a rule to derive sections of TS along directions of TE. The same 
coefficients provide a rule to derive sections of UugE T U TM along directions of 
TM, but we can extend the definition of the covariant derivative acting on sections 
of TTM. 

A related issue is the following. On the unit hyperboloid recall that T l jky^y k = 
0. Therefore, if one considers instead an alternative connection V defined by the 
connection coefficients 

1 jk 1 jk i ^ jki 

the corresponding geodesic equation (parameterized by the proper time of the Lorentzian 
metric rj) is again the Lorentz force equation [49]. Therefore we see that both L V 
and V reproduce the Lorentz force and both are torsion-free (the connection coef- 
ficients are symmetric in the lower indices). This is in contradiction with the fact 
that a connection of Berwald type (linear or non-linear) is determined by the set of 
all geodesies as parameterized curves on the base manifold M and the torsion tensor 
(for linear connections the procedure can be seen in [35]; for non-linear connections, 
a procedure to define the connection is described for example in [34]. 

A solution to this dilemma comes from the fact that the Lorentz force equation 
applies only to time-like trajectories for which the tangent velocity vector fields live 
on the unit hyperboloid XI; that is, we cannot extend the operator 

L D : T( |J TJM) x T( |J TJM) — > r( [J TJM) 

«es «es «es 

to a genuine operator of the form 

L D : r(TM) x r(TM) — ► T(TM). 

In other words, we cannot extract enough information from the Lorentz force equa- 
tion to determine a projective connection [44], because on the base manifold, we 
do not have information about all the possible geodesies. This is why we cannot 
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determine the connection. 

One can be tempted to extend the Lorentz force equation to another equation 
valid for any kind of trajectory. This is partially accomplished by the averaged 
connection. However, these extensions could be non-natural or non unique and 
indeed hide a natural object such like almost- connection. 

Preliminary Definition of Almost- Connect ion 

Let M be a manifold of dimension n. An almost-connection is a spray vector 
field x G TTM defined on a sub-bundle TD ^ TTM, with TD a sub-bundle of 
arbitrary co-dimension. Associated with D is the corresponding almost projective 
covariant derivative on 7r*TM. 

Examples. 

1. A projective connection (in the sense of Cartan [44]) is an almost connection 
such that D = TM. 

2. The Lorentz connection L V provides an example where S = D/ TM. One 
can define a Koszul connection acting on T(|J ugS T U TM). 

3. The averaged covariant derivative < L V > is an affine connection and therefore 
an almost connection in the above sense. 

One can consider the corresponding linear connections L V and V on 7r*TM. 
Generally, their averaged connections < L V > and < V > (introduced in section 
4-1) are not the same. However, we would like to have an averaged operation which 
is well defined for the objects in a given category. By definition this will be the 
category of almost- connections and the corresponding morphisms. We require that 
the result of the averaging operation be the same for each representative belonging 
to the same almost connection, 

Definition A. 4.1 Let M be a manifold of dimension n. A non-linear almost con- 
nection is the maximal set of semi-sprays \ defined on a sub-bundle TD C TTM 
such that they have the same averaged linear covariant derivative < V > and the 
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same torsion tensor 

f(X, Y) := L DyY - L D t X - [X, Y], X, Y G TTD. (A.4.1) 

Associated with \ over D is the corresponding almost-covariant derivative on the 
bundle 7r*TM — > D. Any of the connections in the same almost-connection has 
the same averaged connection. The distance function (4.6.3) can also being defined 
for almost-connections. 

The general properties of almost connections are being explored in a separate 
work. Some of these properties are based on straightforward generalizations of the 
quantities associated with Koszul connections. For instance, the generalization of 
the curvature tensor is 

R(X, Y, Z) := L D jt L D t Z - L D f L D ji Z - L D { ^ f] Z, X,Y,Z G TTD. (A.4.2) 

It is not easy to handle a notion of parallel transport for almost-connections, since 
in general there will be initial conditions such that the auto-parallel curve goes out 
from the sub-bundle D, even for arbitrary short-time parallel transports. Indeed, for 
auto-parallel curves whose initial velocity vector is not on TD, the trajectory goes 
out from D after any finite time. This point is related with the notion of general 
connection [62]. However, the notion of almost-connection is even more general, 
since the solution of the projections on N condition V/j(£)7r*(£) = could not exists. 
For a generalized connection, this projection always exist. 
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A. 5 Basic notions of Sobolev spaces 



Sobolev spaces are complete normed vector spaces (therefore Banach spaces) where 
the norms measure also the derivatives of the function. We have used Sobolev norms 
in chapter 5 to introduce the bounds on some differential expressions appearing in 
the averaged Vlasov model. In this appendix we provide the basic notions of Sobolev 
norms and some additional notions of Sobolev spaces to understand the meaning of 
these expressions and its implications for further generalizations. We will follow 
references [41] and [57] because of their clarity in exposition. We assume that the 
theory of Lebesgue's integral holds. 

Let fl be an open set of a manifold. There are some basic definitions: 

Definition A. 5.1 Two functions f,g:Q — > RU {±00} are equivalent iff they are 
equal almost everywhere on fl, that is, the sub-set A C fl where g 7^ / is a null set. 



Definition A. 5. 2 Let Q be open, p > 1 (p e R). L P (Q) is the set of all Lebesgue 
measurable functions f : fl — > RU {±00} for which \ f\ p is integrable over fi. For 
f G LP(Sl) we set 



In order to define Sobolev norms, we introduce weak derivatives. 

Definition A. 5. 3 Let f G L 1 (ft). A function v G L 1 (J7) is called the weak deriva- 
tive of f in the direction x 1 if 



Definition A.5.4 Let f G L 1 (f2) ; p := (ft, ...,&) with ft > (i = l,...,d), |ft := 
EtiA>0. Then 



Definition A. 5. 5 A function v G L 1 (J1) is called (3-weak derivative of f and written 





d(j)(x) 
dx 1 



dx. 




(11). 
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v = Dpf if 

[ v{x)<j>{x) dx = / f{x)Dp(f)(x)dx. 

We can now define Sobolev spaces and Sobolev norms: 

Definition A. 5. 6 For k G N a natural number, 1 < p < oo, we define the Sobolev 
space W k > p (n) byW k > p (n) := {/ G D>(Q) \ D p f exists and is in U>{n)}. The Sobolev 
norms are defined by 

Wf\\K P --={ E / W(x)\dx)K (A.5.1) 
\/3\<k n 

and by 

\\fh, P ■= ^2 maxsup xe n\Dpf(x)\. (A.5.2) 

\P\<k 

A.5.1 Basic properties of Sobolev spaces 

Some properties of Sobolev spaces are the following: 

Theorem A. 5. 7 Let k G N be a natural number and 1 < p < oo. Then the 
following is true 

1. The normed space (W k ' p (Q), ||/||fc, p ) is a Banach space. 

2. Let f G W k ' p (fl) and 9 G C 1 (R). Then Poincare's inequality holds: 

\f\ {0tP) <C(Sl,p)\f\ hp , f G (W c ^(n)), (A.5.3) 

with C(fl,p) a constant and Wc' p (fi) the completion of the space of smooth 
functions on CI with compact support. 

3. WH ,p {fl) is a Hilbert space. 

There is a relevant result (Sobolev embedding theorem) which gives sufficient con- 
ditions to embed Sobolev spaces in L q spaces, in spaces of continuous functions or 
in spaces with some regularity conditions [41,57,58]. We do not need this theorem 
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here, but it could be important for of further generalizations of the results of chapter 
5. 



A. 5. 2 Sobolev spaces of functions defined on manifolds 

The discussion before of Sobolev spaces is restricted to open domains. One can 
extend the definition to manifolds. First one reduces to an open domain U C M. If 
the manifold is differentiable, there are partitions of the unity [32]. First, one can 
speak of Wc' p (U). Using an atlas of the manifold and a associated partition of the 
unity, we can define the Sobolev norms (W fc ' p (M), ||/||fc, p ) from (W fc ' p (fi), ||/||fc, p ) in 
the usual way as the integrals are defined over manifolds from the local description 
using coordinates neighborhoods [32, chapter 4}- 

A. 5. 3 Examples of Sobolev spaces 

There are some Sobolev norms that we have used in chapter 5. These are: 

1. (W 1 ' 1 (M), ||/||i,i). This is a Hilbert space, with a norm defined by the function 



II/IIwm :=( E / W(x)\dx) = [ (\f(x)\ +\Y J d l f\)dx 
i£rr, J si Jsi ; 



l/?l<i 



2. (W°' 2 (M), ||/||o,2)- The norm of a function is defined as 
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